Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 860,151 A→G 16.3% intergenic (‑123/+23) 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,1510AG16.3% 346.2 / 56.1 163intergenic (‑123/+23)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/66);  new base G (14/13);  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 = 4.66e-01

Only 100 of 180 total aligned reads displayed.

GAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAG  >  NC_000913/860052‑860249
                                                                                                   |                                                                                                  
gAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGcga                                                                                                    >  1:3643041/1‑100 (MQ=255)
gAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGcga                                                                                                    >  1:2148530/1‑100 (MQ=255)
 aaaCCGTAACAGTGTAATAACAATGTGACGCAGAGTATAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGcgac                                                                                                   <  1:558821/100‑1 (MQ=255)
  aaCCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGcgacg                                                                                                  <  1:2720691/100‑1 (MQ=255)
    ccGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGcgacga                                                                                                 >  2:2956064/1‑99 (MQ=255)
    ccGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCgacgac                                                                                                <  2:2065567/100‑1 (MQ=255)
     cGTAACAGTGTAATAACAATGTTACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGacgacg                                                                                               <  1:1146106/100‑1 (MQ=255)
     cGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAATGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGacgacg                                                                                               <  2:2166005/100‑1 (MQ=255)
      gTAACAGTGTAATAACAATGTGACGCAGAGCACAAAATATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGAcgc                                                                                              >  1:2761388/1‑100 (MQ=255)
        aaCAGAGTAATAACAATATGACGCAGAGCACAAATTATATTTAGAATGAAAGTAAGGATTAAATTGATTATGTGTATGATTTAGCCCGGCGacgacg                                                                                               <  2:3021167/97‑1 (MQ=255)
          cAGTTTAATAACAATGTGACGTAGAGCACAAATTATATTTCTAATGAAATTAATGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACgccgcc                                                                                          <  1:2175309/100‑1 (MQ=255)
          cAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACgccgcc                                                                                          <  2:1616322/100‑1 (MQ=255)
            gtgtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGCTTTAGCCCGGCGACGACGCCGCCgg                                                                                        >  2:739689/1‑100 (MQ=255)
             tgtAATAAGAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTCGCCCGGCGACGACGCCGCCggg                                                                                       <  2:1186105/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)
               taataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGacgacgtgctgctgcc                                                                                     >  1:1674469/1‑90 (MQ=255)
               taataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGct                                                                                     <  1:3869218/100‑2 (MQ=255)
               taataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGcc                                                                                     >  2:1406868/1‑100 (MQ=255)
                 ataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGCTga                                                                                   <  2:402043/100‑1 (MQ=255)
                 ataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCga                                                                                   <  2:159796/100‑1 (MQ=255)
                  taaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgag                                                                                  >  1:1869051/1‑100 (MQ=255)
                   aaCAATGTGACGCTGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGTGCCgag                                                                                  <  2:3925037/99‑1 (MQ=255)
                   aaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgagg                                                                                 >  1:3753260/1‑100 (MQ=255)
                    aCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgagga                                                                                <  2:1141324/100‑1 (MQ=255)
                      aaTGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATTATTTAGCCCGGCGGCGTCGTCGCCGGGctaaatca                                                                               <  2:1207285/99‑8 (MQ=37)
                      aaTGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGctaaatcat                                                                              <  2:182359/100‑9 (MQ=37)
                         gtgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTa                                                                           >  1:1676241/1‑100 (MQ=255)
                         gtgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGACGAGGAGATTa                                                                           >  1:1803664/1‑100 (MQ=255)
                          tgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTAc                                                                          >  2:1032527/1‑100 (MQ=255)
                          tgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTAc                                                                          >  2:774781/1‑100 (MQ=255)
                          tgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGAGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTAc                                                                          >  1:2922837/1‑100 (MQ=255)
                            aCGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAg                                                                        <  1:3576249/100‑1 (MQ=255)
                             cGCAGAGAAAAAATTATATTTCGAATGATAGTATGGATGAAATGGGTGATGTGAATGATTTAGCCCGGCGACGATGCCGCCGGGCCGAGGAGATTACAGc                                                                       <  2:1120451/100‑1 (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)
                                      caAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCAGCCGGGCCGAGGAGATTACAGCATATGTTCa                                                              >  1:23341/1‑100 (MQ=255)
                                       aaaTTCTATTCCGAATGCAATTAAGGATGAAATGTATGATTTCAATGATTTTGCCCGGCGACGACGCCGCCGTGCCGAGCAGATTACAGc                                                                       <  2:2415476/90‑1 (MQ=255)
                                       aaaTTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAg                                                             <  1:2127101/100‑1 (MQ=255)
                                       aaaTTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAg                                                             <  2:234993/100‑1 (MQ=255)
                                         aTTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTa                                                           <  2:3355888/100‑1 (MQ=255)
                                            atatTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAACCCGGCGACGACGCCGCCGGGCCGAGGAGATTAAAGCATATGTTCAGTACgg                                                        >  2:3509421/1‑100 (MQ=255)
                                             tatTTCTAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACggg                                                       <  1:3774884/100‑1 (MQ=255)
                                              atTCCGAATGCAAGTAAGGATAAAATTTATGATGTGAATAATTTAGCCCGGCGACGACGACGCCGGGCCTAGGATATTACAGCATATGTTCAGTACGGGc                                                      <  2:469297/100‑1 (MQ=255)
                                               tttCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCg                                                     >  1:2469894/1‑100 (MQ=255)
                                                 tCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCga                                                    >  2:1738821/1‑99 (MQ=255)
                                                  cGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCgatg                                                  >  2:3036399/1‑100 (MQ=255)
                                                    aaTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCgatgat                                                <  2:1257365/100‑1 (MQ=255)
                                                     aTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGata                                               <  2:89144/100‑1 (MQ=255)
                                                      tctaaGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGatat                                              <  2:828268/97‑1 (MQ=255)
                                                       gAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGAAGATTACAGCATATGTTCAGTACGGGCGATGATATc                                             <  2:1561400/100‑1 (MQ=255)
                                                       caaaGTAAGGAGCAAAGTGATGAGGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGTGATTATAGAATATGTTTAGTACGGGCgatgat                                                <  2:2679954/96‑1 (MQ=255)
                                                        aaaGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCg                                            <  1:3850239/100‑1 (MQ=255)
                                                         aaGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGatataga                                           >  1:828278/1‑97 (MQ=255)
                                                         aaGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGt                                           >  2:2596330/1‑100 (MQ=255)
                                                         aaGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGt                                           <  1:3664184/100‑1 (MQ=255)
                                                           gTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTcc                                         <  2:3753503/100‑1 (MQ=255)
                                                              aGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGa                                      <  2:1123418/100‑1 (MQ=255)
                                                               ggATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAg                                     <  1:2665084/100‑1 (MQ=255)
                                                               ggATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAg                                     >  1:1701296/1‑100 (MQ=255)
                                                               ggAAGAAATTTATGATCTGAATGATATAGCCGGGCGACGACGCCGCCGGGGCGAGGAGATTACAGCACATATTGCATACGGGCGATAAGATCGTctaaag                                     >  2:2162607/1‑95 (MQ=255)
                                                                gATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGc                                    <  1:3144564/100‑1 (MQ=255)
                                                                 aTGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCg                                   >  1:1980846/1‑100 (MQ=255)
                                                                  tGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGt                                  >  2:3137663/1‑100 (MQ=255)
                                                                    aaaTTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCt                                >  1:3551390/1‑100 (MQ=255)
                                                                        tgatgaTGTGAATGACTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGgag                            <  1:2533534/100‑1 (MQ=255)
                                                                            gattTGAATGATTGAGCCCGGCGATGACGCCGCCTGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGagagag                        <  2:2464658/100‑1 (MQ=255)
                                                                              tgtgAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGGGCGTCTGGAGAGAgc                       >  2:2291241/1‑99 (MQ=255)
                                                                                 gAATGATTTAGCCCGGCGACGACGCCGCCGGTCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCgg                   <  1:221399/100‑1 (MQ=255)
                                                                                 gAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGcgcg                    >  1:772563/1‑99 (MQ=255)
                                                                                  gcgtcctcggCCCGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGt                  <  1:2781377/91‑1 (MQ=38)
                                                                                   aTGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTg                 >  1:1845916/1‑100 (MQ=255)
                                                                                     gATTTAGCACGGCGACGACGCTGGCGGGCGGAGGAGATTACAGCATATGTTAAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaa               <  2:1047791/100‑1 (MQ=255)
                                                                                      cctcggCCCGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaag              >  2:3403223/6‑100 (MQ=38)
                                                                                      aTTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaag              <  1:3829273/100‑1 (MQ=255)
                                                                                       tttAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaaga             >  1:723281/1‑100 (MQ=255)
                                                                                       ttaaGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaaga             <  2:430786/97‑1 (MQ=255)
                                                                                       ctcggCCCGGCGGCGTCGTCGCCGGGCCGAGGAGCTTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaaga             >  1:3101536/5‑100 (MQ=38)
                                                                                        ttAGCCCGGCGATGACGCCGGCGGGCCGAGAAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaag              <  2:3583908/98‑1 (MQ=255)
                                                                                         cggCCCGGCGGCGTCGTCGCCGGGCCGAGGAGCTTACAGCATATGTTCAGTACGGGCGATGATATCATCCTGAACGTCTGAAGAGAGCGAGGTaaagaa            >  1:1778391/3‑99 (MQ=37)
                                                                                         cggCCCGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTAAGCGTCTGGAGAGAACGCGGTgaagaa            >  1:2433570/3‑99 (MQ=38)
                                                                                           ncccGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcgc         >  1:249174/2‑100 (MQ=255)
                                                                                           gCCCGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTCCGGGCGATGATATCGTCCTGAGCGTCTGTAGCGCGCGCGGTGAAGAAcgc         >  2:594388/1‑100 (MQ=255)
                                                                                            cccGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcgcg        >  1:2042144/1‑100 (MQ=255)
                                                                                             cccgCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCTTCTGGAGAGAGCGCGGTGAATAAcgc         <  2:3221289/95‑1 (MQ=255)
                                                                                             ccGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCgg       >  2:2798253/1‑100 (MQ=255)
                                                                                             ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCgg       >  1:1155483/1‑100 (MQ=255)
                                                                                             ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCgg       <  2:1308961/100‑1 (MQ=255)
                                                                                             ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCgg       <  1:2309681/100‑1 (MQ=255)
                                                                                             ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCgg       <  1:2157057/100‑1 (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)
                                                                                                   |                                                                                                  
GAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAG  >  NC_000913/860052‑860249

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