Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,971 C→A 29.5% intergenic (‑181/+127) 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,9710CA29.5% 553.0 / 251.1 276intergenic (‑181/+127)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (97/97);  new base A (76/5);  total (174/102)
Fisher's exact test for biased strand distribution p-value = 1.71e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-02

Only 100 of 293 total aligned reads displayed.

GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070
                                                                                                  |                                                                                                         
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                          >  2:3560538/1‑100 (MQ=255)
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                          >  2:1421419/1‑100 (MQ=255)
 ngAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGGCAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                         >  1:606291/2‑100 (MQ=255)
   aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                       <  2:895158/100‑1 (MQ=255)
     tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                     >  1:3499661/1‑100 (MQ=255)
     tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                     <  2:148511/100‑1 (MQ=255)
      gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                    <  2:271862/100‑1 (MQ=255)
      gAAGGATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                    <  1:1359184/100‑1 (MQ=255)
       aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                   <  2:1540494/100‑1 (MQ=255)
        aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                  >  2:3772074/1‑100 (MQ=255)
        aGTATTGATGCCGGCTGGCGCTGCACGCAGAGGCACAGCGCCCCCGGGTATACAGATTGTTGACAAGGATAGCTTTGTTTATACGGACTGCGGCGTTaa                                                                                                   >  1:1876403/1‑99 (MQ=255)
         gTATTGCTGCCGGATGGCGCTGCACGTAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                 <  1:1825075/100‑1 (MQ=255)
         gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                 <  2:3619328/100‑1 (MQ=255)
           aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                               <  2:1947613/100‑1 (MQ=255)
              gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                            >  1:1662104/1‑100 (MQ=255)
               ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                           >  2:2558601/1‑100 (MQ=255)
                 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                         >  2:2894380/1‑100 (MQ=255)
                    ggATGGCGCTGTACGCAGAGGCACAGCGCCACCGTGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                      <  1:1011624/100‑1 (MQ=25)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                     <  2:1316653/100‑2 (MQ=37)
                      aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                     <  2:4163479/99‑1 (MQ=255)
                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                >  2:552468/1‑100 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                               >  2:2459558/1‑100 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTCATCCGGCCTAC‑aa                                                                               <  1:2714990/100‑1 (MQ=255)
                           gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                              >  1:3112902/1‑100 (MQ=255)
                           gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                              >  1:1625780/1‑100 (MQ=255)
                             tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                            <  1:2884159/100‑1 (MQ=255)
                             tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTAGGCCGGATAAGGCGTTAACGCCGCATCCGGcataaacaaa                                                                             <  2:4276910/100‑10 (MQ=21)
                              gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                            <  2:4015599/99‑1 (MQ=255)
                              gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                           >  1:437807/1‑100 (MQ=255)
                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                          <  1:410782/100‑1 (MQ=255)
                                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                       <  1:3776775/100‑1 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                <  2:2999329/100‑1 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                >  2:1553579/1‑100 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                >  2:2124495/1‑100 (MQ=255)
                                           caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                              <  2:2402383/100‑1 (MQ=255)
                                            aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATa                                                             <  1:1556991/100‑1 (MQ=255)
                                             gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                            <  1:1349932/100‑1 (MQ=255)
                                              cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACa                                                           <  2:1835161/100‑1 (MQ=255)
                                                ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                         >  2:3730608/1‑100 (MQ=255)
                                                      ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGaa                                                   >  1:1493618/1‑100 (MQ=255)
                                                      ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAAATTGCAGaa                                                   >  2:3682449/1‑100 (MQ=255)
                                                          taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑CAAATCGTGCAAATTCAATACATTGCAGAAttt                                                >  2:2488582/1‑99 (MQ=255)
                                                             aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑Gt                                       >  1:2812120/1‑100 (MQ=255)
                                                                 tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGc                                   <  2:3929540/100‑1 (MQ=255)
                                                                 tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGc                                   >  2:2514801/1‑100 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGcc                                  <  1:867488/100‑1 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGcc                                  >  2:2644443/1‑100 (MQ=255)
                                                                   ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGcc                                  <  1:2013523/99‑1 (MQ=255)
                                                                    tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTg                                >  1:2002146/1‑100 (MQ=255)
                                                                      aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGAt                              <  1:259877/100‑1 (MQ=255)
                                                                       cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATa                             <  2:1672650/100‑1 (MQ=255)
                                                                            tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCgg                        >  1:4209601/1‑100 (MQ=255)
                                                                              ctcTTTGTTTATGCCAGATGCGGCTTCAACGTCTTATCCTGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAg                      <  1:407485/98‑1 (MQ=17)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATATGCGGAgc                     <  2:2314723/100‑1 (MQ=255)
                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                    >  1:532550/1‑100 (MQ=255)
                                                                                cTTTGTTTATGCCGAATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCTCGCAAATTCAATACATTGCCGAATTT‑‑‑TC‑‑GTAGGGCCGATAAGccg                        >  1:3558499/1‑94 (MQ=17)
                                                                                  ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCa                  <  1:2258535/100‑1 (MQ=255)
                                                                                   tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAGTTTAT‑‑‑TG‑‑GTAGGCCTGATAAGCGCAGCGCAt                 >  2:449426/8‑100 (MQ=21)
                                                                                   tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCAGGAGTT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATc                >  2:1251666/8‑100 (MQ=21)
                                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATc                >  2:1471074/1‑100 (MQ=255)
                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               <  2:977857/100‑1 (MQ=255)
                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               <  2:3407641/100‑1 (MQ=255)
                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               >  2:1127174/1‑100 (MQ=255)
                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               <  1:2281558/100‑1 (MQ=255)
                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               <  1:754853/100‑1 (MQ=255)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCA‑‑‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  1:3657391/5‑100 (MQ=18)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCAAGATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  2:1053409/5‑100 (MQ=17)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCAGGAAAC‑‑‑AC‑‑GTAGGCCTGATAAGCGAAGCGCATCAg              >  2:4147909/5‑100 (MQ=18)
                                                                                      gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCATAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              >  1:1720797/5‑100 (MQ=18)
                                                                                       taggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAATCGTGCTAATTCAATATATTGCAGAAACC‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:4184301/4‑100 (MQ=17)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGCAATCTCC‑‑‑GT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:1048740/4‑100 (MQ=17)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAAACTTTGCAAATTCAATATATTGC‑‑AATCT‑‑‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:3139000/4‑100 (MQ=18)
                                                                                       tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg             >  2:93083/1‑100 (MQ=255)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGATCGATCC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:2767858/3‑100 (MQ=21)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCAGGAGTT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:1335584/3‑100 (MQ=21)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCAGGAGTT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:1998899/3‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCAACCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCAGAGATT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:2673004/3‑100 (MQ=17)
                                                                                         tGCCGGATGCGGCGTTAACTACATATCCGGCCTAC‑AAAATCGAGCAAATTCAATACATTGCAGAACTT‑‑‑TC‑‑GTACGCCTGATAAGCGAAGCGCAACAGGc            >  1:3194132/1‑99 (MQ=17)
                                                                                         tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  2:519332/1‑100 (MQ=255)
                                                                                         tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  2:3325062/1‑100 (MQ=255)
                                                                                         tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  1:1295317/1‑100 (MQ=255)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCA‑‑‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:3323310/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAGTTTAT‑‑‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa           >  1:2019798/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAGTTTAT‑‑‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa           >  2:2243457/2‑100 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCACCATCC‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:2437896/2‑100 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCACCATCC‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:753226/2‑100 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCAGGAGCT‑‑‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:3464451/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCAGGAGCT‑‑‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:63589/2‑100 (MQ=17)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑AATCT‑‑‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:2769529/2‑100 (MQ=18)
                                                                                          gtcGGATAAGGCGTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCAGGAACG‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCCTCAGGcca          >  1:2146003/3‑98 (MQ=21)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAGAAACC‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:3923365/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑AATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa           >  2:3042605/1‑100 (MQ=21)
                                                                                          gCCGGATAAGGCGTTCACTCCGCATCCGGCCTAC‑AAAAGCACGCAATTTCAATATATTGCAGAGATT‑‑‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa          >  1:3056236/1‑99 (MQ=14)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTGGGCTCC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  2:4198592/1‑99 (MQ=21)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAGTTTAT‑‑‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa          >  2:4355208/1‑99 (MQ=25)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCACCATCC‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa          >  1:1085725/1‑99 (MQ=21)
                                                                                            cggctacggcgTTAACGCAGCATCCGGACTAC‑ACAATCGTGCAAATTCAATACATTGCAGACTTT‑‑‑TC‑‑GTCGGCCTGATAAGCGGAGCGCATCAGGCAttt        >  1:1805084/7‑100 (MQ=17)
                                                                                             ggATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt       <  1:3231395/100‑1 (MQ=255)
                                                                                                  cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttttttt  >  1:685449/1‑96 (MQ=255)
                                                                                                  cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  >  1:3437667/1‑100 (MQ=255)
                                                                                                  |                                                                                                         
GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070

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