Predicted mutation
evidence seq id position mutation annotation gene description
RA CP028101 693,938 C→T pseudogene (890/1244 nt) C7Y58_03400 ← IS3 family transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP028101693,9380CT100.0% 101.0 / NA 42pseudogene (890/1244 nt)C7Y58_03400IS3 family transposase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/13);  total (29/13)

ATCTAAACTATTTCCTTTTCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGCTTTTCTCCTCTTAAATTAAATTCTGTTACATCTGTAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAATATGATTA  >  CP028101/693695‑694025
                                                                                                                                                                                                                                                   |                                                                                       
aTCTAAACTATTTCCTTTTCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCTAGTTATGACTGCGCGAAATATCATACGAAACTATAGATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTtct                                                                                   >  1:702695/1‑250 (MQ=255)
aTCTAAACTATTTCCTTTTCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTtct                                                                                   >  1:364502/1‑250 (MQ=25)
aTCTAAACTATTTCCTTTTCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTtct                                                                                   >  2:529935/1‑250 (MQ=25)
       cTATTTCCTTTTCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCttaa                                                                           >  1:307056/1‑251 (MQ=21)
       cTATTTCCTTTTCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCttaa                                                                           >  2:631964/1‑251 (MQ=21)
         aTTTCCTTTTCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCGTAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCttaaat                                                                         >  1:750961/1‑251 (MQ=21)
              cTTTTCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTtaaattaaat                                                                    >  1:869405/1‑251 (MQ=21)
                tttCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTc                                                                  <  2:648304/251‑1 (MQ=21)
                tttCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATATTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTc                                                                  <  1:494909/251‑1 (MQ=21)
                           aCTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTc                                                                  <  1:458400/240‑1 (MQ=25)
                           aCTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTc                                                                  >  2:458400/1‑240 (MQ=25)
                                gAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATc                                                         <  1:69225/244‑1 (MQ=25)
                                gAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATc                                                         >  2:69225/1‑244 (MQ=25)
                                 aGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAAc                                                  <  1:73741/250‑1 (MQ=25)
                                  gTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCa                                                >  1:130298/1‑251 (MQ=21)
                                          cttctCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGAt                                        <  2:123400/251‑1 (MQ=21)
                                              tctTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTgg                                    >  1:500974/1‑251 (MQ=21)
                                               ctTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGa                                   >  2:144479/1‑251 (MQ=21)
                                                 ttCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGc                                 <  2:364502/251‑1 (MQ=21)
                                                     gTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGc                                 >  1:376658/1‑247 (MQ=18)
                                                     gTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGc                                 <  2:376658/247‑1 (MQ=18)
                                                       ctttcttGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATAAATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTg                            >  2:907001/1‑250 (MQ=18)
                                                         ttcttGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCtt                         >  2:372891/1‑251 (MQ=18)
                                                              gATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCaaaa                    >  2:749543/1‑251 (MQ=18)
                                                              gATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCaaaa                    >  1:224380/1‑251 (MQ=18)
                                                              gATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCAAAATATCATACGAAACTATATATCTGCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCaaaa                    >  1:755467/1‑251 (MQ=16)
                                                                   gAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAAtct                 >  2:566973/1‑249 (MQ=17)
                                                                   gAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAAtct                 >  1:202005/1‑249 (MQ=17)
                                                                   gAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAAtct                 >  2:639825/1‑249 (MQ=17)
                                                                   gAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAAtct                 >  1:583429/1‑249 (MQ=17)
                                                                     ataatGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAAtctct               >  1:535992/1‑249 (MQ=17)
                                                                        atGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTaa          >  2:48025/1‑251 (MQ=14)
                                                                         tGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGTACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAat         >  2:596995/1‑251 (MQ=14)
                                                                         tGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAat         >  2:553352/1‑251 (MQ=14)
                                                                         tGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAat         >  2:601584/1‑251 (MQ=14)
                                                                         tGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAat         >  1:6185/1‑251 (MQ=14)
                                                                         tGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTCGCTTCAAAATCTCTCTTAat         >  1:622708/1‑251 (MQ=14)
                                                                           cTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAata        <  1:536439/250‑1 (MQ=14)
                                                                           cTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAata        <  2:755467/250‑1 (MQ=14)
                                                                           cTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCACTTTCTTTAAATGCTAAATTTATCATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAata        <  1:639825/250‑1 (MQ=9)
                                                                              aTACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAATATGAt    <  2:202005/251‑1 (MQ=14)
                                                                                aCTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAATATGATTa  <  1:453302/251‑1 (MQ=14)
                                                                                                                                                                                                                                                   |                                                                                       
ATCTAAACTATTTCCTTTTCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGCTTTTCTCCTCTTAAATTAAATTCTGTTACATCTGTAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAATATGATTA  >  CP028101/693695‑694025

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑