Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 23.3% 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,0480GT23.3% 595.5 / 164.9 246intergenic (‑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 (70/114);  new base T (49/6);  total (126/120)
Fisher's exact test for biased strand distribution p-value = 5.62e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.88e-01

Only 100 of 249 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑3706146
                                                                                                      |                                                                                                  
  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAg                                                                                                  <  1:3028648/100‑1 (MQ=255)
  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAg                                                                                                  <  1:2768602/100‑1 (MQ=255)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAgcg                                                                                                >  2:2536623/1‑100 (MQ=255)
      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCa                                                                                              <  1:2319709‑M2/100‑2 (MQ=255)
        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCatc                                                                                            <  1:1243170‑M2/100‑4 (MQ=255)
         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCatca                                                                                           >  1:2578396‑M2/1‑96 (MQ=255)
           taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCACGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                         >  2:141498/4‑100 (MQ=17)
           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTTATGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  2:2064445/4‑100 (MQ=21)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAgg                                                                                         >  1:1066213/3‑99 (MQ=21)
            aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATGCGCGGCGCTTATCAGGc                                                                                        <  1:2180039/100‑1 (MQ=21)
             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                       >  1:1574772/1‑100 (MQ=255)
             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                       >  1:1677161/1‑100 (MQ=255)
             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                       >  1:346974/1‑100 (MQ=255)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                                                      >  1:1546718/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                                                      >  1:1544029/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                                                      >  1:758786/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:628195/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:2496663/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCa                                                                                       >  1:2719297/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                       >  2:3142817/2‑100 (MQ=18)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTCCACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                      >  1:2426375/1‑99 (MQ=18)
               ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTGCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                                                     <  1:73055/100‑1 (MQ=255)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  2:3245867/7‑100 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  1:2054670/7‑100 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  2:125044/7‑100 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                                                   >  2:2154284/7‑100 (MQ=14)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                                                   >  2:2741940/7‑100 (MQ=17)
                  gataaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTCCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                                                  <  1:478085/95‑1 (MQ=21)
                   ataaggcgTTTACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAatttta                                                                                 >  2:852679/5‑99 (MQ=18)
                   aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  2:1203661/100‑1 (MQ=255)
                   aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  1:493807/100‑1 (MQ=255)
                   aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  1:2862721/100‑1 (MQ=255)
                   aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  1:3156448/100‑1 (MQ=255)
                    tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                                                <  2:592659/100‑1 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                              >  1:1487863/1‑100 (MQ=255)
                       ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                                                             >  2:513034/1‑100 (MQ=255)
                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                           <  2:1574763/100‑1 (MQ=255)
                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt                                                                            >  1:1136216/1‑99 (MQ=255)
                          gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                                          <  2:1807535/100‑1 (MQ=255)
                            tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                                        >  2:1455542/1‑100 (MQ=255)
                               cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                                                     >  2:3089330/1‑100 (MQ=255)
                                gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                                                    <  1:2458775/100‑1 (MQ=255)
                                     aTCCGGCCTATAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                               <  1:42761/100‑1 (MQ=255)
                                        cGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                                                            <  2:1188646/100‑1 (MQ=255)
                                            cTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                                                        >  1:1621840/1‑100 (MQ=255)
                                              aCAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                      >  1:1397043/1‑100 (MQ=255)
                                                 aaaaTCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                                                    <  2:1651247/100‑1 (MQ=255)
                                                        tGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  1:407171/100‑1 (MQ=255)
                                                        tGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCGTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  2:3036559/100‑1 (MQ=255)
                                                         gCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTagatcggaa                                            >  2:2420668/1‑92 (MQ=255)
                                                         gCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                            >  1:3237323/1‑100 (MQ=255)
                                                           aaaTTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                                           <  1:877622/99‑1 (MQ=255)
                                                           aaaTTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                                          <  2:107808/100‑1 (MQ=255)
                                                            aaTTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa                                         <  2:877360/100‑1 (MQ=255)
                                                             aTTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag                                        <  1:2536641/100‑1 (MQ=255)
                                                               tCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag                                      <  2:771543/100‑1 (MQ=255)
                                                                cAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga                                     <  1:898533/100‑1 (MQ=255)
                                                                cAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga                                     <  1:234274/100‑1 (MQ=255)
                                                                 aaTACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCTATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGAc                                    >  2:1836600/1‑100 (MQ=255)
                                                                         gCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            <  2:2947838/100‑1 (MQ=255)
                                                                          cAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa                           >  2:407211/1‑100 (MQ=255)
                                                                          cAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa                           >  1:845539/1‑100 (MQ=255)
                                                                             gAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                         <  2:2027922/100‑1 (MQ=255)
                                                                              ccgatctcgtaggCCTGATAAGCGTAGCGCATCAGGCctacgaaaatcttgcaatgtattgaattttcatgattttgtaggccggataaggcgttcacgc                         <  2:1513200‑M1/94‑64 (MQ=255)
                                                                              aaTTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg                        <  2:222898/100‑1 (MQ=255)
                                                                              aaTTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg                        <  1:178610/100‑1 (MQ=255)
                                                                                ttttCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAt                      <  1:1658622/100‑1 (MQ=255)
                                                                                  tttCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATTGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          <  2:1879458/95‑1 (MQ=255)
                                                                                    tCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGtt                   <  1:1917657/100‑1 (MQ=255)
                                                                                     cGTAGGCCTGATAAGCGGAGCGCatcaggcctacgaaaattctgcaatgtattgaatttgcacgattttgtaggccggataaggcgttaacgccgcatcc                  <  1:1738871‑M2/100‑78 (MQ=255)
                                                                                      gTAGGCCTGATAAGCGTAGCGCATCAGGCctacgaaaatcttgcaatgtattgaattttcatgattttgtaggccggataaggcgttcacgccgcatccg                 >  2:1654213‑M1/1‑29 (MQ=255)
                                                                                      gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                 >  1:667050/1‑100 (MQ=255)
                                                                                      gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                 <  1:1087684/100‑1 (MQ=255)
                                                                                       tAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGTCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                 <  2:1287815/99‑1 (MQ=255)
                                                                                        aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  1:96943/100‑1 (MQ=255)
                                                                                        aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  2:1691153/100‑1 (MQ=255)
                                                                                        aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  2:1544183/100‑1 (MQ=255)
                                                                                        aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  1:1112885/100‑1 (MQ=255)
                                                                                         ttccgatctaGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt              <  2:803631/91‑1 (MQ=255)
                                                                                         ggCCTGATAAGCGTAGCGCATCAGGCctacatgatttctgcaatatattgaatttggaagaatttgtaggccggataaggcgtttacgccgcatctggca              >  2:1489341‑M1/1‑26 (MQ=255)
                                                                                         ggCCTGATAAGCGTAGCGCATCAGGCctacatgatttctgcaatatattgaatttggaagaatttgtaggccggataaggcgtttacgccgcatctggca              >  2:322014‑M1/1‑26 (MQ=255)
                                                                                         ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTATTGTATGGCTTTTAATTATTCGATAGAGACGATTTCGAAGTGATGTTTGCCt              >  2:543353/1‑100 (MQ=255)
                                                                                          gCCTGATAAGCGTAGCGCatcaggcctacaaggtttctgcaatatattgaattagcacgattttgtaggccggataaggcgtttacgctgcatccggcat             <  1:594502‑M2/100‑83 (MQ=255)
                                                                                          gCCTGATAAGCGTAGCGCatcaggcctacaaggtttctgcaatatattgaattagcacgattttgtaggccggataaggcgtttacgccgcatccggcat             <  2:90364‑M2/100‑83 (MQ=255)
                                                                                          gCCTGATAAGCGGAGCGCatcaggcctacgaaaattctgcaatgtattgaatttgcacgattttgtaggccggataaggcgttaacgccgcatccggcat             <  2:96290‑M2/100‑83 (MQ=255)
                                                                                          gCCTGATAAGCGGAGCGCatcaggcctacgaaaattctgcaatgtattgaatttgcacgattttgtaggccggataaggcgttaacgccgcatccggcat             <  1:1186655‑M2/100‑83 (MQ=255)
                                                                                          gCCTGATAAGCGGAGCGCatcaggcctacgaaaattctgcaatgtattgaatttgcacgattttgtaggccggataaggcgttaacgccgcatccggcat             <  1:1992686‑M2/100‑83 (MQ=255)
                                                                                           ccTGATAAGCGTAGCGCatcaggcctacaaggtttctgcaatatattgaattagcacgattttgtaggccggataaggcgtttacgccgcatccggcata            <  1:2119614‑M2/100‑84 (MQ=255)
                                                                                           ccTGATAAGCGTAGCGCATCAGGCctacatgatttctgcaatatattgaatttggaagaatttgtaggccggataaggcgtttacgccgcatctggcata            >  2:173377‑M1/1‑24 (MQ=255)
                                                                                           ccTGATAAGCGTAGCGCATCAGGCctacatgatttctgcaatatattgaatttggaagaatttgtaggccggataaggcgtttacgccgcatctggcata            >  2:1954015‑M1/1‑24 (MQ=255)
                                                                                           ccTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTaa            >  1:2673882/1‑100 (MQ=255)
                                                                                           ccTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTaa            >  1:2979237/1‑100 (MQ=255)
                                                                                           ccTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTaa            >  1:3180543/1‑100 (MQ=255)
                                                                                            cTGATAAGCGGAGCGCATCAGGCATTTTTTTTTTTTTTTTTTTTCTTGTATGGCTGTTAATTATTTGATAGAGACGTTTTGGAAGTGATGTTTGCCttct           >  2:3153354/1‑97 (MQ=255)
                                                                                             tGATAAGCGGAGCGTATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg          <  2:1605472/100‑1 (MQ=255)
                                                                                                tAAGCGTAGCGCATCAGGCctacatgatttctgcaatatattgaatttggaagaatttgtaggccggataaggcgtttacgccgcatctggcataaacaa       >  2:360169‑M1/1‑19 (MQ=255)
                                                                                                  aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATTGATCa     >  2:2968442/1‑100 (MQ=255)
                                                                                                  aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa     >  1:1277038/1‑100 (MQ=255)
                                                                                                    cGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAc   >  1:227621/1‑100 (MQ=255)
                                                                                                     ggAGCGCATCAGGCATTTTTTTTTTTTTTTTCATTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc  >  2:1361042/1‑100 (MQ=255)
                                                                                                      |                                                                                                  
TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑3706146

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