Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 27.0% intergenic (‑258/+50) 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,706,0480GT27.0% 573.4 / 186.3 263intergenic (‑258/+50)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (90/97);  new base T (64/3);  total (163/100)
Fisher's exact test for biased strand distribution p-value = 2.11e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

Only 100 of 266 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑3706146
                                                                                                         |                                                                                                  
 gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGa                                                                                                   >  1:2010450/1‑100 (MQ=255)
  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAg                                                                                                  <  1:2709829/100‑1 (MQ=255)
   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                 >  2:2063265/1‑100 (MQ=255)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                <  1:1419714/100‑1 (MQ=255)
      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                                              >  2:3356334/1‑100 (MQ=255)
       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                                                                                             >  1:1198358/1‑100 (MQ=255)
        ttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATc                                                                                            >  2:2307935/7‑100 (MQ=18)
        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                                                            <  2:1256428/100‑1 (MQ=255)
         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                                                                                           >  2:2341635/1‑100 (MQ=255)
           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  1:2069677/4‑100 (MQ=18)
            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                         >  2:1063249/3‑99 (MQ=17)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGc                                                                                        >  2:172828/3‑100 (MQ=18)
            aTGCCGGATGCGGCGTTAACGCCGTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                                         <  1:1470/99‑1 (MQ=255)
             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                       >  1:1945101/1‑100 (MQ=255)
             ggtcggataaggcgTTCACGCCGCATCCGGCCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:2767944/11‑100 (MQ=18)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  1:386876/2‑100 (MQ=14)
              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                      >  1:2885851/1‑100 (MQ=255)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTCCAC‑‑GTAGGCCTGATAAGCGTAGCTCATCAGGCag                                                                                      >  1:2497197/1‑99 (MQ=18)
               ccggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                     >  2:1005942/9‑98 (MQ=18)
               ccggataaggcgTTCACGCCGCACCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTC‑‑‑‑AGGCCTGATAAGCGTAGCGCATCAGGCttt                                                                                     >  1:131964/9‑99 (MQ=25)
                cggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑CTAGGCCTGATAAGCGTAGCGCATCAGGCAtttg                                                                                   >  1:3469626/8‑99 (MQ=18)
                cggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                                                    >  1:1803314/8‑97 (MQ=17)
                 ggataaggcgTTTACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                                                   >  2:148673/7‑100 (MQ=18)
                 ggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                                                  >  2:1381571/7‑98 (MQ=18)
                 ggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                                                  >  1:1032296/7‑98 (MQ=18)
                 ggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                                                  >  1:1381035/7‑98 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                                                   >  1:368198/7‑96 (MQ=25)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  2:557841/7‑100 (MQ=21)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                                                                                   >  2:2193868/7‑96 (MQ=21)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                                                                                   >  2:1211845/7‑96 (MQ=21)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  1:2177461/7‑100 (MQ=17)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCCGGGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                                                   >  1:3378774/7‑96 (MQ=18)
                 ggataaggcgTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCTGAGCGCATCAGGCAtt                                                                                     >  1:3122850/7‑98 (MQ=17)
                  gataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagttct                                                                                  >  1:2310443/6‑95 (MQ=18)
                   ataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttg                                                                                 <  2:2145583/96‑1 (MQ=17)
                   aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  2:2232917/100‑1 (MQ=255)
                   aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  1:3025232/100‑1 (MQ=255)
                   aTGCGGCGTTAACGCCTTATCCGGACTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  1:3610125/100‑1 (MQ=255)
                    tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                                                <  2:1952107/100‑1 (MQ=255)
                       ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                                                             >  2:3596102/1‑100 (MQ=255)
                         cgTTTACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAACATATTGCCACTCACCCAGTAGGCCTGATAAGCGCAGCGCATCAGGCAattttacattt                                                                            >  2:2714125/1‑94 (MQ=25)
                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                           >  2:3673206/1‑100 (MQ=255)
                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                           <  2:2958997/100‑1 (MQ=255)
                           ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcn                                                                         <  1:281638/100‑2 (MQ=255)
                            tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                                        <  1:874478/100‑1 (MQ=255)
                            tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                                        >  2:1786395/1‑100 (MQ=255)
                             aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                                                       >  1:1838723/1‑100 (MQ=255)
                              aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                      >  2:1623522/1‑100 (MQ=255)
                              aCGCCTTATCCGGACTACAAAATCGTGAAAATTAAATACATTGCAGCATTTTA‑‑GTAGGAGTGATAAGAGGAGCGCATCAGGCATTtttgcctttggtatc                                                                      >  1:2364465/1‑95 (MQ=18)
                               cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGACTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                                                     >  1:3549058/1‑100 (MQ=255)
                                 ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                                                   >  1:837584/1‑100 (MQ=255)
                                   ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                                                 <  2:597113/100‑1 (MQ=255)
                                     aTCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAtacaagact                                                               <  2:469193/100‑9 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGTGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                             <  2:2692026/100‑1 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                             >  2:896758/1‑100 (MQ=255)
                                           ccTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                         <  1:2799469/100‑1 (MQ=255)
                                            cTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                                                        >  1:3082533/1‑100 (MQ=255)
                                             tACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                       <  1:2642341/100‑1 (MQ=255)
                                             tACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                       >  1:3267130/1‑100 (MQ=255)
                                                 aaaaTCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                                                    <  1:853209/100‑1 (MQ=255)
                                                  aaaTCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                                                   <  1:1230150/100‑1 (MQ=255)
                                                    aTCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatt                                                 <  2:721508/100‑1 (MQ=255)
                                                      cGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                                               >  1:221091/1‑100 (MQ=255)
                                                        tGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  2:1109801/100‑1 (MQ=255)
                                                         gCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                            >  1:27021/1‑100 (MQ=255)
                                                         gCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                            <  2:2477566/100‑1 (MQ=255)
                                                            aaTTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa                                         >  2:2385970/1‑100 (MQ=255)
                                                             aTTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag                                        <  2:3628064/100‑1 (MQ=255)
                                                               tCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag                                      <  1:1971069/100‑1 (MQ=255)
                                                                cAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga                                     >  1:2422744/1‑100 (MQ=255)
                                                                 aaTACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGAc                                    >  2:2521314/1‑100 (MQ=255)
                                                                 aaTACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGAc                                    >  2:402391/1‑100 (MQ=255)
                                                                     cATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt                                >  2:3212731/1‑100 (MQ=255)
                                                                         gCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            <  2:30052/100‑1 (MQ=255)
                                                                         gCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            >  2:2886699/1‑100 (MQ=255)
                                                                          cAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa                           >  2:19145/1‑100 (MQ=255)
                                                                           aGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGTGACGTTTTCGAAg                          <  2:2588363/100‑1 (MQ=255)
                                                                             gAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTCTGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          >  2:2673547/1‑99 (MQ=255)
                                                                             gAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                         >  2:801042/1‑100 (MQ=255)
                                                                              aaTTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                         >  1:529011/1‑99 (MQ=255)
                                                                              aaTTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg                        <  2:2797521/100‑1 (MQ=255)
                                                                                aTTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa                       >  2:1528450/1‑100 (MQ=255)
                                                                                 ttttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGATTCGTTTTCGAAGTGAt                      >  1:1577528/1‑100 (MQ=255)
                                                                                     tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGtt                   <  1:580638/100‑1 (MQ=255)
                                                                                         gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGTCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                 <  2:949173/100‑1 (MQ=255)
                                                                                         gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                 <  2:13740/100‑1 (MQ=255)
                                                                                           aGGCCTGATAAGCGGAGCGCATCAGGTATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  1:1905825/100‑1 (MQ=255)
                                                                                            ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt              >  2:3262512/1‑100 (MQ=255)
                                                                                            ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt              >  1:2767700/1‑100 (MQ=255)
                                                                                             gCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTa             <  1:1528441/100‑1 (MQ=255)
                                                                                              ccTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTaa            >  2:3016328/1‑100 (MQ=255)
                                                                                                   tAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt       >  1:3754554/1‑100 (MQ=255)
                                                                                                   tAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt       <  1:2759358/100‑1 (MQ=255)
                                                                                                   tAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt       <  1:107219/100‑1 (MQ=255)
                                                                                                   tAAGCGGAGCGCACCAGGTATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt       <  1:445300/100‑1 (MQ=255)
                                                                                                    aaGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATTAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATc      >  2:3277672/1‑100 (MQ=255)
                                                                                                     aGCGGAGCGCATCAGGCATTTTTTTTTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa     >  1:2588637/1‑100 (MQ=255)
                                                                                                     aGCGGAGCGCATCAGGCATTTTTTCTTTTTTTTTCCGTCTAGTATGGCTTATAATTATTCAATAGAAAAGTTTTCGTAGTGATGTTTGCCTCTTGGATCa     >  1:799969/1‑100 (MQ=255)
                                                                                                      gCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCaa    >  2:1717800/1‑100 (MQ=255)
                                                                                                      gCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCaa    >  2:1448122/1‑100 (MQ=255)
                                                                                                         |                                                                                                  
TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑3706146

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: