Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,973,096 Δ1 bp 15.3% coding (575/1068 nt) rffG → dTDP‑glucose 4,6‑dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,973,0950G.15.3% 281.2 / 48.2 98coding (574/1068 nt)rffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base G (44/39);  new base . (9/6);  total (53/45)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

CCTTATTCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGCAACGGGCAGCAAATCCGT  >  NC_000913/3972996‑3973193
                                                                                                   |                                                                                                  
ccTTATTCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACg                                                                                                    >  1:2525072/1‑100 (MQ=255)
  ttATTCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGc                                                                                                  <  2:2106551/100‑1 (MQ=255)
    aTTCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGccc                                                                                                >  1:1012257/1‑100 (MQ=255)
     ttCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCt                                                                                               >  2:855230/1‑100 (MQ=255)
      tCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCtt                                                                                              <  2:1343957/100‑1 (MQ=255)
      tCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCtt                                                                                              >  1:2257513/1‑100 (MQ=255)
      tCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGCCCTTa                                                                                             >  2:631009/1‑100 (MQ=255)
      tCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCGATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGCCCtt                                                                                              >  1:1863687/1‑99 (MQ=255)
       ccGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTa                                                                                             >  2:890662/1‑100 (MQ=255)
       ccGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTa                                                                                             >  1:709275/1‑100 (MQ=255)
       ccGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTa                                                                                             >  1:1827398/1‑100 (MQ=255)
       ccGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTa                                                                                             >  1:1980730/1‑100 (MQ=255)
         gcgTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTAcc                                                                                           <  2:1726108/100‑1 (MQ=255)
          cgTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGCCCTTACCAc                                                                                         <  1:78319/100‑1 (MQ=255)
           gTCAAAAGCCAGCGGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTAccgc                                                                                         <  2:697975/100‑3 (MQ=255)
           gTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCAc                                                                                         <  2:957288/100‑1 (MQ=255)
            tCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCGTACCAc                                                                                         >  1:1860549/1‑99 (MQ=255)
              aaaaGCCAGCAGCGACCATCTGGTGAGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGCCCTTACCACTTTc                                                                                     >  1:2178739/1‑100 (MQ=255)
               aaaGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTTCTCTAATAACTACGGCCCGTACCACGTTc                                                                                     >  1:859062/1‑100 (MQ=255)
                naGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACCTTcc                                                                                    >  1:100/2‑100 (MQ=255)
                aaGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTcc                                                                                    <  2:2257528/100‑1 (MQ=255)
                aaGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTcc                                                                                    >  2:1636408/1‑100 (MQ=255)
                 aGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCg                                                                                   <  1:1659865/100‑1 (MQ=255)
                  gCCAGCAGCGACCATCTGGTGCGCGCCTGGTTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGCCCTTACCACTTTCCGGa                                                                                 >  1:2320856/1‑100 (MQ=255)
                  gCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCgg                                                                                  >  2:2534181/1‑100 (MQ=255)
                   ccAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGa                                                                                 >  2:810293/1‑100 (MQ=255)
                   ccAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGa                                                                                 >  2:1692764/1‑100 (MQ=255)
                    cagcagCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGaa                                                                                >  1:603035/1‑100 (MQ=255)
                    cagcagCGACCATCTGGTGCGCGCCTGGCTGCAGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGCTCTTACCACTTTCCGGaaa                                                                               >  2:2457251/1‑100 (MQ=255)
                     agcagcGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGaaa                                                                               >  1:489028/1‑100 (MQ=255)
                     agcagcGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGaaa                                                                               <  1:224566/100‑1 (MQ=255)
                      gcagcGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGaaaa                                                                              >  2:1405877/1‑100 (MQ=255)
                      gcagcGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGaaaa                                                                              >  1:606284/1‑100 (MQ=255)
                         gcGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACt                                                                           >  2:1606714/1‑100 (MQ=255)
                          cGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTg                                                                          <  2:510720/100‑1 (MQ=255)
                           gACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCATACCACTTTCCGGAAAAACt                                                                           >  1:1991834/1‑98 (MQ=255)
                            aCCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGan                                                                        <  2:1243156/100‑2 (MQ=255)
                               aTCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATccc                                                                     <  2:2475972/100‑1 (MQ=255)
                                ncTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCg                                                                    >  1:2184254/2‑100 (MQ=255)
                                  tGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACTAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCt                                                                  <  1:2453365/100‑1 (MQ=255)
                                  tGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCt                                                                  >  2:414392/1‑100 (MQ=255)
                                    gTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCtga                                                                <  2:2031389/100‑1 (MQ=255)
                                      gcgcccCTGGCTGCGGACATATGGTCTGCAGACGCTTACCACCAACTGATCGCATGACTACGGACCTTAACCCTTTCCGGAAAAACTGCTCCAGCtgct                                                               <  2:764913/99‑3 (MQ=255)
                                          gcCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATcc                                                          <  2:1012260/100‑1 (MQ=255)
                                          gcCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCt                                                         <  2:2108483/100‑1 (MQ=255)
                                           ccTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCCCATTCCGGGAAAACTGATCCCGCTGATGATc                                                           >  1:1006208/1‑98 (MQ=255)
                                              ggCTGCGGACCTATAGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCaa                                                      <  1:1719918/100‑1 (MQ=255)
                                               gCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAAcg                                                    >  1:268132/1‑100 (MQ=255)
                                                cTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAAcg                                                    <  1:1830024/100‑1 (MQ=255)
                                                cTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAAcg                                                    <  1:855227/100‑1 (MQ=255)
                                                cTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAAcg                                                    <  2:1192409/100‑1 (MQ=255)
                                                 tGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCTAATAACTACGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAAcgcg                                                  <  1:1738573/100‑1 (MQ=255)
                                                 tGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAAcgc                                                   <  1:716733/100‑1 (MQ=255)
                                                  gCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGCTGATCCTCAAcgcg                                                  >  1:2539046/1‑100 (MQ=255)
                                                     gACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTgg                                              >  2:873754/1‑100 (MQ=255)
                                                      aCCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTgg                                              <  1:226378/100‑1 (MQ=255)
                                                      aCCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGc                                             <  2:1365830/100‑1 (MQ=255)
                                                       ccTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGc                                             >  2:647592/1‑100 (MQ=255)
                                                         tATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCCCAACGCGCTg                                               >  1:1723675/1‑96 (MQ=255)
                                                         tATGGTCTGCCGACGCTTATCACCAACTGCCCGAATAACTACGGCCCTTACCACTTTCCTGAAAAACTGATCCCGCTGCTGATCCTCAACGCGCTGGCgg                                           >  1:45525/1‑100 (MQ=255)
                                                          aTGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCggg                                          <  2:574166/100‑1 (MQ=255)
                                                           tGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGt                                         <  2:2525091/100‑1 (MQ=255)
                                                            tctCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTa                                        <  2:1149944/98‑1 (MQ=255)
                                                             gtttggcGATGCTTATCAAAAAGTGCTGGAATAATTACGCCCTTACAACTTTCCGGAAAAACTGATCCCGTTGATGATCCTCAACGCGCTTGCGGGTaaa                                      >  1:2424588/7‑100 (MQ=255)
                                                             gTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTaa                                       >  1:2225968/1‑100 (MQ=255)
                                                             gTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTaaa                                      >  2:317326/1‑100 (MQ=255)
                                                               cTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTaaga                                     >  1:1149940/1‑98 (MQ=255)
                                                                 gCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCg                                   <  2:1264248/100‑1 (MQ=255)
                                                                 gCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCg                                   <  2:140858/100‑1 (MQ=255)
                                                                   cGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCt                                 >  1:1836953/1‑100 (MQ=255)
                                                                        cTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCgg                            >  2:2160915/1‑100 (MQ=255)
                                                                        cTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCgg                            <  2:1419478/100‑1 (MQ=255)
                                                                         ttATCACCAACTGCTCGAATAACTACGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGta                          <  2:2521521/100‑1 (MQ=255)
                                                                          tATCACCAACTGCTCGAATAACTACGGCCCTTACCTCTTTCCGGAAAAACTGATCCCGCTGATGATCCCCAACGCGCTGGCGGGTCAATCGCTccc                              >  1:48826/1‑96 (MQ=255)
                                                                          tATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGta                          >  1:449504/1‑100 (MQ=255)
                                                                          tATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGta                          <  1:1174122/100‑1 (MQ=255)
                                                                              aCCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATg                      >  2:2185176/1‑100 (MQ=255)
                                                                              aCCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATg                      <  2:652710/100‑1 (MQ=255)
                                                                              aCCAACTGCTCGAATAACTACGGCCCTTACCACGTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATg                      <  1:2271088/100‑1 (MQ=255)
                                                                               ccAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATgg                     >  2:757726/1‑100 (MQ=255)
                                                                                cAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGc                    >  1:1732405/1‑100 (MQ=255)
                                                                                    tGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGCAACg                >  2:312788/1‑100 (MQ=255)
                                                                                    tGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGCAACg                <  2:412462/100‑1 (MQ=255)
                                                                                     gCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCTCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGCAACg                <  1:1105632/99‑1 (MQ=255)
                                                                                     gCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGCAACgg               <  1:2538322/100‑1 (MQ=255)
                                                                                        cGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGCAACGGgca            <  2:589137/100‑1 (MQ=255)
                                                                                        cGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGCAACGGgca            >  2:897124/1‑100 (MQ=255)
                                                                                        cGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGCAACGGgc             >  1:1888492/1‑99 (MQ=255)
                                                                                        cGAATAACTACGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGCAACGGgcag           <  2:2245354/100‑1 (MQ=255)
                                                                                          aaTAACTACGGCCCTTACCACTTTCCGGAAAAGCTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGCAACGGgcagc          <  1:371926/100‑1 (MQ=255)
                                                                                          aaTAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGCAACGGgcag           >  1:406866/1‑99 (MQ=255)
                                                                                          aaTAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGCAACGGgcag           >  1:407816/1‑99 (MQ=255)
                                                                                           aTAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGCAACGGgcagca         <  2:1343771/100‑1 (MQ=255)
                                                                                            tAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTTGCGGGGAAATGGCTGCCGGTAAATGGCAACgg               >  1:1407692/1‑93 (MQ=255)
                                                                                             aaCTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGCAACGGGCAGCaaa       >  2:2243588/1‑100 (MQ=255)
                                                                                             aaCTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGCAACGGGCAGCaaa       <  1:375226/100‑1 (MQ=255)
                                                                                              aCTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGCAACGGGCAGCAAAt      >  1:1780920/1‑100 (MQ=255)
                                                                                                tACGGCCCTTACCACAAACAGAAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGCAACGGGCAGCAAATcc    <  1:1756196/100‑1 (MQ=255)
                                                                                                  cGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCTGGTAAATCGCTGCCGGTATATGGCAACGCGCAGCAAATccat  <  2:1752769/100‑3 (MQ=255)
                                                                                                   |                                                                                                  
CCTTATTCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTATATGGCAACGGGCAGCAAATCCGT  >  NC_000913/3972996‑3973193

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