Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,921,420 C→A 7.9% M227I (ATG→ATT)  gudP ← galactarate/glucarate/glycerate transporter GudP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,921,4200CA 7.9% 201.6 / 5.0 63M227I (ATG→ATT) gudPgalactarate/glucarate/glycerate transporter GudP
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base A (2/3);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

ACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAAAGAAGAAGACGTGTGACCAGC  >  NC_000913/2921271‑2921568
                                                                                                                                                     |                                                                                                                                                    
aCCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATcc                                                                                                                                                      <  2:1795008/150‑1 (MQ=255)
    ggAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCatat                                                                                                                                                  <  1:1157177/150‑1 (MQ=255)
     gAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATAtt                                                                                                                                                 <  1:381510/150‑1 (MQ=255)
      aaCCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGAGATATATAAACGCCGATCATCATCCGTGACCCTAGGAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTg                                                                                                                                                <  2:942499/150‑1 (MQ=255)
        ccAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGAt                                                                                                                                              >  1:1081528/1‑150 (MQ=255)
            gTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAgc                                                                                                                                          <  1:924093/150‑1 (MQ=255)
             taataaAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAgcg                                                                                                                                         >  1:939471/1‑150 (MQ=255)
               ataaAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCAATATGGATCAgcaaa                                                                                                                                       <  2:460101/150‑4 (MQ=255)
                 aaaGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGcacc                                                                                                                                     >  2:963479/1‑150 (MQ=255)
                  aagaagTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGcacca                                                                                                                                    >  1:711104/1‑150 (MQ=255)
                   agaagTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAACGcaccac                                                                                                                                   >  1:1673474/1‑150 (MQ=255)
                    gaagTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCaccacc                                                                                                                                  >  1:1638343/1‑150 (MQ=255)
                    gaagTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAACGCaccacc                                                                                                                                  <  2:1673474/150‑1 (MQ=255)
                         aagtCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCgcgg                                                                                                                             >  1:772931/1‑150 (MQ=255)
                             cAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGAt                                                                                                                         >  2:176503/1‑150 (MQ=255)
                               gggCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGt                                                                                                                       <  2:840861/150‑1 (MQ=255)
                                ggCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTa                                                                                                                      >  1:1393640/1‑150 (MQ=255)
                                ggCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTa                                                                                                                      >  1:1137572/1‑150 (MQ=255)
                                    ttGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTcc                                                                                                                  <  1:279219/150‑1 (MQ=255)
                                     tGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGATTATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCa                                                                                                                 >  2:1828882/1‑150 (MQ=255)
                                              aCTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCtctttct                                                                                                        <  1:438514/150‑1 (MQ=255)
                                              aCTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCtctttct                                                                                                        <  1:394882/150‑1 (MQ=255)
                                                  aCCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTAtt                                                                                                    >  2:361065/1‑150 (MQ=255)
                                                             cGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGAt                                                                                         >  2:1429939/1‑150 (MQ=255)
                                                                cGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGtt                                                                                      <  1:1191623/150‑1 (MQ=255)
                                                                cGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGtt                                                                                      <  1:963479/150‑1 (MQ=255)
                                                                cGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGtt                                                                                      <  2:1638343/150‑1 (MQ=255)
                                                                cGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGtt                                                                                      >  1:1010643/1‑150 (MQ=255)
                                                                cGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGtt                                                                                      <  1:789484/150‑1 (MQ=255)
                                                                 gATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTA‑TTACCCCCGGATGTTGa                                                                                    <  1:1136039/150‑1 (MQ=255)
                                                                    catcatcCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGAtt                                                                                  <  1:12530/150‑1 (MQ=255)
                                                                      tcatcCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTg                                                                                <  2:77531/150‑1 (MQ=255)
                                                                       catcCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTA‑CCCCGGATGTTGATTTGGc                                                                              >  1:965967/1‑150 (MQ=255)
                                                                          ccGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCAATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTc                                                                            >  2:362880/1‑150 (MQ=255)
                                                                           cGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCg                                                                           <  1:1564642/150‑1 (MQ=255)
                                                                           cGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCg                                                                           >  1:1685679/1‑150 (MQ=255)
                                                                               ccccTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCCGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGa                                                                       >  1:1843054/2‑150 (MQ=255)
                                                                                 ccTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATg                                                                     <  2:1081528/150‑1 (MQ=255)
                                                                                 ccTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATg                                                                     >  2:1297577/1‑150 (MQ=255)
                                                                                      cagcTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTcc                                                                                                                  <  1:1163978/100‑1 (MQ=255)
                                                                                      cagcTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTcc                                                                                                                  >  2:1163978/1‑100 (MQ=255)
                                                                                           gTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAAcc                                                           <  2:1010643/150‑1 (MQ=255)
                                                                                             ttGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAg                                                         <  2:1137572/150‑1 (MQ=255)
                                                                                              tGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGa                                                        <  2:669122/150‑1 (MQ=255)
                                                                                               gATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGAt                                                       >  1:1333730/1‑150 (MQ=255)
                                                                                                    gCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGc                                                  >  2:878113/1‑150 (MQ=255)
                                                                                                    gCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGc                                                  <  1:1429939/150‑1 (MQ=255)
                                                                                                       cccACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCtga                                               >  1:1816971/1‑150 (MQ=255)
                                                                                                         cACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCtgatg                                             <  2:1142688/150‑1 (MQ=255)
                                                                                                          aCTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCtgatga                                            <  2:1816971/150‑1 (MQ=255)
                                                                                                           cTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGAc                                           >  2:1203289/1‑150 (MQ=255)
                                                                                                             tCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACaa                                         <  1:361065/150‑1 (MQ=255)
                                                                                                                         gCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGc                             <  2:772931/150‑1 (MQ=255)
                                                                                                                          cACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGcc                            <  2:939471/150‑1 (MQ=255)
                                                                                                                           aCTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGccc                           <  1:1203289/150‑1 (MQ=255)
                                                                                                                                 aCTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAaag                     >  2:84760/1‑150 (MQ=255)
                                                                                                                                  cTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAaaga                    >  1:565498/1‑150 (MQ=255)
                                                                                                                                      ggTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAaagaagaa                <  1:176503/150‑1 (MQ=255)
                                                                                                                                          gtgcgCTGATCAATATGGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATa                                                      >  1:167219/5‑108 (MQ=255)
                                                                                                                                          gtgcgCTGATCAATATGGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATa                                                      <  2:167219/104‑1 (MQ=255)
                                                                                                                                            ttGCTGATCCATATTGATCAGCGCACCACCAGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAAAGAAGAAGACgtg          >  2:814008/1‑150 (MQ=255)
                                                                                                                                            ttGCTGATCAATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCAAGACCGCCCATAAAGAAGAAGACgtg          <  1:362880/150‑1 (MQ=255)
                                                                                                                                               cTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAAAGAAGAAGACGTGTGa       >  2:1029628/1‑150 (MQ=255)
                                                                                                                                                 gATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAAAGAAGAAGACGTGTGAcc     >  1:1433080/1‑150 (MQ=255)
                                                                                                                                                    ccATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAAAGAAGAAGACGTGTGACCAGc  >  2:475703/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGTGTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAAAGAAGAAGACGTGTGACCAGC  >  NC_000913/2921271‑2921568

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