Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,709,547 G→A 6.5% intergenic (+67/+45) PP_5006 → / ← PP_5007 TetR family transcriptional regulator/polyhydroxyalkanoate granule‑associated protein GA2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,709,5470GA 6.5% 889.3 / 32.3 290intergenic (+67/+45)PP_5006/PP_5007TetR family transcriptional regulator/polyhydroxyalkanoate granule‑associated protein GA2
Reads supporting (aligned to +/- strand):  ref base G (130/140);  new base A (11/8);  total (141/149)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 297 total aligned reads displayed.

GATGATGCTGGTGGCGCCGCATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCAGGCGCAGCGCTGGCAGCTGG  >  NC_002947/5709399‑5709696
                                                                                                                                                    |                                                                                                                                                     
gatgatGCTGGTGGCGCCGCATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGgcgc                                                                                                                                                      >  1:2234095/1‑150 (MQ=255)
  tgatgCTGGTGGCGCCGCATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGTGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCAt                                                                                                                                                    <  1:626028/150‑1 (MQ=255)
   gatgCTGGTGGCGCCGCATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATc                                                                                                                                                   >  2:3775260/1‑150 (MQ=255)
   gatgCTGGTGGCGCCGCATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATc                                                                                                                                                   >  1:2205781/1‑150 (MQ=255)
   gatgCTGGTGGCGCCGCATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAAACCGCTCCTACACGGTTTGGCGCaaa                                                                                                                                                   >  1:5770975/1‑148 (MQ=255)
        tggtggCGCCGCATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTa                                                                                                                                              >  2:2865599/1‑150 (MQ=255)
        tggtggCGCCGCATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTa                                                                                                                                              >  2:909927/1‑150 (MQ=255)
        gggtggCGCCGCATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTa                                                                                                                                              <  2:6436664/149‑1 (MQ=255)
         ggtggCGCCGCATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGGAg                                                                                                                                             <  1:6316841/150‑1 (MQ=255)
          gtggCGCCGCATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAgg                                                                                                                                            >  2:1315772/1‑150 (MQ=255)
            ggCGCCGCATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAg                                                                                                                                          >  2:4264687/1‑150 (MQ=255)
             gcgcCGCATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGc                                                                                                                                         <  1:934877/150‑1 (MQ=255)
             gcgcCGCATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGc                                                                                                                                         <  1:2694549/150‑1 (MQ=255)
             gcgcCGCATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGATGCGCCAAACCTTGTAGGAGc                                                                                                                                         <  1:6346270/150‑1 (MQ=38)
             gcgcCGCATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGAGGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGGTTGGCGCCTCTTGTAGGAGc                                                                                                                                         >  2:4034729/1‑150 (MQ=255)
               gccgcATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTCCCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCgg                                                                                                                                       <  1:4183093/150‑1 (MQ=255)
               gccgcATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACCAGGTTTGGCGCATCTTGTAGGAGCgg                                                                                                                                       >  2:6248461/1‑150 (MQ=255)
               gccgcATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCgg                                                                                                                                       <  1:4343090/150‑1 (MQ=255)
               gccgcATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGATGCGCCAAACCTTGTAGGAGCgg                                                                                                                                       <  1:4636296/150‑1 (MQ=37)
                ccgcATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGGGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCggg                                                                                                                                      >  1:1721073/1‑150 (MQ=255)
                ccgcATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCggg                                                                                                                                      <  1:6173728/150‑1 (MQ=255)
                ccgcATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCggg                                                                                                                                      >  1:4086857/1‑150 (MQ=255)
                ccgcATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCggg                                                                                                                                      <  1:4825268/150‑1 (MQ=255)
                  gcctCTGCAGGCCCCGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGGACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGtt                                                                                                                                    <  2:6740592/147‑1 (MQ=255)
                       tGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACCAGGTTTGGCGCATCTTGTCGGAGCGGGTTTAccc                                                                                                                               >  2:5795521/1‑150 (MQ=255)
                       tGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTAccc                                                                                                                               <  1:6620338/150‑1 (MQ=255)
                       tGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCACCAACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTAccc                                                                                                                               >  1:3644617/1‑150 (MQ=255)
                        gCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCcg                                                                                                                              >  1:5710246/1‑150 (MQ=255)
                         caggcccaggccCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCcgc                                                                                                                             >  1:6215334/1‑150 (MQ=255)
                         caggcccaggccCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCcgc                                                                                                                             >  1:5041152/1‑150 (MQ=255)
                          aggcccaggcccGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCcgcg                                                                                                                            >  2:3084253/1‑150 (MQ=255)
                          aggcccaggcccGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCcgcg                                                                                                                            >  2:6820622/1‑150 (MQ=255)
                           ggccgcggcccGGGCGGGGGGGGAGCTGCGGGCTGTGAAGTACCTGGTGGGGTTGGCCTGCTGGAGATGTCGGGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCGTACAAGGTTTGGCGCATCTTGTAGGTGGGGGTTTACCCGCga                                                                                                                           <  1:6634351/144‑1 (MQ=255)
                             cccaggcccGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCgaag                                                                                                                         >  2:3288961/1‑150 (MQ=255)
                             cccaggcccGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCgaaga                                                                                                                        >  2:1607754/1‑150 (MQ=255)
                              ccaggcccGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCgaaga                                                                                                                        >  1:3977440/1‑150 (MQ=255)
                               caggcccGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACACGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCgaagaa                                                                                                                       <  1:5644067/150‑1 (MQ=255)
                               caggcccGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCgaagaa                                                                                                                       <  2:4365448/150‑1 (MQ=255)
                                    ccGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGa                                                                                                                  <  2:6174428/150‑1 (MQ=255)
                                            cggcggAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTat                                                                                                          >  1:6729218/1‑150 (MQ=255)
                                                  agcagcTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACATGGTTTGGCGCATCTTGTAGGAGGGGGTTTACCCGCGAAGAAGGCGACGCGGTatcagatc                                                                                                    <  1:5990948/150‑1 (MQ=255)
                                                   gcagcTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTatcagatca                                                                                                   <  1:4549511/150‑1 (MQ=255)
                                                   gcagcTGGCAGTGAAGTACAGGGAGGGTTAGGCCTGATGTAGATGTAGTGTTGCAGCCTACGCATTCGCGGGAAAACCCTCACCTACAAGTTTTGCCTCAACTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTatcagatca                                                                                                   <  1:158431/150‑1 (MQ=37)
                                                           cAGTTACGTACCTGGAGGGGTAGGCCTGCTGGAGTTGTAGTTTTGCAGCCGACGAATTCGCGGTTAAACCCGCTCCTACAAGGTTTGGCGCATATTGTCGGTGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCg                                                                                           <  1:5892092/150‑1 (MQ=255)
                                                           cAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACACCCAACAAAAACGTTTGGCGCCAAATTTATGAGCCGGTTTAACCCCGAAGAAGGCGACCCCGGAAAAAATCAGGGTaacc                                                                                           >  1:3000657/1‑147 (MQ=255)
                                                           cAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAAACCGCTACTAAAAGGTTTGGCGGATCTTGTAGGATCGGGTTTAACCGCGAAGAAGGCGAAGCGGTATCAGATAATGGTACCg                                                                                           >  2:2174048/1‑150 (MQ=255)
                                                            aGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCgg                                                                                          >  2:4923875/1‑150 (MQ=255)
                                                            aGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGt                                                                                         >  2:2463499/1‑150 (MQ=255)
                                                             gTGAAGTCCCTGGAGGGGTAGGCCTGCTGGAGATGCAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGt                                                                                         <  1:3839014/150‑1 (MQ=255)
                                                             gTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGt                                                                                         <  2:3693812/150‑1 (MQ=255)
                                                                   tACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTatca                                                                                                        >  2:6035363/1‑129 (MQ=255)
                                                                      cTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCgg                                                                                >  2:4141215/1‑150 (MQ=255)
                                                                              gTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTAACGGTGCCAGTCGGTGCCGACg                                                                        >  2:4894634/1‑150 (MQ=255)
                                                                              gTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCAACAAGGTTTGGCGCATCGTGTAGGAGCGGGTTTACCCGCGACGAAGGCGACGCGGGATCAGATCAGGGTACCGGTGCCCGTCGGTGCCCACg                                                                        >  2:4695214/1‑150 (MQ=255)
                                                                               tAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCAACTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTAACAGATCAGGGTACCGGTGCCAGTCGGTGCCGACgg                                                                       >  1:2761092/1‑150 (MQ=255)
                                                                                 ggCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCg                                                                     <  2:3158517/150‑1 (MQ=255)
                                                                                      gctgGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGtgctg                                                                <  2:2330287/150‑1 (MQ=255)
                                                                                        tgGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCCCACCAACAAGGTTTGTCGCAACATGTTGGAGCGGGTTTAACCCCCAAAAAAGCGACCCGGTATCAGATAAAGGTACACGAGCCACTCCGTGCCCACGGCGTtcagct                                                              >  1:688860/1‑150 (MQ=255)
                                                                                          gagaTGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCCCTCCTACAAGGTTTTGCGCATCTTGTGGGAGCGGGTTTAACCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGGTGCTgg                                                            >  1:1396521/1‑150 (MQ=255)
                                                                                               gTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCggg                                                       <  1:3459839/150‑1 (MQ=255)
                                                                                               gTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCACATCGTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCggg                                                       >  2:1195618/1‑150 (MQ=255)
                                                                                                tAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCgggg                                                      <  1:6826857/150‑1 (MQ=255)
                                                                                                tAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCGTGTAGGTGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCgggg                                                      >  1:374319/1‑150 (MQ=255)
                                                                                                   tgtTGCAGCCGACGCATTCGCGGGTAAAACCCCGCCAACAAGGTTTGGCGCATCTTGTAGGAGCGGGGTTACCCGCGAAGAAGGCGCCGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGGGCTGCTGGCCGGGGTgg                                                   >  1:2411549/1‑150 (MQ=255)
                                                                                                     ttGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACACGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTggcg                                                 <  1:1307726/150‑1 (MQ=255)
                                                                                                       gCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTggcggc                                               >  2:5139310/1‑150 (MQ=255)
                                                                                                       gCAGCCGACGCATTCGCGGGTAAACCCGCTCCGCCCAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTggcggc                                               <  2:3583112/150‑1 (MQ=255)
                                                                                                        cAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGTTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCa                                              <  1:4941203/150‑1 (MQ=255)
                                                                                                        cAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCa                                              <  1:3084253/150‑1 (MQ=255)
                                                                                                        cAGCCGACGCATTCGCGGGTAAACCCGCTACTACAAGGTTTGTCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTAACAGATCAGGGTACCGGTTCCACTCGGTGCCGACCGCGTGCTGCTGGCCGGGGTGGCGGCa                                              >  1:1064128/1‑150 (MQ=39)
                                                                                                         aGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAg                                             <  1:2529398/150‑1 (MQ=255)
                                                                                                           ccGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTGGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCCGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGtt                                           >  2:6592217/1‑150 (MQ=255)
                                                                                                                   ttCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGc                                   >  1:5086451/1‑150 (MQ=255)
                                                                                                                       cgGGTAACCCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCtgtg                               <  1:6108947/150‑1 (MQ=255)
                                                                                                                        gggTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTgg                              >  1:3274994/1‑150 (MQ=255)
                                                                                                                              aCCCGCTCCTACAAGATGCGCCAAACCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGc                                                                                       >  2:3442143/1‑87 (MQ=17)
                                                                                                                                ccGCTCCTACTAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCa                      <  2:5330578/150‑1 (MQ=255)
                                                                                                                                ccGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCa                      >  1:6513667/1‑150 (MQ=255)
                                                                                                                                ccGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCa                      >  1:5404719/1‑150 (MQ=255)
                                                                                                                                ccGCTCCTACAAGATGCGCCAAACCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTCCCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCa                      >  2:101673/1‑150 (MQ=37)
                                                                                                                                ccGCTCCTACAAGATGCGCCAAACCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGctgctg                                                             <  2:4991022/111‑1 (MQ=21)
                                                                                                                                ccGCTCCTACAAGATGCGCCAAACCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGct                      >  1:4212728/1‑149 (MQ=37)
                                                                                                                                ccGCTCCTACAAGATGCGCCAAACCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCa                      >  2:1721258/1‑150 (MQ=37)
                                                                                                                                 cGCTCCTACACGGTTTGGCGCATCTTGTAGGAGCGGGATTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCAg                     <  1:868320/150‑1 (MQ=255)
                                                                                                                                 cGCTCCTACAAGATGCGCCAAACCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGc                                                                                       >  2:5053099/1‑84 (MQ=17)
                                                                                                                                 cGCTCCTACAAGATGCGCCAAACCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCAg                     >  1:1160812/1‑150 (MQ=38)
                                                                                                                                  gCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCAgg                    >  1:5563869/1‑150 (MQ=255)
                                                                                                                                   cTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTAGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCAGgc                   <  2:256499/150‑1 (MQ=255)
                                                                                                                                    tCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCAGgcg                  >  2:214024/1‑150 (MQ=255)
                                                                                                                                       tACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCAGGCGCAg               >  2:5565875/1‑150 (MQ=255)
                                                                                                                                       tACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCAGGCGCAg               >  1:6744763/1‑150 (MQ=255)
                                                                                                                                       tACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCCTGCTGCTGGCCGGGGTGGCGGCAGt                                            <  2:40128/121‑1 (MQ=255)
                                                                                                                                       tACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCCTGCTGCTGGCCGGGGTGGCGGCAGt                                            >  1:40128/1‑121 (MQ=255)
                                                                                                                                        aCAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCAGGCGCAgc              >  2:4580293/1‑150 (MQ=255)
                                                                                                                                         cAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCAGGCGCAgcg             >  2:1959282/1‑150 (MQ=255)
                                                                                                                                         cAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCAGGCGCAgcg             >  2:1539907/1‑150 (MQ=255)
                                                                                                                                         cAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCAGGCGCAgcg             >  1:2294357/1‑150 (MQ=255)
                                                                                                                                            ggTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCAGGCGCAGCGCTg          >  2:1414728/1‑150 (MQ=255)
                                                                                                                                             gTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGGTGCCGGTGCTGTGGTCGCTGCAGGCGCAGCGCTgg         <  2:5509838/150‑1 (MQ=255)
                                                                                                                                                    gcATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCAGGCGCAGCGCTGGCAGCTgg  >  1:6071505/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATGATGCTGGTGGCGCCGCATCTGCAGGCCCAGGCCCGGGCGGCGGCGGAGCAGCTGGCAGTGAAGTACCTGGAGGGGTAGGCCTGCTGGAGATGTAGTGTTGCAGCCGACGCATTCGCGGGTAAACCCGCTCCTACAAGGTTTGGCGCATCTTGTAGGAGCGGGTTTACCCGCGAAGAAGGCGACGCGGTATCAGATCAGGGTACCGGTGCCAGTCGGTGCCGACGGCGTGCTGCTGGCCGGGGTGGCGGCAGTTGCCGGTGCTGTGGTCGCTGCAGGCGCAGCGCTGGCAGCTGG  >  NC_002947/5709399‑5709696

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