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.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑3706146
                                                                                                      |                                                                                                  
  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAg                                                                                                  <  1:3028648/100‑1 (MQ=255)
  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAg                                                                                                  <  1:2768602/100‑1 (MQ=255)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAgcg                                                                                                >  2:2536623/1‑100 (MQ=255)
      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCa                                                                                              <  1:2319709‑M2/100‑2 (MQ=255)
        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCatc                                                                                            <  1:1243170‑M2/100‑4 (MQ=255)
         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCatca                                                                                           >  1:2578396‑M2/1‑96 (MQ=255)
           taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                         >  2:141498/4‑100 (MQ=17)
           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  2:2064445/4‑100 (MQ=21)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAgg                                                                                         >  1:1066213/3‑99 (MQ=21)
            aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑ATTGGTAGGCCTGATGCGCGGCGCTTATCAGGc                                                                                        <  1:2180039/100‑1 (MQ=21)
             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                       >  1:1574772/1‑100 (MQ=255)
             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                       >  1:1677161/1‑100 (MQ=255)
             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                       >  1:346974/1‑100 (MQ=255)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                                                      >  1:1546718/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                                                      >  1:1544029/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                                                      >  1:758786/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:628195/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:2496663/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GGAA‑ACACGTAGGCCTGATAAGCGAAGCGCATCAGGCa                                                                                       >  1:2719297/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                       >  2:3142817/2‑100 (MQ=18)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GGCT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                      >  1:2426375/1‑99 (MQ=18)
               ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTGCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                                                     <  1:73055/100‑1 (MQ=255)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  2:3245867/7‑100 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GGAG‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  1:2054670/7‑100 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GGAG‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  2:125044/7‑100 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGA‑TCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                                                   >  2:2154284/7‑100 (MQ=14)
                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGA‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                                                   >  2:2741940/7‑100 (MQ=17)
                  gataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GGCT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                                                  <  1:478085/95‑1 (MQ=21)
                   ataaggcgTTTACGCCTTATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATGGGC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttta                                                                                 >  2:852679/5‑99 (MQ=18)
                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  2:1203661/100‑1 (MQ=255)
                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  1:493807/100‑1 (MQ=255)
                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  1:2862721/100‑1 (MQ=255)
                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  1:3156448/100‑1 (MQ=255)
                    tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                                                <  2:592659/100‑1 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                              >  1:1487863/1‑100 (MQ=255)
                       ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                                                             >  2:513034/1‑100 (MQ=255)
                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                           <  2:1574763/100‑1 (MQ=255)
                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt                                                                            >  1:1136216/1‑99 (MQ=255)
                          gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                                          <  2:1807535/100‑1 (MQ=255)
                            tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                                        >  2:1455542/1‑100 (MQ=255)
                               cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                                                     >  2:3089330/1‑100 (MQ=255)
                                gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                                                    <  1:2458775/100‑1 (MQ=255)
                                     aTCCGGCCTAT‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                               <  1:42761/100‑1 (MQ=255)
                                        cGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                                                            <  2:1188646/100‑1 (MQ=255)
                                            cTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                                                        >  1:1621840/1‑100 (MQ=255)
                                              aC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                      >  1:1397043/1‑100 (MQ=255)
                                                 aaaaTCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                                                    <  2:1651247/100‑1 (MQ=255)
                                                        tGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  1:407171/100‑1 (MQ=255)
                                                        tGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCGTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  2:3036559/100‑1 (MQ=255)
                                                         gCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTagatcggaa                                            >  2:2420668/1‑92 (MQ=255)
                                                         gCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                            >  1:3237323/1‑100 (MQ=255)
                                                           aaaTTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                                           <  1:877622/99‑1 (MQ=255)
                                                           aaaTTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                                          <  2:107808/100‑1 (MQ=255)
                                                            aaTTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa                                         <  2:877360/100‑1 (MQ=255)
                                                             aTTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag                                        <  1:2536641/100‑1 (MQ=255)
                                                               tCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag                                      <  2:771543/100‑1 (MQ=255)
                                                                cAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga                                     <  1:898533/100‑1 (MQ=255)
                                                                cAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga                                     <  1:234274/100‑1 (MQ=255)
                                                                 aaTACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCTATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGAc                                    >  2:1836600/1‑100 (MQ=255)
                                                                         gCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            <  2:2947838/100‑1 (MQ=255)
                                                                          cA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa                           >  2:407211/1‑100 (MQ=255)
                                                                          cA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa                           >  1:845539/1‑100 (MQ=255)
                                                                             gAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                         <  2:2027922/100‑1 (MQ=255)
                                                                              ccgatctcgtaggCCTGATAAGCGTAGCGCATCAGGCctacgaaaatcttgcaatgtattgaattttcatgattttgtaggccggataaggcgttcacgc                         <  2:1513200‑M1/94‑64 (MQ=255)
                                                                              aaT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg                        <  2:222898/100‑1 (MQ=255)
                                                                              aaT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg                        <  1:178610/100‑1 (MQ=255)
                                                                                t‑tttCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAt                      <  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)
                                                                                                      |                                                                                                  
TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  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: ‑