Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,017 C→T 18.6% intergenic (‑227/+81) 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,0170CT18.6% 524.3 / 91.5 213intergenic (‑227/+81)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (66/100);  new base T (35/3);  total (110/103)
Fisher's exact test for biased strand distribution p-value = 1.25e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 216 total aligned reads displayed.

GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116
                                                                                                    |                                                                                                      
gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                                                                        <  2:756339/100‑1 (MQ=255)
gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                                                                        <  2:1717408/100‑1 (MQ=255)
      ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCa                                                                                                  <  1:1722107/100‑1 (MQ=255)
        gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                                                   <  1:2540662/97‑1 (MQ=255)
        gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Ga                                                                                               <  1:2726085/100‑1 (MQ=255)
          gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGACCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Gaa                                                                                              >  2:1000886/1‑99 (MQ=255)
             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAt‑t‑tt                                                                                        >  1:96292/1‑100 (MQ=255)
              aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTc                                                                                       <  2:2553549/100‑1 (MQ=255)
              aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTc                                                                                       >  2:445149/1‑100 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCg                                                                                      <  1:554843/100‑1 (MQ=255)
                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTa                                                                                    >  1:71821/1‑100 (MQ=255)
                  ggcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAg                                                                                   >  1:2269458/8‑100 (MQ=18)
                  aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAg                                                                                   <  2:56786/100‑1 (MQ=255)
                  aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAg                                                                                   <  1:3241177/100‑1 (MQ=255)
                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGG                                                                                  <  1:120160‑M1/100‑1 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCC                                                                                <  2:3057609‑M1/100‑1 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCC                                                                                <  1:2453805‑M1/100‑1 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCC                                                                                >  1:2298712‑M1/1‑100 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGA                                                                             >  2:1112880‑M1/1‑100 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTG                                                                              <  1:1419994‑M1/99‑1 (MQ=255)
                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATa                                                                           <  2:496698/100‑1 (MQ=255)
                           aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCtaataa                                                                          >  2:2980698/1‑100 (MQ=255)
                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAg                                                                    <  1:3028648/100‑1 (MQ=255)
                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAg                                                                    <  1:2768602/100‑1 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAgcg                                                                  >  2:2536623/1‑100 (MQ=255)
                                    gttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGCA‑ATCT‑C‑CGTGTAGGCCTGATAAGCGTAgcgc                                                                 >  1:1381591/10‑100 (MQ=18)
                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCa                                                                <  1:2319709‑M2/100‑2 (MQ=255)
                                     atttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCC‑CTTGTAGGCCTGATAAGCGTAGCGCAt                                                               >  1:459037/9‑100 (MQ=18)
                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCatc                                                              <  1:1243170‑M2/100‑4 (MQ=255)
                                        tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCatca                                                             >  1:2578396‑M2/1‑96 (MQ=255)
                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCatcag                                                            <  1:2027936‑M2/100‑6 (MQ=255)
                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑T‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  2:2064445/4‑100 (MQ=21)
                                           ntgCCGGAAGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATAGTGCAAATACAATACATTGCA‑GAAT‑C‑TCCGTAGGCCGGATAAGCGGAGCGCatcaggc                                                          <  2:1296808‑M2/99‑8 (MQ=255)
                                           aggccGGATAAGGCGTTAACGACGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCatcaggc                                                          <  2:1491592‑M2/98‑8 (MQ=255)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑C‑CTTGTAGGCCTGATAAGCGTACCGCATCAGGc                                                          >  2:2132410/3‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCGTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GGCT‑C‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  1:1827863/3‑100 (MQ=21)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑A‑TTGGTAGGCCTGATAAGCGCAGCGCATCAgg                                                           >  1:1066213/3‑99 (MQ=21)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑A‑TTGGTAGGCCTGATAAGCGCAGCGCATCAGGc                                                          >  2:1433463/3‑100 (MQ=25)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGA‑T‑CATGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  1:2929041/3‑100 (MQ=14)
                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATgcgctccgcttatcaggc                                                          <  1:739240‑M1/100‑19 (MQ=255)
                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATgcgctccgcttatcaggc                                                          <  2:96290‑M1/100‑19 (MQ=255)
                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATgcgctccgcttatcaggc                                                          <  1:1186655‑M1/100‑19 (MQ=255)
                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATgcgctccgcttatcaggc                                                          <  1:1992686‑M1/100‑19 (MQ=255)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                         >  1:651762/1‑100 (MQ=255)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                         >  1:346974/1‑100 (MQ=255)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                         >  2:334523/1‑100 (MQ=255)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                         >  1:1574772/1‑100 (MQ=255)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                         >  1:2514116/1‑100 (MQ=255)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCC‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                        >  1:758786/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑C‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  2:628195/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCA‑GGAT‑A‑CGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  1:1097111/2‑100 (MQ=17)
                                            ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                         >  2:3142817/2‑100 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑C‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                        >  2:2926657/1‑99 (MQ=18)
                                              ccGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑C‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                       >  2:2161661/1‑98 (MQ=18)
                                               cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATgcgctccgcttatcaggcctac                                                      <  1:2065379‑M1/100‑23 (MQ=255)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑C‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  2:3245867/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GGAA‑A‑CACGTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                                                     >  2:3173169/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑T‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                                                    >  1:150825/7‑98 (MQ=21)
                                                ggATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATgcgctccgcttatcaggcctacg                                                     <  1:1738871‑M1/100‑24 (MQ=255)
                                                 gataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAAT‑CCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgcg                                                   >  2:1654229/6‑97 (MQ=21)
                                                 gataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GGCT‑C‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                    <  1:478085/95‑1 (MQ=21)
                                                 gATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGACgcgctccgcttatcaggcctacga                                                    <  1:41579‑M1/100‑25 (MQ=255)
                                                  ataaggcgTTTACGCCTTATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATGGGC‑C‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttta                                                   >  2:852679/5‑99 (MQ=18)
                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                   <  1:493807/100‑1 (MQ=255)
                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                   >  2:898528/1‑100 (MQ=255)
                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                   <  2:1203661/100‑1 (MQ=255)
                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                   <  1:2862721/100‑1 (MQ=255)
                                                   tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                  <  2:592659/100‑1 (MQ=255)
                                                   tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                  >  1:1655618/1‑100 (MQ=255)
                                                   tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                  >  1:1544192/1‑100 (MQ=255)
                                                       gcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATgcgctccgcttatcaggcctacgaaaattc                                              <  1:3235040‑M1/100‑31 (MQ=255)
                                                        cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                             <  2:540356/100‑1 (MQ=255)
                                                        cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                             <  2:138552/100‑1 (MQ=255)
                                                         gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                            >  2:1570183/1‑100 (MQ=255)
                                                           tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                          >  2:1455542/1‑100 (MQ=255)
                                                            aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑TAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                          >  1:317301/1‑99 (MQ=255)
                                                              cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                       >  2:244358/1‑100 (MQ=255)
                                                               gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                      <  2:1595774/100‑1 (MQ=255)
                                                                ccTTATCCGGCCTAC‑AAAATCGATCAAATACAATAAATTGCA‑GAAT‑A‑TTCGTAGGCCTGATAAGCGGAGCGCATCAAGCATTTTTGCCTTTGCTAtta                                       >  2:2597756/1‑96 (MQ=18)
                                                                    aTCCGGCCTAT‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                 <  1:42761/100‑1 (MQ=255)
                                                                       cGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                              <  2:2269442/100‑1 (MQ=255)
                                                                       cGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                              <  2:1188646/100‑1 (MQ=255)
                                                                        ggCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                             <  2:284799/100‑1 (MQ=255)
                                                                         gCCTAT‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATgg                            <  2:3126902/100‑1 (MQ=255)
                                                                           cTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                          >  1:1621840/1‑100 (MQ=255)
                                                                              c‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTTGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTCTTTTTCTATCTGTCTTGTATCGCtt                         >  2:2004695/1‑98 (MQ=21)
                                                                               ttccgatctGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                       <  1:2420684/92‑1 (MQ=255)
                                                                                aaaaTCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                      <  2:1651247/100‑1 (MQ=255)
                                                                                  aaTCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                    <  2:684105/100‑1 (MQ=255)
                                                                                     cGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                 <  2:176973/100‑1 (MQ=255)
                                                                                       tGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               >  2:1214875/1‑100 (MQ=255)
                                                                                         cAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTTTTTTTTTTTCAATCTTGTAAGGCTTTTAATTATTCGa             >  2:798048/1‑100 (MQ=21)
                                                                                          aaaTTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa             <  1:877622/99‑1 (MQ=255)
                                                                                          aaaTTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt            <  2:107808/100‑1 (MQ=255)
                                                                                           aaTTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa           <  2:877360/100‑1 (MQ=255)
                                                                                            aTTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag          <  1:1132154/100‑1 (MQ=255)
                                                                                             ttCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga         >  1:892122/1‑100 (MQ=255)
                                                                                               cAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga       <  1:234274/100‑1 (MQ=255)
                                                                                                aaTACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCTATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGAc      >  2:1836600/1‑100 (MQ=255)
                                                                                                 aTACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACg     >  1:3094236/1‑100 (MQ=255)
                                                                                                    |                                                                                                      
GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116

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: ‑