Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,507,106 A→G 7.4% G263G (GGT→GGC fryC ← putative PTS enzyme IIC component FryC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,507,1060AG 7.4% 431.2 / 18.4 149G263G (GGT→GGCfryCputative PTS enzyme IIC component FryC
Reads supporting (aligned to +/- strand):  ref base A (67/70);  new base G (5/6);  total (73/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-01

Only 100 of 155 total aligned reads displayed.

GGATAAAGTTGAGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACGTA  >  NC_000913/2506957‑2507255
                                                                                                                                                     |                                                                                                                                                     
ggATAAAGTTGAGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGCCCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACaaaa                                                                                                                                                       >  2:1889194/1‑150 (MQ=255)
      aGTTGAGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCg                                                                                                                                                 <  1:392536/150‑1 (MQ=255)
         tGAGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCt                                                                                                                                              >  1:1023410/1‑150 (MQ=255)
         tGAGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCt                                                                                                                                              >  2:265446/1‑150 (MQ=255)
           aGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCttt                                                                                                                                            >  1:720056/1‑150 (MQ=255)
           aGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCttt                                                                                                                                            <  1:308858/150‑1 (MQ=255)
           aGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCttt                                                                                                                                            <  2:2348291/150‑1 (MQ=255)
            gCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTttg                                                                                                                                           <  2:152751/150‑1 (MQ=255)
            gCGTTGAACAGGCGTTTGCCGGTTAAACGACAGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTttg                                                                                                                                           <  2:1688786/150‑1 (MQ=255)
              gTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTttgtt                                                                                                                                         <  1:278123/150‑1 (MQ=255)
              gTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTttgtt                                                                                                                                         <  1:282343/150‑1 (MQ=255)
              gTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTttgtt                                                                                                                                         <  2:1931065/150‑1 (MQ=255)
               ttGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAGCCGCAGGTTttgttg                                                                                                                                        <  2:1186189/150‑1 (MQ=255)
                tGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGa                                                                                                                                       >  1:1954472/1‑150 (MQ=255)
                 gAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAGCCGCAGGTTttgttgcc                                                                                                                                      <  2:1238223/150‑3 (MQ=255)
                  aaCAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAGCCGCAGGTTttgttgcct                                                                                                                                     <  2:2047114/150‑4 (MQ=255)
                  aaCAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAGCCGCAGGTTttgttgcct                                                                                                                                     <  2:888439/150‑4 (MQ=255)
                     aGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCggg                                                                                                                                  <  1:526439/150‑1 (MQ=255)
                         gTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGccac                                                                                                                              >  1:2322418/1‑150 (MQ=255)
                          tttGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGccacc                                                                                                                             <  1:487591/150‑1 (MQ=255)
                                 gTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCg                                                                                                                      <  1:489081/150‑1 (MQ=255)
                                     aaCGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTg                                                                                                                  <  2:1023410/150‑1 (MQ=255)
                                      acgacgGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGc                                                                                                                 <  1:2387093/150‑1 (MQ=255)
                                       cgacgGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCt                                                                                                                >  2:179278/1‑150 (MQ=255)
                                        gacgGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTg                                                                                                               <  1:376223/150‑1 (MQ=255)
                                             tCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTcgcc                                                                                                          <  1:1540221/150‑1 (MQ=255)
                                                 tAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCgcgg                                                                                                      <  1:17764/150‑1 (MQ=255)
                                                  aaTGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCgcggc                                                                                                     <  2:239480/150‑1 (MQ=255)
                                                      gTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATc                                                                                                 >  2:2088687/1‑150 (MQ=255)
                                                      gTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATc                                                                                                 >  2:166481/1‑150 (MQ=255)
                                                         gCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATcccc                                                                                              >  1:1450506/1‑150 (MQ=255)
                                                          ccagccccagccCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCa                                                                                             >  1:331316/1‑150 (MQ=255)
                                                             gCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTg                                                                                          <  2:1608416/150‑1 (MQ=255)
                                                              ccccAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGCTAATACGAGGTCAGTGGTAAACCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTAGCCGCGGCGATACCCATTGc                                                                                         >  2:129103/1‑150 (MQ=255)
                                                                   gACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACa                                                                                    <  2:1928215/150‑1 (MQ=255)
                                                                      cAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCa                                                                                 >  2:142302/1‑150 (MQ=255)
                                                                      cAATCGGCGGAATAACGATAGCGACAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCa                                                                                 <  1:967679/150‑1 (MQ=255)
                                                                       aaTCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAg                                                                                <  2:2441313/150‑1 (MQ=255)
                                                                           ggcggAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAgcgcg                                                                            >  1:1344276/1‑150 (MQ=255)
                                                                            gcggAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAgcgcgc                                                                           <  2:1129877/150‑1 (MQ=255)
                                                                                 aTAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTtc                                                                      <  2:789523/150‑1 (MQ=255)
                                                                                 aTAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGGCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTtc                                                                      <  1:2050955/150‑1 (MQ=255)
                                                                                   aaCGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTtctc                                                                    >  2:726444/1‑150 (MQ=255)
                                                                                   aaCGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTtctc                                                                    <  1:158195/150‑1 (MQ=255)
                                                                                    aCGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAACCCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCAGTCGCCGCGGCGAGCCCCATTGCGTCCCTCAGCGCGCCCTTCTCa                                                                   >  2:950714/1‑150 (MQ=255)
                                                                                    aCGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCa                                                                   <  2:297952/150‑1 (MQ=255)
                                                                                    aCGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCa                                                                   >  1:1652349/1‑150 (MQ=255)
                                                                                           cgatagAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCgg                                                            >  2:507489/1‑150 (MQ=255)
                                                                                            gatagAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGc                                                           <  2:1652525/150‑1 (MQ=255)
                                                                                            gatagAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGc                                                           <  1:1688331/150‑1 (MQ=255)
                                                                                            gatagAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGc                                                           >  1:2356047/1‑150 (MQ=255)
                                                                                            gatagAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGc                                                           <  1:166481/150‑1 (MQ=255)
                                                                                            gatagAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTAACCTGCGGc                                                           >  2:906211/1‑150 (MQ=255)
                                                                                              tagAACGAGCTGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAg                                                         <  1:1253243/150‑1 (MQ=255)
                                                                                              tagAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAg                                                         <  2:619189/150‑1 (MQ=255)
                                                                                                taaCGAGCGGTGACCTGTAATTCGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCTCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTTCGGCAGTc                                                       <  2:1954472/149‑1 (MQ=255)
                                                                                                gAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTc                                                       <  2:571105/150‑1 (MQ=255)
                                                                                                gAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTc                                                       >  2:1611867/1‑150 (MQ=255)
                                                                                                gAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTc                                                       <  2:2419327/150‑1 (MQ=255)
                                                                                                   cGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGc                                                    >  1:884303/1‑150 (MQ=255)
                                                                                                     aGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGcac                                                  <  1:146101/150‑1 (MQ=255)
                                                                                                     aGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGcac                                                  <  1:1409372/150‑1 (MQ=255)
                                                                                                        ggTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGt                                               >  2:439307/1‑150 (MQ=255)
                                                                                                          tGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTAc                                             <  1:644277/150‑1 (MQ=255)
                                                                                                            aCCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACgg                                           <  1:439307/150‑1 (MQ=255)
                                                                                                               ggTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCa                                                                   >  1:1535853/1‑123 (MQ=255)
                                                                                                               ggTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCa                                                                   <  2:1535853/123‑1 (MQ=255)
                                                                                                               ggTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATa                                        >  1:2164577/1‑150 (MQ=255)
                                                                                                                    tACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGACGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATAccgcc                                   >  1:336571/1‑150 (MQ=255)
                                                                                                                      cGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGt                                 >  2:426037/1‑150 (MQ=255)
                                                                                                                       gTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCACGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGtt                                <  1:663063/150‑1 (MQ=255)
                                                                                                                        tGGTCAGGGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTg                               <  2:2248097/150‑1 (MQ=255)
                                                                                                                         ggTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGa                              >  1:1292059/1‑150 (MQ=255)
                                                                                                                          gTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGAt                             <  2:2164577/150‑1 (MQ=255)
                                                                                                                            cAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATcc                           <  1:2372403/150‑1 (MQ=255)
                                                                                                                              gTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAg                         <  1:310527/150‑1 (MQ=255)
                                                                                                                                ggTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGcc                       >  1:1769591/1‑150 (MQ=255)
                                                                                                                                 gTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGcca                      >  2:1273943/1‑150 (MQ=255)
                                                                                                                                  tAAAGCTGAAGGCAACAAAACCTGCGGCGTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGccac                     >  2:1097862/1‑150 (MQ=255)
                                                                                                                                   aaaGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGccacc                    <  2:331316/150‑1 (MQ=255)
                                                                                                                                      gCTGAAGGCAACAAAACCTGCGGCTTTGTTGATAGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCaaa                 >  1:2463287/1‑150 (MQ=255)
                                                                                                                                       cagAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAg                >  2:1166905/3‑150 (MQ=255)
                                                                                                                                        tGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAgg               <  2:1505337/150‑1 (MQ=255)
                                                                                                                                        cccgatcaacaaAGCCGCAGGTTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAgg               >  1:823390/7‑150 (MQ=255)
                                                                                                                                         ccgatcaacaaAGCCGCAGGTTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAggg              >  2:2071323/6‑150 (MQ=255)
                                                                                                                                          aaGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAgggg             <  1:1866535/150‑1 (MQ=255)
                                                                                                                                          aaGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAgggg             >  1:1397057/1‑150 (MQ=255)
                                                                                                                                           gatcaacaaAGCCGCAGGTTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAggggg            >  1:2068497/4‑150 (MQ=255)
                                                                                                                                           gatcaacaaAGCCGCAGGTTTTGTTGATCGGGACACCAAGGTCGATTGATGTCGCCGCGGCGATACCCATTGCGTACATCAGCGAGCCCTTATCAAATTCGGCAGAAAGCACAGTACGGATACCGCCGTTGATCCAGCCCCCAAAggggg            >  2:1560661/4‑150 (MQ=255)
                                                                                                                                           aGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAggggg            >  1:1100206/1‑150 (MQ=255)
                                                                                                                                           aGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAggggg            >  1:2342889/1‑150 (MQ=255)
                                                                                                                                             tcaacaaAGCCGCAGGTTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGAAACCGCCGTTGATCCAGCCACCAAAGGCGGtg          >  2:1977185/2‑150 (MQ=255)
                                                                                                                                              caacaaAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGtga         >  1:541793/1‑150 (MQ=255)
                                                                                                                                               aacaaAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGtgat        >  2:1244010/1‑150 (MQ=255)
                                                                                                                                                acaaAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGtgatg       >  2:337521/1‑150 (MQ=255)
                                                                                                                                                acaaAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGtgatg       >  1:1733461/1‑150 (MQ=255)
                                                                                                                                                 caaAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGtgctga      >  2:1129658/1‑150 (MQ=255)
                                                                                                                                                 caaAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGtgatga      <  1:1203892/150‑1 (MQ=255)
                                                                                                                                                    aaCCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACg  <  1:715941/150‑2 (MQ=255)
                                                                                                                                                     aCCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACgta  >  1:196953/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGATAAAGTTGAGCGTTGAACAGGCGTTTGCCGGTTAAACGACGGTCGATAATGGTCGCCAGCCCCAGACCAATCGGCGGAATAACGATAGCGATAGAACGAGCGGTGACCGGTAATACGTGGTCAGTGGTAAAGCTGAAGGCAACAAAACCTGCGGCTTTGTTGATCGGGCCACCAAGGTCGATTGCTGTCGCCGCGGCGATCCCCATTGCGTACATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACGTA  >  NC_000913/2506957‑2507255

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: