Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 29.9% intergenic (‑194/+114) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9840TG29.9% 352.6 / 186.2 188intergenic (‑194/+114)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (69/63);  new base G (49/7);  total (118/70)
Fisher's exact test for biased strand distribution p-value = 2.67e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-02

Only 100 of 192 total aligned reads displayed.

TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTC  >  NC_000913/3705885‑3706082
                                                                                                   |                                                                                                           
 tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                            >  2:281038/1‑100 (MQ=255)
  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                           >  2:2739913/1‑100 (MQ=255)
  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                           >  2:2684643/1‑100 (MQ=255)
     gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                        <  1:2782064/100‑1 (MQ=255)
       cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                      <  1:2865785/100‑1 (MQ=255)
       cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                      <  1:984871/100‑1 (MQ=255)
        ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                     <  2:230048/100‑1 (MQ=37)
        ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTATGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                     <  1:1018245/100‑1 (MQ=37)
          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                   <  2:3297967/100‑3 (MQ=37)
          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                   <  2:602555/100‑3 (MQ=37)
              cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑aa                                                                                            >  2:3311698/1‑100 (MQ=255)
                  gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATc                                                                                        >  1:3209056/1‑100 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCg                                                                                       <  2:276078/100‑1 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTg                                                                                     <  2:3148290/100‑1 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGc                                                                                    <  1:2910478/100‑1 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGc                                                                                    <  1:2990643/100‑1 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGc                                                                                    <  1:3030007/100‑1 (MQ=255)
                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCa                                                                                   >  2:3006712/1‑100 (MQ=255)
                        agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCaa                                                                                  <  2:2808960/100‑1 (MQ=255)
                         gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCaa                                                                                  <  2:2148281/99‑1 (MQ=255)
                         gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCAGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCaaa                                                                                 >  2:512624/1‑100 (MQ=255)
                           ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAAtt                                                                               >  2:2543383/1‑100 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAAt                                                                           <  1:1981949/100‑1 (MQ=255)
                                aGCGCCACCGTGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGTCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATa                                                                          <  2:362836/100‑1 (MQ=255)
                                 gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATAc                                                                         <  1:406756/100‑1 (MQ=255)
                                 gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATAc                                                                         <  1:3554152/100‑1 (MQ=255)
                                    ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACAtt                                                                      <  2:1750981/100‑1 (MQ=255)
                                      aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGc                                                                    <  2:406812/100‑1 (MQ=255)
                                            tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑t                                                           >  1:1918168/1‑100 (MQ=255)
                                             ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑tt                                                          <  2:1780143/100‑1 (MQ=255)
                                                 aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑Gt                                                   <  2:2162275/100‑1 (MQ=255)
                                                  gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTa                                                  >  1:2973093/1‑100 (MQ=255)
                                                    ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAgg                                                >  1:494795/1‑100 (MQ=255)
                                                       ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCt                                             <  1:2013829/100‑1 (MQ=255)
                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTg                                            >  1:75724/1‑100 (MQ=255)
                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTg                                            <  1:2940/100‑1 (MQ=255)
                                                           cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCTTTCTCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATa                                         <  1:2492857/100‑1 (MQ=255)
                                                              aGTGCGCTTTGGTTATGCCGTATGCGGCGTTAACGCCTTAACCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAg                                       <  1:2134437/99‑1 (MQ=255)
                                                               gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCg                                     >  2:3193447/1‑100 (MQ=255)
                                                                tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCgg                                    <  2:2713378/100‑1 (MQ=255)
                                                                  ctcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAg                                  <  1:533548/98‑1 (MQ=255)
                                                                  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAg                                  >  1:599879/1‑100 (MQ=255)
                                                                   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAgc                                 <  1:2422677/100‑1 (MQ=255)
                                                                    tgatttgtaggCCGGATAACGCATTCACGCCGCATCCGGCCT‑‑AC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑TT‑AT‑‑GTAGTCCTGATAAGCGTAgcgc                               >  1:1612746/11‑100 (MQ=18)
                                                                    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                >  2:72301/1‑100 (MQ=255)
                                                                        ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CC‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATc                            >  2:1442898/7‑100 (MQ=18)
                                                                        ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCAT‑‑AC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGAT‑TA‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATc                            >  2:3109539/7‑100 (MQ=17)
                                                                        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATc                            <  1:1136920/100‑1 (MQ=255)
                                                                         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                           >  2:2466740/1‑100 (MQ=255)
                                                                         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                           <  2:1421305/100‑1 (MQ=255)
                                                                         ngtaggCCGGATAAGGCGTTCACGCCGCATCCGGCAT‑‑AC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AA‑AC‑AC‑‑GTAGGCCTGATAAGCGAAGCGCATCa                           >  2:163321/6‑100 (MQ=18)
                                                                           taggCCGGATAAGGCGTTTACGTCGCATCCGGCCT‑‑AC‑AAAAGAGTGCAAATTCAATACATTGCA‑TG‑GG‑CC‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                         >  1:1260320/4‑100 (MQ=17)
                                                                           taggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAACATATTGC‑‑‑A‑AT‑TC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAgg                         >  2:2022583/4‑100 (MQ=18)
                                                                           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                         <  1:788065/100‑1 (MQ=255)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCC‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  1:2724005/3‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑A‑TA‑TG‑AC‑‑GTAGGCCTGATAAGCGTAGCGCaagatcgg                       >  2:2526873/3‑93 (MQ=18)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGCA‑AT‑CT‑CC‑GT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                        >  2:955032/3‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCAT‑‑AC‑AAAATTGTGTAAATTCAATAAATTGAT‑CG‑AT‑CC‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                        >  2:1950722/3‑100 (MQ=21)
                                                                             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                       >  1:3083606/1‑100 (MQ=255)
                                                                             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATCCAATACATCGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGACAAGCGGAGCGCATCAGGCa                       >  1:694511/1‑100 (MQ=255)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑CC‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  2:623107/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CC‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  2:3343404/2‑100 (MQ=21)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑AT‑CT‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  2:1269527/2‑100 (MQ=18)
                                                                              gtcggataaggcgTTCACGCCGCATCCGACCC‑‑AC‑AAAGTCTTGCAAAATCAATAAATTGCAGGA‑AC‑GA‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                      >  2:152731/10‑99 (MQ=21)
                                                                              gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAACTTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                        >  2:3389846/1‑98 (MQ=255)
                                                                              gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                      >  1:1659209/1‑100 (MQ=255)
                                                                              gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                      >  1:1510454/1‑100 (MQ=255)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCC‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                     >  2:2264169/1‑100 (MQ=18)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCC‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                     >  1:1808540/1‑100 (MQ=18)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CC‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                      >  2:1643232/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TC‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                      >  2:2846143/1‑99 (MQ=17)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CA‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                      >  1:2365248/1‑99 (MQ=18)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CA‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                     >  1:334901/9‑98 (MQ=18)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacg                     >  2:1411698/9‑98 (MQ=17)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑AT‑CT‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                     <  2:2101170/92‑3 (MQ=18)
                                                                                cggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CC‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                    >  1:666677/8‑100 (MQ=18)
                                                                                cGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                    <  2:1825694/100‑1 (MQ=255)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑AT‑CCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                  >  2:1348226/7‑98 (MQ=21)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CC‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                   >  2:437911/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CC‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                   >  2:870255/7‑96 (MQ=21)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CC‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                   >  1:2585068/7‑96 (MQ=21)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AA‑AC‑AC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                   >  1:1997370/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑TT‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                  >  1:327839/7‑98 (MQ=21)
                                                                                  gataaggcgTTTACGCCGCATCCGGCATAAAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CC‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                    >  1:3251720/6‑100 (MQ=14)
                                                                                  gataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTGTGCAAATTCAATCAATTGCA‑GG‑AA‑AC‑AC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                   >  2:849191/6‑99 (MQ=17)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATTCATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                 <  2:1380178/100‑1 (MQ=255)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                 <  1:1101177/100‑1 (MQ=255)
                                                                                      cggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct              >  2:1019939/1‑100 (MQ=255)
                                                                                       ggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt             >  2:2890662/1‑100 (MQ=255)
                                                                                        gcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGTCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt            <  2:2507712/100‑1 (MQ=255)
                                                                                          tctAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc          <  1:1201389/98‑1 (MQ=255)
                                                                                          gTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc          >  2:1228273/1‑100 (MQ=255)
                                                                                            tAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa        >  1:2317624/1‑100 (MQ=255)
                                                                                             aaCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcaga       >  2:1201396/1‑97 (MQ=255)
                                                                                             aaCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt       >  1:1349859/1‑100 (MQ=255)
                                                                                              aCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc      >  2:2956892/1‑100 (MQ=255)
                                                                                              aCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc      >  2:940539/1‑100 (MQ=255)
                                                                                               cGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc      <  2:3421840/99‑1 (MQ=255)
                                                                                                gCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg    >  1:1948386/1‑100 (MQ=255)
                                                                                                gCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg    >  2:2216038/1‑100 (MQ=255)
                                                                                                   |                                                                                                           
TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTC  >  NC_000913/3705885‑3706082

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑