Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,074,029 C→G 7.9% I97M (ATC→ATG)  yihV → 6‑deoxy‑6‑sulfofructose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,074,0290CG 7.9% 204.8 / 6.8 63I97M (ATC→ATG) yihV6‑deoxy‑6‑sulfofructose kinase
Reads supporting (aligned to +/- strand):  ref base C (27/31);  new base G (2/3);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCC  >  NC_000913/4073881‑4074176
                                                                                                                                                    |                                                                                                                                                   
cggcAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCcatca                                                                                                                                                    >  2:1266654/1‑150 (MQ=255)
 ggcAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGCCACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAAACTTCGCAATCCGCcatcat                                                                                                                                                   >  2:495815/1‑150 (MQ=255)
    aaGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCAtggt                                                                                                                                                <  2:1085614/150‑1 (MQ=255)
     aGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCACCATGAtggcg                                                                                                                                               <  2:1518792/150‑3 (MQ=255)
      ggCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCAtggtgg                                                                                                                                              >  1:114361/1‑150 (MQ=255)
      ggCGGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGGGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCAtggtgg                                                                                                                                              <  1:1519899/150‑1 (MQ=255)
        cTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGAt                                                                                                                                            <  2:1510597/150‑1 (MQ=255)
          gggggCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCACCATGATGGCGGAttg                                                                                                                                          <  2:6666/150‑3 (MQ=255)
           ggggCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCACCATGATGGCGGAttgc                                                                                                                                         <  1:356958/150‑4 (MQ=255)
             ggCGCAGGTCGATTTTATTGGTCGCATAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATAc                                                                                                                                          <  2:974310/147‑1 (MQ=255)
             ggCGCAGGTCGATTTTATTGGTCGCATAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATAc                                                                                                                                          >  1:974310/1‑147 (MQ=255)
               cgcAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAg                                                                                                                                     >  1:976644/1‑150 (MQ=255)
                  aGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACGAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCg                                                                                                                                  <  2:1393917/150‑1 (MQ=255)
                      cGATTTTATTGGTCGCGTAGGTGATGACGAGACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACATCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCg                                                                                                                              <  2:1383777/150‑1 (MQ=255)
                      cGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCg                                                                                                                              <  1:495070/150‑1 (MQ=255)
                      cGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCg                                                                                                                              <  2:550120/150‑1 (MQ=255)
                       gATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAAACAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCgg                                                                                                                             >  2:163290/1‑150 (MQ=255)
                         ttttATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGAt                                                                                                                           <  2:494454/150‑1 (MQ=255)
                            tATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAAt                                                                                                                        <  2:1590188/150‑1 (MQ=255)
                                ggTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCAtt                                                                                                                    <  2:1503396/150‑1 (MQ=255)
                                 gTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCCGGCAGAGCTAGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTa                                                                                                                   <  2:1587089/150‑1 (MQ=255)
                                         ggTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTAcccc                                                                                                           >  2:627092/1‑150 (MQ=255)
                                             atgaCGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGcc                                                                                                       <  1:644849/150‑1 (MQ=255)
                                               gacgacACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCg                                                                                                     >  1:1573813/1‑150 (MQ=255)
                                                    cacGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATct                                                                                                <  2:1623952/150‑1 (MQ=255)
                                                    cacGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATct                                                                                                >  1:1138932/1‑150 (MQ=255)
                                                         gCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCtgctgc                                                                                           >  2:1254970/1‑150 (MQ=255)
                                                                 ctgctgGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCa                                                                                   <  2:1414915/150‑1 (MQ=255)
                                                                      ggCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTg                                                                              >  2:1133256/1‑150 (MQ=255)
                                                                       gCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTgg                                                                             >  2:78682/1‑150 (MQ=255)
                                                                           agCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGCCGCAGAGTGGTTgg                                                                         >  2:1433631/1‑150 (MQ=255)
                                                                             cTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTggag                                                                       <  1:929244/150‑1 (MQ=255)
                                                                                     cTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGa                                                               <  2:347239/150‑1 (MQ=255)
                                                                                      tGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGAt                                                              <  1:424010/150‑1 (MQ=255)
                                                                                       ggggCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGAtt                                                             >  1:101473/1‑150 (MQ=255)
                                                                                           cGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTtctc                                                         >  1:1034949/1‑150 (MQ=255)
                                                                                               aaCACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTa                                                                                                                   >  1:1444264/1‑88 (MQ=255)
                                                                                               aaCACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTa                                                                                                                   <  2:1444264/88‑1 (MQ=255)
                                                                                               aaCACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCCGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAg                                                     >  2:968227/1‑150 (MQ=255)
                                                                                               aaCACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAg                                                     >  2:412076/1‑150 (MQ=255)
                                                                                                acacCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGt                                                    <  2:1371986/150‑1 (MQ=255)
                                                                                                 cacCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTg                                                   >  2:1409089/1‑150 (MQ=255)
                                                                                                 cacCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTg                                                   <  2:475903/150‑1 (MQ=255)
                                                                                                        tACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGAtgtt                                            >  1:947857/1‑150 (MQ=255)
                                                                                                         acacCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGAtgttg                                           >  2:433318/1‑150 (MQ=255)
                                                                                                             cAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCt                                       <  2:1189445/150‑1 (MQ=255)
                                                                                                             cAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCt                                       <  2:964659/150‑1 (MQ=255)
                                                                                                               aaCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTgg                                     <  2:520006/150‑1 (MQ=255)
                                                                                                               aaCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTgg                                     >  2:649486/1‑150 (MQ=255)
                                                                                                                aCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGc                                    <  2:1316119/150‑1 (MQ=255)
                                                                                                                aCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGc                                    >  2:1436268/1‑150 (MQ=255)
                                                                                                                aCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGAGACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGc                                    <  1:751956/150‑1 (MQ=255)
                                                                                                                 cGGAATAACCAGGCGAAATCTTCGCAATCCGCCATCATGCTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCa                                   >  2:1193195/1‑150 (MQ=255)
                                                                                                                          cAGGCGAAATCTTCGCAACCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGc                          >  1:561343/1‑150 (MQ=255)
                                                                                                                            ggCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTg                        <  1:1160661/150‑1 (MQ=255)
                                                                                                                            ggCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTg                        >  1:811918/1‑150 (MQ=255)
                                                                                                                             gCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTgg                       <  1:1266654/150‑1 (MQ=255)
                                                                                                                             gCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTgg                       >  2:613634/1‑150 (MQ=255)
                                                                                                                              cGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTgcc                      >  2:946365/1‑148 (MQ=255)
                                                                                                                                    cTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCacgacg                <  1:1361931/150‑1 (MQ=255)
                                                                                                                                        gCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCa                                   >  2:1487069/1‑127 (MQ=255)
                                                                                                                                        gCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCa                                   <  1:1487069/127‑1 (MQ=255)
                                                                                                                                         gtaTCCACCATGATGGCGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCt           >  1:1496030/3‑150 (MQ=255)
                                                                                                                                         gtaTCCACCATGATGGCGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCt           >  2:1004132/3‑150 (MQ=255)
                                                                                                                                          aaTCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTa          <  2:1138932/150‑1 (MQ=255)
                                                                                                                                                  atcatGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGcc  <  2:1573813/150‑1 (MQ=255)
                                                                                                                                                  atcatGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGcc  <  2:1161751/150‑1 (MQ=255)
                                                                                                                                                  atcatGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGcc  <  2:947857/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCC  >  NC_000913/4073881‑4074176

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