Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,968,163 Δ1 bp coding (272/2247 nt) ptsP ← PEP‑protein phosphotransferase enzyme I; GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,968,1600T.99.0% 464.1 / ‑3.4 97coding (275/2247 nt)ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (48/48);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCG  >  NC_000913/2968011‑2968306
                                                                                                                                                     |                                                                                                                                                  
aggaTTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCtt                                                                                                                                                   >  2:833486/1‑151 (MQ=255)
     tcttcACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCtttgcg                                                                                                                                               <  2:1367179/150‑4 (MQ=255)
      cttcACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCtttgcgc                                                                                                                                              >  2:1368510/1‑146 (MQ=255)
          aCTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTg                                                                                                                                        <  2:659269/151‑1 (MQ=255)
           cTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTg                                                                                                                                        <  2:131938/150‑1 (MQ=255)
            tttCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGc                                                                                                                                       >  1:29702/1‑150 (MQ=255)
              tCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCaa                                                                                                                                     <  2:917528/150‑1 (MQ=255)
              tCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCa                                                                                                                                      >  1:87294/1‑149 (MQ=255)
                 tCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGtt                                                                                                                                 >  1:1207600/1‑151 (MQ=255)
                  cATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGtt                                                                                                                                 <  2:426994/150‑1 (MQ=255)
                  cATACTGGCGCAACTCTCGTTGCTGGACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGtt                                                                                                                                 <  2:577960/150‑1 (MQ=255)
                   aTACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGttt                                                                                                                                >  2:165217/1‑150 (MQ=255)
                   aTACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGttt                                                                                                                                <  1:52036/150‑1 (MQ=255)
                   aTACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGttt                                                                                                                                >  2:561781/1‑150 (MQ=255)
                     aCTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTAt                                                                                                                              <  2:1395983/150‑1 (MQ=255)
                          cgcAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTc                                                                                                                        <  1:963853/151‑1 (MQ=255)
                             aaCTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGc                                                                                                                     >  1:855094/1‑151 (MQ=255)
                             aaCTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGc                                                                                                                     >  2:1096627/1‑151 (MQ=255)
                                tctcGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGCAGCTGGGGTGCTTTGTGCATCTGCAAGGTTTATCGGTTCCGcc                                                                                                                    >  2:992565/1‑149 (MQ=255)
                                 ctcGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCa                                                                                                                   <  1:340913/149‑1 (MQ=255)
                                   cGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCt                                                                                                               >  2:340214/1‑151 (MQ=255)
                                       gCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGcc                                                                                                            >  1:1465714/1‑150 (MQ=255)
                                       gCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGcc                                                                                                            >  1:67298/1‑150 (MQ=255)
                                       gCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGcc                                                                                                            >  1:522963/1‑150 (MQ=255)
                                          gTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAAcc                                                                                                        <  1:392098/151‑1 (MQ=255)
                                          gTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAAc                                                                                                         <  1:1314797/150‑1 (MQ=255)
                                            accaccAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCa                                                                                                       <  2:1163974/150‑1 (MQ=255)
                                            accaccAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAg                                                                                                      <  1:1288330/151‑1 (MQ=255)
                                             ccaccaGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAgg                                                                                                     <  2:877627/151‑1 (MQ=255)
                                                ccaGTACACCAAGCAACTGGCGACGTTGAAGAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCcg                                                                                                  <  1:811532/151‑1 (MQ=255)
                                                 caGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTTTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCcg                                                                                                  >  1:426690/1‑150 (MQ=255)
                                                  aGTACACCAAGCAACTGGCGACGTTGAATAATTGGTGCGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCcgac                                                                                                <  2:1368971/151‑1 (MQ=255)
                                                  aGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCcgac                                                                                                <  1:1045088/151‑1 (MQ=255)
                                                  aGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCcgac                                                                                                <  2:472826/151‑1 (MQ=255)
                                                     acacCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCcgacga                                                                                              <  2:1323238/151‑1 (MQ=255)
                                                     acacCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCCAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGAt                                                                                             >  2:303330/1‑151 (MQ=255)
                                                     acacCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGAt                                                                                             <  2:518036/151‑1 (MQ=255)
                                                     acacCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGAt                                                                                             <  2:1137080/151‑1 (MQ=255)
                                                        ccAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATccc                                                                                          >  1:1452881/1‑151 (MQ=255)
                                                             cAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTtcat                                                                                     <  1:1400382/151‑1 (MQ=255)
                                                                  ggCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCaaa                                                                                 <  2:228526/150‑1 (MQ=255)
                                                                   gCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAAcg                                                                               <  2:802198/151‑1 (MQ=255)
                                                                    cGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAAcgc                                                                              <  2:1308626/151‑1 (MQ=255)
                                                                     gACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAAcgcg                                                                             <  1:1451840/151‑1 (MQ=255)
                                                                     gACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAAcgcg                                                                             >  2:476093/1‑151 (MQ=255)
                                                                       cGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAAcgcg                                                                             >  2:318512/1‑149 (MQ=255)
                                                                          tGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGt                                                                        >  1:1410176/1‑151 (MQ=255)
                                                                             ataatTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTAc                                                                     >  2:92074/1‑151 (MQ=255)
                                                                               aatTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACa                                                                    <  2:1446694/150‑1 (MQ=255)
                                                                               aatTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAg                                                                   <  2:324811/151‑1 (MQ=255)
                                                                                atTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGt                                                                  >  1:87601/1‑151 (MQ=255)
                                                                                   ggTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGc                                                                <  1:1467725/150‑1 (MQ=255)
                                                                                   ggTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCg                                                               <  1:512815/151‑1 (MQ=255)
                                                                                    gTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTg                                                                 <  2:979935/148‑1 (MQ=255)
                                                                                        gCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACcgc                                                          >  1:600885/1‑151 (MQ=255)
                                                                                        gCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACcgc                                                          >  2:1479461/1‑151 (MQ=255)
                                                                                          cTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACcgcg                                                         >  1:459279/1‑150 (MQ=255)
                                                                                          cTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACcgcg                                                         >  2:901917/1‑150 (MQ=255)
                                                                                          cTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGt                                                        >  2:1208481/1‑151 (MQ=255)
                                                                                          cTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGt                                                        >  2:1343087/1‑151 (MQ=255)
                                                                                            aaaaaCGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGt                                                        <  1:1126012/149‑1 (MQ=255)
                                                                                            aaaaaCGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTgg                                                      >  2:280742/1‑151 (MQ=255)
                                                                                             aaaaCGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGt                                                     <  1:321291/151‑1 (MQ=255)
                                                                                              aaaCGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCg                                                               >  2:490834/1‑140 (MQ=255)
                                                                                              aaaCGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCg                                                               <  1:490841/140‑1 (MQ=255)
                                                                                              aaaCGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGtt                                                    <  2:410100/151‑1 (MQ=255)
                                                                                               aaCGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGttt                                                   >  2:1236255/1‑151 (MQ=255)
                                                                                                  gcgcgGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGtttttt                                                <  2:927126/151‑1 (MQ=255)
                                                                                                  gcgcgGAAACGTTCTTCTTTGACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGtttttt                                                <  1:26853/151‑1 (MQ=255)
                                                                                                   cgcgGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTc                                               <  1:229496/151‑1 (MQ=255)
                                                                                                   cgcgGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTc                                               >  2:161875/1‑151 (MQ=255)
                                                                                                      ggAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGc                                            >  2:1411386/1‑151 (MQ=255)
                                                                                                        aaaCGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGccc                                          >  1:352864/1‑151 (MQ=255)
                                                                                                        aaaCGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGccc                                          >  1:654500/1‑151 (MQ=255)
                                                                                                         caCGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGcccc                                         >  1:318155/2‑151 (MQ=255)
                                                                                                           cGTTCTTCTTTTACGGCGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCg                                       >  2:860583/1‑151 (MQ=255)
                                                                                                           cGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGccccc                                        <  2:884143/150‑1 (MQ=255)
                                                                                                           cGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCg                                       <  2:1015201/151‑1 (MQ=255)
                                                                                                              tcttctTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCggg                                     <  1:75250/150‑1 (MQ=255)
                                                                                                                ttcttTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCg                                  >  1:993408/1‑151 (MQ=255)
                                                                                                                  cttTTACGGAGGGGATGTATTTGAAACTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGcc                                >  1:192591/1‑151 (MQ=255)
                                                                                                                   ttttACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGcc                                >  2:205991/1‑150 (MQ=255)
                                                                                                                   ttttACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGcc                                >  1:795255/1‑150 (MQ=255)
                                                                                                                      tACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATc                             >  1:183681/1‑150 (MQ=255)
                                                                                                                      tACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATc                             >  1:1362843/1‑150 (MQ=255)
                                                                                                                       aCGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATc                             >  1:673647/1‑149 (MQ=255)
                                                                                                                       aCGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATc                             >  1:1441253/1‑149 (MQ=255)
                                                                                                                       aCGGAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAg                           <  1:1456827/151‑1 (MQ=255)
                                                                                                                         ggAGGGGATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGgt                         <  2:1433812/151‑1 (MQ=255)
                                                                                                                            ggggATGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCATGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGgt                         >  2:1352310/1‑148 (MQ=255)
                                                                                                                                aTGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTaaca                  <  1:28379/151‑1 (MQ=255)
                                                                                                                                aTGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTaac                   >  1:778325/1‑150 (MQ=255)
                                                                                                                                 tGTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTaaca                  >  1:1292401/1‑150 (MQ=255)
                                                                                                                                  gTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTaacaac                >  2:158549/1‑151 (MQ=255)
                                                                                                                                  gTATTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTaacaa                 <  2:1023813/150‑1 (MQ=255)
                                                                                                                                    aTTTGAAGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGt              <  1:243798/151‑1 (MQ=255)
                                                                                                                                         aaGCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATc         <  1:1173445/151‑1 (MQ=255)
                                                                                                                                           gCTGGGGTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCAt       >  2:827771/1‑151 (MQ=255)
                                                                                                                                               gggTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATc   <  1:429751/151‑1 (MQ=255)
                                                                                                                                                ggTGCTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCg  <  1:340751/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCG  >  NC_000913/2968011‑2968306

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: