Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 30.3% 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,9840TG30.3% 579.0 / 306.1 313intergenic (‑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 (109/109);  new base G (86/9);  total (195/118)
Fisher's exact test for biased strand distribution p-value = 7.56e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-02

Only 100 of 315 total aligned reads displayed.

TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTC  >  NC_000913/3705886‑3706082
                                                                                                  |                                                                                                           
 gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                           >  1:2252936/1‑100 (MQ=255)
  ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                          >  1:3246108/1‑100 (MQ=255)
   tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCTTTAACGCCTTacc                                                                                                         <  2:2898226/100‑3 (MQ=255)
   tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                         >  2:2335483/1‑100 (MQ=255)
    gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGAGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                        >  1:3039823/1‑100 (MQ=255)
      cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                      >  2:2086148/1‑100 (MQ=255)
      cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                      <  1:189310/100‑1 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                    <  2:3968923/100‑2 (MQ=37)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                    <  1:3732519/100‑2 (MQ=37)
         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                   <  2:4047473/100‑1 (MQ=255)
         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                   <  2:646497/100‑3 (MQ=37)
            gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑a                                                                                             <  2:3414368/100‑1 (MQ=255)
             cgttGCACGCAGAGGTACAGCGCCCCCGGGTATATAGATTGTTGTCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑a                                                                                             <  1:3858042/96‑1 (MQ=21)
             cgcTGTTCGCAGAGGCACAGCGCTACCGGGTATACAGATTGGTGACAAAGTGCGCTTTGTTTATTCCGGATAAGGCGTTATCGCCGCATCCGGCAT‑‑aaac                                                                                             <  1:1002361/100‑4 (MQ=21)
              gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑aaa                                                                                           >  2:1877786/1‑100 (MQ=255)
              gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑aaa                                                                                           >  1:1268005/1‑100 (MQ=255)
              gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑aaa                                                                                           >  2:1919390/1‑100 (MQ=255)
               cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑aaaa                                                                                          <  2:3086585/100‑1 (MQ=255)
                tGCAGGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAAt                                                                                         <  1:2319343/100‑1 (MQ=255)
                tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAAt                                                                                         <  2:1124123/100‑1 (MQ=255)
                tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAAt                                                                                         <  2:1767500/100‑1 (MQ=255)
                   aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGt                                                                                      >  1:4011915/1‑100 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTg                                                                                     <  1:2980139/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTg                                                                                     <  1:147500/100‑1 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGc                                                                                    <  2:1994244/100‑1 (MQ=255)
                          ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAAtt                                                                               >  1:2217753/1‑100 (MQ=255)
                          ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAAtt                                                                               >  2:2262042/1‑100 (MQ=255)
                          gacacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAAtt                                                                               >  2:3654339/3‑100 (MQ=255)
                           gtgcaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGTTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAAGCGTGCa                                                                                   <  1:3399583/92‑1 (MQ=255)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCa                                                                             >  2:4241280/1‑100 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAAt                                                                           <  2:1035893/100‑1 (MQ=255)
                               aGCGCCACCGGGTCTACAGATTGTTGACAAAGTGCGCCTTGTTTATACCAGATTCGGCGTTCACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTTAATa                                                                          >  1:462295/1‑100 (MQ=17)
                                     aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑c                                                                   <  2:2502954/100‑1 (MQ=255)
                                             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑At‑ttt                                                        <  2:3308850/100‑1 (MQ=255)
                                              aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTc                                                       >  1:577399/1‑100 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑g                                                    <  2:3196924/100‑1 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑g                                                    <  2:2109352/100‑1 (MQ=255)
                                                aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGAGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑At‑ttt                                                        >  1:407070/1‑97 (MQ=255)
                                                   gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAg                                                 >  2:2848951/7‑99 (MQ=17)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGcc                                              <  1:644276/100‑1 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGcc                                              <  1:2871173/100‑1 (MQ=255)
                                                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAggagat                                            >  2:1686246/1‑96 (MQ=255)
                                                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTg                                            <  2:1370459/99‑1 (MQ=255)
                                                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGa                                           >  1:206378/1‑100 (MQ=255)
                                                         aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGAt                                          <  1:2038855/100‑1 (MQ=255)
                                                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATaa                                        >  1:3012156/1‑98 (MQ=255)
                                                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCg                                     >  2:2175300/1‑100 (MQ=255)
                                                               tGCGCGTTGTTTATGCCGGATTCGGAGTTATCGGCTTATCTGGCCT‑‑AC‑AAAATCGTGCAAACTCAAGACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                                    <  1:2624705/100‑1 (MQ=18)
                                                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAg                                  <  2:3472903/100‑1 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                 >  2:4162666/1‑100 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                 <  1:1567754/100‑1 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                >  2:547377/1‑100 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                <  1:1531116/100‑1 (MQ=255)
                                                                    nctttgtaggCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAC‑AAAGTCCTGCAGATTCAATATATTG‑CA‑GG‑TA‑CTGC‑‑GTAGGCCTGATAAGCGTAgc                                 >  1:3130058/10‑100 (MQ=9)
                                                                    attttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCATAAAC‑AAAATCGTGCAAATTCAATAAATTG‑CA‑GA‑AT‑TTCC‑‑GTAGGCATGCTAAGCGGa                                   >  1:3543556/10‑98 (MQ=2)
                                                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                              <  1:3805391/100‑1 (MQ=255)
                                                                      tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAACATATTG‑CA‑AT‑TC‑TCTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                             >  2:27686/8‑100 (MQ=18)
                                                                      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                             >  1:2553322/1‑100 (MQ=255)
                                                                      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                             >  1:1600225/1‑100 (MQ=255)
                                                                        tgtaggCCGGATAAGGCGCTTACGCCGCATCCGGCCT‑‑AC‑AAGAACGTGCAAATTCAACAAATTGTGATGC‑AT‑CCT‑‑‑GTAGGCCTGATAAGCGTAGCGCATc                            >  2:115490/6‑100 (MQ=21)
                                                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                          >  1:1781256/1‑100 (MQ=255)
                                                                         gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTG‑CA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                          >  2:2557763/5‑100 (MQ=17)
                                                                         gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCAT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                         >  2:679164/5‑100 (MQ=18)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                         >  2:251099/4‑100 (MQ=18)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTG‑CA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                         >  2:3502878/4‑100 (MQ=21)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑‑AAATCATGCAAATCCAAAATATTG‑CA‑GG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                        >  1:3457958/4‑100 (MQ=21)
                                                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAGTCGTGCAAATTCACTACATTG‑AG‑GA‑AT‑TTAC‑‑GTAGGCCTGATAAGGGGAGCGCATCAgg                         >  1:861178/1‑100 (MQ=17)
                                                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                         >  2:184129/1‑100 (MQ=255)
                                                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                         <  1:1429360/100‑1 (MQ=255)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                        >  2:3253923/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                        >  2:384718/3‑100 (MQ=18)
                                                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                        <  1:1191859/100‑1 (MQ=255)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GA‑AA‑CCTT‑‑CTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  1:2564641/2‑100 (MQ=17)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTG‑C‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  2:2640927/2‑100 (MQ=18)
                                                                             gtcggataaggcgTTCACGCCGCATCCGACCC‑‑AC‑AAAGTCTTGCAAAATCAATAAATTG‑CA‑GG‑AA‑CGAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCaa                      >  1:3884023/10‑99 (MQ=21)
                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                      >  2:4083402/1‑100 (MQ=255)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                     >  1:1646114/1‑100 (MQ=18)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                      >  2:2704278/1‑99 (MQ=18)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                  >  2:2957040/7‑98 (MQ=21)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTG‑CA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                   >  1:3491711/7‑100 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                   >  1:1127275/7‑100 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTG‑CA‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                   >  2:1312661/7‑96 (MQ=25)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTG‑CA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                   >  1:280916/7‑96 (MQ=21)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAAACATGCAAATTCAATAGATTG‑CA‑GAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                   >  2:494128/7‑100 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAAACATGCAAATTCAATAGATTG‑CA‑GA‑GA‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                   >  1:1320148/7‑100 (MQ=17)
                                                                                ggataaggcgTTAACGCCGCACCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                   >  2:3754633/7‑100 (MQ=18)
                                                                                ggATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                   <  1:1960470/100‑1 (MQ=255)
                                                                                  aTGCGGCGTTAACGTCTTATCCGGCCT‑‑AC‑AAAATCGGGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                 <  1:2836547/100‑1 (MQ=255)
                                                                                  aTGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                 >  1:1841852/1‑100 (MQ=255)
                                                                                   taaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAGATCGTGCAAATTCAATATATTG‑CA‑TG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAAGCAGTtttgc                >  2:1971434/4‑100 (MQ=17)
                                                                                   tGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                >  2:2425261/1‑100 (MQ=255)
                                                                                   tGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                <  1:813497/100‑1 (MQ=255)
                                                                                     cggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATACGCGGAGCGCATCAGGCATTtttgcct              <  2:2605978/100‑1 (MQ=255)
                                                                                        cgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg           <  1:2889977/100‑1 (MQ=255)
                                                                                         gTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc          >  2:771379/1‑100 (MQ=255)
                                                                                          ttAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc          <  2:871169/99‑1 (MQ=255)
                                                                                            accGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt       >  2:2888554/3‑100 (MQ=255)
                                                                                              cGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa     >  2:1487354/1‑100 (MQ=255)
                                                                                               gCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg    >  2:1515926/1‑100 (MQ=255)
                                                                                                 cTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc  <  2:1605103/100‑1 (MQ=255)
                                                                                                  |                                                                                                           
TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTC  >  NC_000913/3705886‑3706082

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