Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,507,114 C→T 5.4% A261T (GCC→ACC)  fryC ← putative enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,507,1140CT 5.4% 416.8 / 7.7 129A261T (GCC→ACC) fryCputative enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base C (57/65);  new base T (3/4);  total (60/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.77e-01

Only 100 of 130 total aligned reads displayed.

TGAGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACGTAGTA  >  NC_000913/2506966‑2507258
                                                                                                                                                    |                                                                                                                                                
tGAGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACATGCGGCt                                                                                                                                                 >  2:333961/1‑150 (MQ=255)
  aGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCttt                                                                                                                                               >  1:1275985/1‑150 (MQ=255)
   gCGTTGAACAGGCGTTTGCCGGTTAAACGACGGGCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAAAAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTttg                                                                                                                                              <  2:1631978/150‑1 (MQ=255)
     gTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGCACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTttgtt                                                                                                                                            >  2:2109507/1‑150 (MQ=255)
     gTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGGCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTttgtt                                                                                                                                            <  1:1941838/150‑1 (MQ=255)
      ttGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTttgttg                                                                                                                                           <  2:1880384/150‑1 (MQ=255)
      ttGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTttgttg                                                                                                                                           <  1:718666/150‑1 (MQ=255)
       tGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGa                                                                                                                                          <  1:2085914/150‑1 (MQ=255)
        gAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATACCGCTAGCGATAGAACGAGCGGGGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTGTGTTGAt                                                                                                                                         >  1:309619/1‑150 (MQ=255)
         aaCAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATATCGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTCTTGaac                                                                                                                                        >  1:2050929/1‑148 (MQ=255)
         aaCAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAGCCGCAGGTTttgttgcct                                                                                                                                        <  2:2065715/150‑4 (MQ=255)
         aaCAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATc                                                                                                                                        <  2:908618/150‑1 (MQ=255)
           cAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCgg                                                                                                                                      >  2:1763444/1‑150 (MQ=255)
             ggCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTAGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGc                                                                                                                                    >  1:1342548/1‑150 (MQ=255)
             ggCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTAAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGc                                                                                                                                    >  1:209760/1‑150 (MQ=255)
                  ttGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGccacca                                                                                                                               <  1:752628/150‑1 (MQ=255)
                    gCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAg                                                                                                                             >  2:1038772/1‑150 (MQ=255)
                     ccGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAgg                                                                                                                            >  2:929156/1‑150 (MQ=255)
                           aaaCGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGAtt                                                                                                                      >  1:1523017/1‑150 (MQ=255)
                            aaCGACGGTCGATTATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGGGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTTAATCTGAAGTCAACATAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTg                                                                                                                     <  2:1945839/150‑1 (MQ=255)
                            aaCGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTg                                                                                                                     >  2:1561437/1‑150 (MQ=255)
                            aaCGACGGTCGATAATGGTCGCCAGCCACAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTg                                                                                                                     >  1:551563/1‑150 (MQ=255)
                             acgacgGTCGATAATGGTCGCCAGCCCCAGACGAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTTGTCTGTGGTAAAGCTGAAGGCAACATTACCTGCGGCTTGGTTGATCGGGCCACCAAGGTCGATTGc                                                                                                                    <  2:306394/150‑1 (MQ=255)
                             acgacgGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGc                                                                                                                    <  1:2201552/150‑1 (MQ=255)
                             acgacgGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGc                                                                                                                    <  1:1616806/150‑1 (MQ=255)
                               gacgGTCGCTAATGGGCGCCAGCCCCAGACCAATCGGCGGCATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTg                                                                                                                  <  1:1393120/150‑1 (MQ=255)
                               gacgGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTg                                                                                                                  <  2:2086273/150‑1 (MQ=255)
                                acgGTCGATAATGGTAGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGt                                                                                                                 >  1:134081/1‑150 (MQ=255)
                                  ggTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATACCGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTcg                                                                                                               <  2:85218/150‑1 (MQ=255)
                                   gTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGAAAGCGATAGAACGAGCGGTGACCGGTAATACGTTGTCAGTGGTAAAGCTGAAGGCAACAAAACCTTCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTcgc                                                                                                              <  1:1317562/150‑1 (MQ=255)
                                       aTAATGGTCTCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCcgcg                                                                                                          <  2:1414094/150‑1 (MQ=255)
                                         aaTGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCgcggc                                                                                                        <  2:1751879/150‑1 (MQ=255)
                                          aTGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGCACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCgcggcg                                                                                                       >  1:1485862/1‑150 (MQ=255)
                                          aTGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCgcggcg                                                                                                       <  1:1524055/150‑1 (MQ=255)
                                          aTGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCgcggcg                                                                                                       >  1:662915/1‑150 (MQ=255)
                                              tCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGGCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAACGGCGATTGCTGTCGCCGCGGCGATcc                                                                                                   >  1:384868/1‑150 (MQ=255)
                                              gcGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATACAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATcc                                                                                                   <  1:1880164/149‑1 (MQ=255)
                                                 ccAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGCAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCGTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCa                                                                                                >  1:469206/1‑150 (MQ=255)
                                                      cccAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCg                                                                                           >  2:745905/1‑150 (MQ=255)
                                                       ccAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTGTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGCTCCCCATTGCGt                                                                                          >  1:1974858/1‑150 (MQ=255)
                                                        cAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTa                                                                                         >  2:1516380/1‑150 (MQ=255)
                                                         aGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTAc                                                                                        <  2:1015812/150‑1 (MQ=255)
                                                         aGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTAc                                                                                        <  2:2162235/150‑1 (MQ=255)
                                                         aGAACAATCGGCGGAATAACGATAGCGATAGAACGAGCGGAGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGGCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTAc                                                                                        <  2:1392507/150‑1 (MQ=255)
                                                          gcccAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACa                                                                                       <  1:1375038/148‑1 (MQ=255)
                                                            ccAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATc                                                                                     >  2:64427/1‑150 (MQ=255)
                                                              aaTCGGCGGAATTACGATAGCGATATAACGAGCGGTGACCGGTAATTCGTGGTCAGTGGTAAAGCTGTAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAg                                                                                   <  1:857495/150‑1 (MQ=255)
                                                                   gcggAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAgcgcgc                                                                              >  1:1892394/1‑150 (MQ=255)
                                                                    cggAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGcc                                                                             >  1:1688632/1‑150 (MQ=255)
                                                                     ggAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGccc                                                                            <  1:859769/150‑1 (MQ=255)
                                                                         tAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTtct                                                                        <  2:134081/150‑1 (MQ=255)
                                                                         tAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTtct                                                                        >  1:829941/1‑150 (MQ=255)
                                                                         aaaCGATAGCGATAGAACGAGCGGTGACCGGTAATATGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTtct                                                                        <  2:1636775/149‑1 (MQ=255)
                                                                          aaCGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTtctc                                                                       >  2:1039748/1‑150 (MQ=255)
                                                                           aCGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCCCCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCa                                                                      >  2:592225/1‑150 (MQ=255)
                                                                            cgatagcgatagAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCAc                                                                     >  1:1946338/1‑150 (MQ=255)
                                                                             gatagcgatagAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCAcc                                                                    <  1:226353/150‑1 (MQ=255)
                                                                             gatagcgatagAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCAcc                                                                    <  2:2119328/150‑1 (MQ=255)
                                                                                agcgatagAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGc                                                                 <  2:2102876/150‑1 (MQ=255)
                                                                                  cgatagAACGAGCGGGGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCGTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCgg                                                               <  1:927246/150‑1 (MQ=255)
                                                                                   gatagAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGc                                                              <  1:1038772/150‑1 (MQ=255)
                                                                                    atagAACGAGCGGTGACCGGTATTACGTGGTCAGTGGTAACGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCa                                                             <  1:499189/150‑1 (MQ=255)
                                                                                    atagAACGAGCGGTGACCGGTAATACGTGGGCAGTGGGAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCa                                                             <  2:924403/150‑1 (MQ=255)
                                                                                      agAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGt                                                           <  2:2034212/150‑1 (MQ=255)
                                                                                       gAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTc                                                          <  2:1548827/150‑1 (MQ=255)
                                                                                        aaCGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCa                                                         <  2:20539/150‑1 (MQ=255)
                                                                                          cGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGc                                                       >  1:497500/1‑150 (MQ=255)
                                                                                                gTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTa                                                 >  1:2149379/1‑150 (MQ=255)
                                                                                                  gACCGGTAATACGTGGTCAGTGGTAAAGCTGACGGCAACAAAACCTGCGGCGTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACg                                               <  1:1098192/150‑1 (MQ=255)
                                                                                                  gACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACg                                               >  1:462514/1‑150 (MQ=255)
                                                                                                      ggTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATa                                           >  1:1735340/1‑150 (MQ=255)
                                                                                                        tAATACGTGGTCAGTGGTAAAGATGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGGCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATAcc                                         <  1:1387462/150‑1 (MQ=255)
                                                                                                             cGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGt                                    >  2:868006/1‑150 (MQ=255)
                                                                                                             cGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGt                                    >  2:1149772/1‑150 (MQ=255)
                                                                                                               tcgTCAGTGGTACAGCTGAAGGCCACACAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCAATGCGTCCATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTg                                  <  2:1374195/148‑1 (MQ=255)
                                                                                                                   cAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGGCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATcc                              <  2:1473932/150‑1 (MQ=255)
                                                                                                                       ggTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGcc                          <  1:2067598/150‑1 (MQ=255)
                                                                                                                         tAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGccac                        >  1:2075972/1‑150 (MQ=255)
                                                                                                                          aaaGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGccacc                       <  1:1516583/150‑1 (MQ=255)
                                                                                                                          aaaGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGccacc                       >  1:721995/1‑150 (MQ=255)
                                                                                                                           aaGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGccacca                      <  1:2108432/150‑1 (MQ=255)
                                                                                                                             gCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCaaa                    >  2:974359/1‑150 (MQ=255)
                                                                                                                               cccgatcaacaaAGCCGCAGGTTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAgg                  >  2:2203344/7‑150 (MQ=255)
                                                                                                                                 cgatcaacaaAGCCGCAGGTTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCCCCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAgggg                >  2:1517423/5‑150 (MQ=255)
                                                                                                                                    gCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGtg             >  2:1894582/1‑150 (MQ=255)
                                                                                                                                     caacaaAGCCGCAGGTTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGACCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGtga            >  2:1695391/1‑150 (MQ=255)
                                                                                                                                     caacaaAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATGCAGCCACCAAAGGGGGtga            >  2:293694/1‑150 (MQ=255)
                                                                                                                                     caacaaAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACAGCCGTTGATACAGCCACCAAAGGGGGtga            >  1:592065/1‑150 (MQ=255)
                                                                                                                                     caacaaAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGGCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCCGCACAGTACGGATACCGCCGTTGCTCCAGCCACCAAAGGGGGtga            >  1:406553/1‑150 (MQ=255)
                                                                                                                                      aacaaAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGtgat           >  2:1733896/1‑150 (MQ=255)
                                                                                                                                       acaaAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGtgatg          <  1:1578029/150‑1 (MQ=255)
                                                                                                                                         aaaaCCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCACCCACCAAAGGGGGggatgcc        >  1:1124023/1‑148 (MQ=255)
                                                                                                                                          aaaCCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGCTCCAGCCACCAAAGGGGGTGATGACg       >  1:1148915/1‑150 (MQ=255)
                                                                                                                                          aaaCCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACg       <  1:1763444/150‑1 (MQ=255)
                                                                                                                                          aaaCCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACg       >  1:1173096/1‑150 (MQ=255)
                                                                                                                                           acccTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCAGGCGCCGCGGCGATCGCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACgt      <  1:333961/148‑1 (MQ=255)
                                                                                                                                              cTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACgtagt   >  2:134152/1‑150 (MQ=255)
                                                                                                                                               tGCGGCTTTGTTGATCGGTCCACCAAGGTCTATTGCTGGCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACgtagta  <  2:2078417/150‑1 (MQ=255)
                                                                                                                                               tGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACgtagta  <  2:686444/150‑1 (MQ=255)
                                                                                                                                               ggcggcgTTGTTGATCGGGCCACCAAGGTCGATTGCTGGCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACgtagta  <  1:1473547/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TGAGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACGTAGTA  >  NC_000913/2506966‑2507258

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