Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,507,106 A→G 5.6% G263G (GGT→GGC 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,1060AG 5.6% 373.3 / 9.3 124G263G (GGT→GGCfryCputative enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base A (55/62);  new base G (3/4);  total (58/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.68e-01

Only 100 of 130 total aligned reads displayed.

ATAAAGTTGAGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACG  >  NC_000913/2506959‑2507253
                                                                                                                                                   |                                                                                                                                                   
aTAAAGTTGAGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAAcc                                                                                                                                                   <  2:2206498/150‑1 (MQ=255)
   aaGTTGAGAGTTGAACAGGCGTTTGACTGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAATCTCAAGGCAACAAAACCTGc                                                                                                                                                >  1:93701/1‑150 (MQ=255)
    aGTTGAGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCg                                                                                                                                               >  2:226353/1‑150 (MQ=255)
     gTTGAGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGGCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCgg                                                                                                                                              <  2:144961/150‑1 (MQ=255)
         aGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCttt                                                                                                                                          >  1:1275985/1‑150 (MQ=255)
          gCGTTGAACAGGCGTTTGCCGGTTAAACGACGGGCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAAAAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTttg                                                                                                                                         <  2:1631978/150‑1 (MQ=255)
           cGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCTATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATTGAGATTGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTttgt                                                                                                                                        <  1:525648/150‑1 (MQ=255)
            gTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGCACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTttgtt                                                                                                                                       >  2:2109507/1‑150 (MQ=255)
            gTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAGCCGCAGGTTttgtt                                                                                                                                       <  2:28162/150‑1 (MQ=255)
             ttGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTttgttg                                                                                                                                      <  1:718666/150‑1 (MQ=255)
              tGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAGCCGCAGGTTttgttgc                                                                                                                                     <  1:962689/150‑2 (MQ=255)
               gAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATACCGCTAGCGATAGAACGAGCGGGGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTGTGTTGAt                                                                                                                                    >  1:309619/1‑150 (MQ=255)
               gAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCTGCTTTGTTCAt                                                                                                                                    >  1:2078417/1‑150 (MQ=255)
                aaCAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAGCCGCAGGTTttgttgcct                                                                                                                                   <  2:2065715/150‑4 (MQ=255)
                aaCAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATc                                                                                                                                   <  2:908618/150‑1 (MQ=255)
                aaCAGGCGTTTGCCGCTTAAACGACGGTCGATAATGGGCGCCAGCCCCAGACCAATCGGCGGAATAACCAAAGCGATAGAACGAGCGTTGACTGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAGCCGCAGGTTttgttgcct                                                                                                                                   <  2:366934/150‑4 (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)
                                   aaCGACGGTCGATAATGGTCGCCAGCCACAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTg                                                                                                                >  1:551563/1‑150 (MQ=255)
                                    acgacgGTCGATAATGGTCGCCAGCCCCAGACGAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTTGTCTGTGGTAAAGCTGAAGGCAACATTACCTGCGGCTTGGTTGATCGGGCCACCAAGGTCGATTGc                                                                                                               <  2:306394/150‑1 (MQ=255)
                                    acgacgGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGc                                                                                                               <  1:1616806/150‑1 (MQ=255)
                                      gacgGTCGCTAATGGGCGCCAGCCCCAGACCAATCGGCGGCATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTg                                                                                                             <  1:1393120/150‑1 (MQ=255)
                                       acgGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGt                                                                                                            >  1:2160790/1‑150 (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:1668660/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)
                                                 aTGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCgcggcg                                                                                                  >  1:2029983/1‑150 (MQ=255)
                                                     gcGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATACAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATcc                                                                                              <  1:1880164/149‑1 (MQ=255)
                                                         cAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCaa                                                                                          >  2:1839098/1‑149 (MQ=255)
                                                         cAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTAGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCAt                                                                                          <  2:1342548/150‑1 (MQ=255)
                                                           gCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGGCAGTGGTAAATCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTg                                                                                        <  2:1509641/150‑1 (MQ=255)
                                                             cccAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCg                                                                                      >  2:745905/1‑150 (MQ=255)
                                                              ccAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTGTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGCTCCCCATTGCGt                                                                                     >  1:1974858/1‑150 (MQ=255)
                                                                 gcccAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACa                                                                                  <  1:1375038/148‑1 (MQ=255)
                                                                   ccAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATc                                                                                >  2:64427/1‑150 (MQ=255)
                                                                          gcggAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAgcgcgc                                                                         >  1:1892394/1‑150 (MQ=255)
                                                                           cggAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGcc                                                                        >  1:1688632/1‑150 (MQ=255)
                                                                            ggAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGccc                                                                       <  1:859769/150‑1 (MQ=255)
                                                                                tAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTtct                                                                   >  1:829941/1‑150 (MQ=255)
                                                                                tAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTtct                                                                   <  2:134081/150‑1 (MQ=255)
                                                                                aaaCGATAGCGATAGAACGAGCGGTGACCGGTAATATGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTtct                                                                   <  2:1636775/149‑1 (MQ=255)
                                                                                 aaCGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTtctc                                                                  <  1:1561437/150‑1 (MQ=255)
                                                                                 aaCGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTtctc                                                                  <  1:1506968/150‑1 (MQ=255)
                                                                                 aaCGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTtctc                                                                  >  2:1039748/1‑150 (MQ=255)
                                                                                  aCGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCCCCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCa                                                                 >  2:592225/1‑150 (MQ=255)
                                                                                  aCGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCa                                                                 <  1:1731800/150‑1 (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)
                                                                                        gcgatagAACGAGCGGGGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCg                                                           <  1:700174/150‑1 (MQ=255)
                                                                                         cgatagAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCgg                                                          <  2:2113577/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)
                                                                                              gAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTc                                                     <  2:1548827/150‑1 (MQ=255)
                                                                                              gAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTc                                                     <  2:2076261/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)
                                                                                                               tAATACGTGGTCAGTGGTAAAGATGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGGCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATAcc                                    <  1:1387462/150‑1 (MQ=255)
                                                                                                                aaTACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATAccg                                   >  1:779864/1‑150 (MQ=255)
                                                                                                                    cGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGt                               >  2:868006/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)
                                                                                                                                 aaaGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGTCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGccacc                  <  2:1556042/150‑1 (MQ=255)
                                                                                                                                  aaGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCCGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGccacca                 >  1:33008/1‑150 (MQ=255)
                                                                                                                                  aaGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGccacca                 <  1:2108432/150‑1 (MQ=255)
                                                                                                                                    gCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCaaa               >  2:974359/1‑150 (MQ=255)
                                                                                                                                    gCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCaaa               >  1:1916505/1‑150 (MQ=255)
                                                                                                                                      tGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAgg             >  1:431077/1‑150 (MQ=255)
                                                                                                                                      cccgatcaacaaAGCCGCAGGTTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAgg             >  2:2203344/7‑150 (MQ=255)
                                                                                                                                        cgatcaacaaAGCCGCAGGTTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCCCCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAgggg           >  2:1517423/5‑150 (MQ=255)
                                                                                                                                           gCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGtg        <  2:2029983/150‑1 (MQ=255)
                                                                                                                                           gCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGtg        >  2:1894582/1‑150 (MQ=255)
                                                                                                                                            caacaaAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATGCAGCCACCAAAGGGGGtga       >  2:293694/1‑150 (MQ=255)
                                                                                                                                            caacaaAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGGCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCCGCACAGTACGGATACCGCCGTTGCTCCAGCCACCAAAGGGGGtga       >  1:406553/1‑150 (MQ=255)
                                                                                                                                             aacaaAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGtgat      >  2:1733896/1‑150 (MQ=255)
                                                                                                                                                aaaaCCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGAc   <  2:1097785/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)
                                                                                                                                                   |                                                                                                                                                   
ATAAAGTTGAGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACG  >  NC_000913/2506959‑2507253

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