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.

GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑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)
                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                >  2:552468/1‑100 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                               >  2:2459558/1‑100 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTCATCCGGCCTACaa                                                                               <  1:2714990/100‑1 (MQ=255)
                           gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                              >  1:3112902/1‑100 (MQ=255)
                           gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                              >  1:1625780/1‑100 (MQ=255)
                             tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                            <  1:2884159/100‑1 (MQ=255)
                             tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTAGGCCGGATAAGGCGTTAACGCCGCATCCGGcataaacaaa                                                                             <  2:4276910/100‑10 (MQ=21)
                              gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                            <  2:4015599/99‑1 (MQ=255)
                              gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                           >  1:437807/1‑100 (MQ=255)
                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                          <  1:410782/100‑1 (MQ=255)
                                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                       <  1:3776775/100‑1 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                <  2:2999329/100‑1 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                >  2:1553579/1‑100 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                >  2:2124495/1‑100 (MQ=255)
                                           caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                              <  2:2402383/100‑1 (MQ=255)
                                            aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                             <  1:1556991/100‑1 (MQ=255)
                                             gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                            <  1:1349932/100‑1 (MQ=255)
                                              cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACa                                                           <  2:1835161/100‑1 (MQ=255)
                                                ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                         >  2:3730608/1‑100 (MQ=255)
                                                      ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaa                                                   >  1:1493618/1‑100 (MQ=255)
                                                      ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAAATTGCAGaa                                                   >  2:3682449/1‑100 (MQ=255)
                                                          taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACCAAATCGTGCAAATTCAATACATTGCAGAAttt                                                >  2:2488582/1‑99 (MQ=255)
                                                             aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑Gt                                       >  1:2812120/1‑100 (MQ=255)
                                                                 tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGc                                   <  2:3929540/100‑1 (MQ=255)
                                                                 tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGc                                   >  2:2514801/1‑100 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGcc                                  <  1:867488/100‑1 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGcc                                  >  2:2644443/1‑100 (MQ=255)
                                                                   ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGcc                                  <  1:2013523/99‑1 (MQ=255)
                                                                    tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTg                                >  1:2002146/1‑100 (MQ=255)
                                                                      aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGAt                              <  1:259877/100‑1 (MQ=255)
                                                                       cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATa                             <  2:1672650/100‑1 (MQ=255)
                                                                            tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCgg                        >  1:4209601/1‑100 (MQ=255)
                                                                              ctcTTTGTTTATGCCAGATGCGGCTTCAACGTCTTATCCTGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAg                      <  1:407485/98‑1 (MQ=17)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATATGCGGAgc                     <  2:2314723/100‑1 (MQ=255)
                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                    >  1:532550/1‑100 (MQ=255)
                                                                                cTTTGTTTATGCCGAATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCTCGCAAATTCAATACATTGCCGAATTT‑‑‑TC‑‑GTAGGGCCGATAAGccg                        >  1:3558499/1‑94 (MQ=17)
                                                                                  ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCa                  <  1:2258535/100‑1 (MQ=255)
                                                                                   tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTAT‑‑‑TG‑‑GTAGGCCTGATAAGCGCAGCGCAt                 >  2:449426/8‑100 (MQ=21)
                                                                                   tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATc                >  2:1251666/8‑100 (MQ=21)
                                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATc                >  2:1471074/1‑100 (MQ=255)
                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               <  2:977857/100‑1 (MQ=255)
                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               <  2:3407641/100‑1 (MQ=255)
                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               >  2:1127174/1‑100 (MQ=255)
                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               <  1:2281558/100‑1 (MQ=255)
                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               <  1:754853/100‑1 (MQ=255)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCA‑‑‑T‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  1:3657391/5‑100 (MQ=18)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  2:1053409/5‑100 (MQ=17)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAAC‑‑‑AC‑‑GTAGGCCTGATAAGCGAAGCGCATCAg              >  2:4147909/5‑100 (MQ=18)
                                                                                      gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCATACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              >  1:1720797/5‑100 (MQ=18)
                                                                                       taggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAAATCGTGCTAATTCAATATATTGCAGAAACC‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:4184301/4‑100 (MQ=17)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTCC‑‑‑GT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:1048740/4‑100 (MQ=17)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAACTTTGCAAATTCAATATATTGC‑‑AATCT‑‑‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:3139000/4‑100 (MQ=18)
                                                                                       tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg             >  2:93083/1‑100 (MQ=255)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATCC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:2767858/3‑100 (MQ=21)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:1335584/3‑100 (MQ=21)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:1998899/3‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCAACCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:2673004/3‑100 (MQ=17)
                                                                                         tGCCGGATGCGGCGTTAACTACATATCCGGCCTACAAAATCGAGCAAATTCAATACATTGCAGAACTT‑‑‑TC‑‑GTACGCCTGATAAGCGAAGCGCAACAGGc            >  1:3194132/1‑99 (MQ=17)
                                                                                         tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  2:519332/1‑100 (MQ=255)
                                                                                         tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  2:3325062/1‑100 (MQ=255)
                                                                                         tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  1:1295317/1‑100 (MQ=255)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCA‑‑‑T‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:3323310/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTAT‑‑‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa           >  1:2019798/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTAT‑‑‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa           >  2:2243457/2‑100 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCC‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:2437896/2‑100 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCC‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:753226/2‑100 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCT‑‑‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:3464451/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCT‑‑‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:63589/2‑100 (MQ=17)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑AATCT‑‑‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:2769529/2‑100 (MQ=18)
                                                                                          gtcGGATAAGGCGTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACG‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCCTCAGGcca          >  1:2146003/3‑98 (MQ=21)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACC‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:3923365/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGC‑‑AATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa           >  2:3042605/1‑100 (MQ=21)
                                                                                          gCCGGATAAGGCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAATTTCAATATATTGCAGAGATT‑‑‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa          >  1:3056236/1‑99 (MQ=14)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTCC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  2:4198592/1‑99 (MQ=21)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTAT‑‑‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa          >  2:4355208/1‑99 (MQ=25)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCC‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa          >  1:1085725/1‑99 (MQ=21)
                                                                                            cggctacggcgTTAACGCAGCATCCGGACTACACAATCGTGCAAATTCAATACATTGCAGACTTT‑‑‑TC‑‑GTCGGCCTGATAAGCGGAGCGCATCAGGCAttt        >  1:1805084/7‑100 (MQ=17)
                                                                                             ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt       <  1:3231395/100‑1 (MQ=255)
                                                                                                  cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttttttt  >  1:685449/1‑96 (MQ=255)
                                                                                                  cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  >  1:3437667/1‑100 (MQ=255)
                                                                                                  |                                                                                                         
GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑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: