Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 860,154 A→T 17.1% intergenic (‑126/+20) ybiV ← / ← ybiW sugar phosphatase; fructose‑1‑P/ribose‑5‑P/glucose‑6‑P phosphatase/putative pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913860,1540AT17.1% 360.2 / 65.9 163intergenic (‑126/+20)ybiV/ybiWsugar phosphatase; fructose‑1‑P/ribose‑5‑P/glucose‑6‑P phosphatase/putative pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base A (70/65);  new base T (14/14);  total (84/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.44e-01

Only 100 of 175 total aligned reads displayed.

CCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCC  >  NC_000913/860056‑860251
                                                                                                  |                                                                                                 
ccGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGcgacga                                                                                                   <  1:3745939/99‑1 (MQ=255)
ccGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGcgacga                                                                                                   >  2:2956064/1‑99 (MQ=255)
 cGTAACAGTGTAATAACAATGTTACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGacgacg                                                                                                 <  1:1146106/100‑1 (MQ=255)
  gTAACAGTGTAATAACAATGTGACGCAGAGCACAAAATATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGAcgc                                                                                                >  1:2761388/1‑100 (MQ=255)
    aaCAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGAcgccg                                                                                              <  2:381768/100‑1 (MQ=255)
    aaCAGAGTAATAACAATATGACGCAGAGCACAAATTATATTTAGAATGAAAGTAAGGATTAAATTGATTATGTGTATGATTTAGCCCGGCGacgacg                                                                                                 <  2:3021167/97‑1 (MQ=255)
      cAGTTTAATAACAATGTGACGTAGAGCACAAATTATATTTCTAATGAAATTAATGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACgccgcc                                                                                            <  1:2175309/100‑1 (MQ=255)
       aGTGTAATAACAATGAGACGCAGAGCACAAATGATATTGCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGccgccg                                                                                           >  2:570890/1‑100 (MQ=255)
         tgtAATAACACTGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATATGAATGATTTAGCCCGGCGACGACGCCGCCggg                                                                                         <  2:2566626/100‑1 (MQ=255)
          gtAATAACAATGTGACGCAGATCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGc                                                                                        <  1:979982/100‑1 (MQ=255)
          gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGc                                                                                        <  2:2889942/100‑1 (MQ=255)
          gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGc                                                                                        <  1:262558/100‑1 (MQ=255)
          gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGc                                                                                        <  1:2504195/100‑1 (MQ=255)
          gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGc                                                                                        >  1:3753542/1‑100 (MQ=255)
          gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGGAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGGCGCCGGGc                                                                                        <  2:1836943/100‑1 (MQ=255)
            aataaCAATGCGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCg                                                                                      <  2:2806394/100‑1 (MQ=255)
             ataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGCTga                                                                                     <  2:402043/100‑1 (MQ=255)
             ataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCga                                                                                     <  2:159796/100‑1 (MQ=255)
              taaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgag                                                                                    >  1:1869051/1‑100 (MQ=255)
              naaCAATGTGACGCACAGCACAAATAATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTCAGCGCGGagacga                                                                                                   >  2:2188528/2‑85 (MQ=255)
               aaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATTAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGctaaat                                                                                   <  2:3653150/100‑6 (MQ=38)
                aCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgagga                                                                                  <  1:2482241/100‑1 (MQ=255)
                  aaTGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATTATTTAGCCCGGCGGCGTCGTCGCCGGGctaaatca                                                                                 <  2:1207285/99‑8 (MQ=37)
                   aTGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGctaaatca                                                                                 <  2:1604366/98‑8 (MQ=37)
                    tgtgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGAt                                                                               <  2:1268864/99‑1 (MQ=255)
                     gtgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTa                                                                             >  2:620851/1‑100 (MQ=255)
                     gtgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTa                                                                             >  1:1676241/1‑100 (MQ=255)
                      tgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTAc                                                                            >  1:2668296/1‑100 (MQ=255)
                      tgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTAc                                                                            >  2:774781/1‑100 (MQ=255)
                         cGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGCTTACAGc                                                                         >  2:1137977/1‑100 (MQ=255)
                         cGCAGAGAAAAAATTATATTTCGAATGATAGTATGGATGAAATGGGTGATGTGAATGATTTAGCCCGGCGACGATGCCGCCGGGCCGAGGAGATTACAGc                                                                         <  2:1120451/100‑1 (MQ=255)
                          gCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCa                                                                        >  2:1822523/1‑100 (MQ=255)
                          gCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGAATTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCa                                                                        >  2:2753079/1‑100 (MQ=255)
                           cAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATTATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGc                                                                         <  1:1411794/98‑1 (MQ=255)
                             gagCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCatat                                                                     <  1:620858/100‑1 (MQ=255)
                              agCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATg                                                                    <  1:3192537/100‑1 (MQ=255)
                              agCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATg                                                                    >  2:2044328/1‑100 (MQ=255)
                                cacaAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGtt                                                                  >  1:1205711/1‑100 (MQ=255)
                                     aTTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTa                                                             <  2:3355888/100‑1 (MQ=255)
                                       tttatTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACg                                                           <  1:2567766/98‑1 (MQ=255)
                                        atatTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAACCCGGCGACGACGCCGCCGGGCCGAGGAGATTAAAGCATATGTTCAGTACgg                                                          >  2:3509421/1‑100 (MQ=255)
                                         tatTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACggg                                                         >  1:3622174/1‑100 (MQ=255)
                                           tttCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCg                                                       >  1:380902/1‑100 (MQ=255)
                                             tCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCgat                                                     >  1:3143921/1‑100 (MQ=255)
                                              cGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCgatg                                                    >  2:3036399/1‑100 (MQ=255)
                                                 aTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGata                                                 <  2:89144/100‑1 (MQ=255)
                                                  tctaaGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGatat                                                <  2:828268/97‑1 (MQ=255)
                                                   caaaGTAAGGAGCAAAGTGATGAGGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGTGATTATAGAATATGTTTAGTACGGGCgatgat                                                  <  2:2679954/96‑1 (MQ=255)
                                                    aaaGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCg                                              <  1:2289127/100‑1 (MQ=255)
                                                    aaaGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCg                                              <  1:3850239/100‑1 (MQ=255)
                                                    aaaGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGAATATGTTCAGTACGGGCGATGATATCg                                              >  1:1655768/1‑100 (MQ=255)
                                                     aaGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGatataga                                             >  1:828278/1‑97 (MQ=255)
                                                     aaGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGt                                             <  1:3664184/100‑1 (MQ=255)
                                                     aaGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGt                                             >  2:2596330/1‑100 (MQ=255)
                                                          aGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGa                                        <  2:1123418/100‑1 (MQ=255)
                                                          aGGATGAAATTGATGATGTAAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGa                                        <  2:1668123/100‑1 (MQ=255)
                                                           ggATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAg                                       >  1:1701296/1‑100 (MQ=255)
                                                           ggATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAg                                       >  2:2455080/1‑100 (MQ=255)
                                                             aTGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCg                                     >  1:1980846/1‑100 (MQ=255)
                                                              tGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGt                                    >  2:3137663/1‑100 (MQ=255)
                                                                aaaTTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCt                                  >  1:3551390/1‑100 (MQ=255)
                                                                     gatgatGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGgaga                             >  1:2105025/1‑100 (MQ=255)
                                                                      atgattTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGgaga                             <  1:1798267/99‑1 (MQ=255)
                                                                       tgatgTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGgagaga                           >  1:2031033/1‑100 (MQ=255)
                                                                        gatgTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATAAGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGagagag                          >  1:2179917/1‑100 (MQ=255)
                                                                          tgtgAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGGGCGTCTGGAGAGAgc                         >  2:2291241/1‑99 (MQ=255)
                                                                            tgAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGcgcg                      <  2:3368069/100‑1 (MQ=255)
                                                                            tgAATGATTTAGCCCGGCGAAGACGCCGCCGGGCCGAGGAAATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAAAGcgcg                      >  2:117816/1‑100 (MQ=255)
                                                                             gAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGcgcg                      >  1:772563/1‑99 (MQ=255)
                                                                             gAAAGATTTAGCCCGGCAACAACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATAGCGTACTGACAGTCTGGAGAGAGCGCgg                     >  2:2163345/1‑100 (MQ=255)
                                                                              gcgtcctcggCCCGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGt                    <  1:2781377/91‑1 (MQ=38)
                                                                                tGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTga                  <  1:2913473/100‑1 (MQ=255)
                                                                                 gATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaa                 >  2:2141381/1‑100 (MQ=255)
                                                                                 gATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaa                 >  2:433288/1‑100 (MQ=255)
                                                                                 gATTTAGCACGGCGACGACGCTGGCGGGCGGAGGAGATTACAGCATATGTTAAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaa                 <  2:1047791/100‑1 (MQ=255)
                                                                                  aTTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaa                 >  1:397922/1‑99 (MQ=255)
                                                                                   tttAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaaga               >  1:723281/1‑100 (MQ=255)
                                                                                   ttaaGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaaga               <  2:430786/97‑1 (MQ=255)
                                                                                    ttAGCCCGGCGATGACGCCGGCGGGCCGAGAAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaag                <  2:3583908/98‑1 (MQ=255)
                                                                                     cggCCCGGCGGCGTCGTCGCCGGGCCGAGGAGCTTACAGCATATGTTCAGTACGGGCGATGATATCATCCTGAACGTCTGAAGAGAGCGAGGTaaagaa              >  1:1778391/3‑99 (MQ=37)
                                                                                     cggCCCGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTAAGCGTCTGGAGAGAACGCGGTgaagaa              >  1:2433570/3‑99 (MQ=38)
                                                                                       gCCCGGCGGCGTCGTCGCCGGGCTGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGAGCGGTGAAGAAcgc           >  2:2391966/1‑100 (MQ=255)
                                                                                       gCCCGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATAATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcgc           >  1:2438529/1‑100 (MQ=255)
                                                                                       gCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcgc           >  1:3190916/1‑100 (MQ=255)
                                                                                        cccGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGGGGTGAAGAAcgcg          >  2:102601/1‑100 (MQ=255)
                                                                                        cccGGCGGCGTCGTCGCCGGGCCGAGGAGATCACAGCATATGTTTAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcg            >  2:1476281/1‑98 (MQ=255)
                                                                                         cccgCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCTTCTGGAGAGAGCGCGGTGAATAAcgc           <  2:3221289/95‑1 (MQ=255)
                                                                                         ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCgg         >  1:1155483/1‑100 (MQ=255)
                                                                                         ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCgg         <  2:1308961/100‑1 (MQ=255)
                                                                                         ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCgg         >  1:3851611/1‑100 (MQ=255)
                                                                                         ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCgg         <  1:2157057/100‑1 (MQ=255)
                                                                                         ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACAGGAGATGATATCAACCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCgg         >  1:231672/1‑100 (MQ=255)
                                                                                          cGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGa        >  2:3408188/1‑100 (MQ=255)
                                                                                          cGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGa        <  1:1554962/100‑1 (MQ=255)
                                                                                            gCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAAt      >  2:969588/1‑100 (MQ=255)
                                                                                            gCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAATACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAAt      <  1:858033/100‑1 (MQ=255)
                                                                                             cgacgaCGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATa     <  2:2914413/100‑1 (MQ=255)
                                                                                              gacgacGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAg    >  1:183918/1‑100 (MQ=255)
                                                                                                cgacgCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGcc  >  2:1330602/1‑100 (MQ=255)
                                                                                                 gacgGCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGcc  <  1:3031088/99‑1 (MQ=255)
                                                                                                  |                                                                                                 
CCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCC  >  NC_000913/860056‑860251

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