Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,165,671 C→T 5.5% G159R (GGG→AGG)  PP_3666 ← MHS family metabolite MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,165,6710CT 5.5% 221.8 / 2.4 73G159R (GGG→AGG) PP_3666MHS family metabolite MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (34/35);  new base T (2/2);  total (36/37)
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.03e-01

CGCCAGCAGGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGATGCAGAAGGTG  >  NC_002947/4165525‑4165819
                                                                                                                                                  |                                                                                                                                                    
cgccAGCAGGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAgg                                                                                                                                                   <  2:1821701/150‑1 (MQ=255)
    agcagGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAggcggc                                                                                                                                               <  1:1627285/150‑1 (MQ=255)
       agGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGTCCAGGCGGcgcg                                                                                                                                            <  2:702081/150‑1 (MQ=255)
        ggAACGGCACGCGCCCGACACAATAATGATGTGTTCCCCTGCTCAGCTGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTGCAAGCCACGCCGCCTGGgcggcgtgg                                                                                                                                           <  2:958680/150‑4 (MQ=255)
           aCGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACt                                                                                                                                        >  2:515798/1‑150 (MQ=255)
               cAGGCGCCAGCCCCACTCATGCAGGGCTTCATCTTTGCGCTTTTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGc                                                                                                                                    <  2:628227/150‑1 (MQ=255)
                      cAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGaa                                                                                                                              >  2:1439988/1‑149 (MQ=255)
                      cAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGaa                                                                                                                              <  1:1439988/149‑1 (MQ=255)
                           ccACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGc                                                                                                                        >  1:1419621/1‑150 (MQ=255)
                           ccACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGc                                                                                                                        >  1:1420226/1‑150 (MQ=255)
                             aCTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCgg                                                                                                                      >  2:213185/1‑150 (MQ=255)
                                 aTGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGTGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGcccc                                                                                                                  >  2:1905949/1‑150 (MQ=255)
                                       ggCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCg                                                                                                                       >  1:1681960/1‑139 (MQ=255)
                                       ggCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCg                                                                                                                       <  2:1681960/139‑1 (MQ=255)
                                       ggCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTa                                                                                                            <  1:1362160/150‑1 (MQ=255)
                                         cTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACt                                                                                                          >  1:847458/1‑150 (MQ=255)
                                          ttCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTc                                                                                                         >  2:1148415/1‑150 (MQ=255)
                                             cccgaagAGCCGCTCATGCAGCACCTCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCAcc                                                                                                      <  1:969645/144‑1 (MQ=255)
                                             aCTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCAcc                                                                                                      >  1:1386925/1‑150 (MQ=255)
                                                gctgAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGa                                                                                                   <  2:761181/150‑1 (MQ=255)
                                                  tgAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAATGCAGCGGCCCCGGCGTACTCACCCGAgg                                                                                                 <  1:1708991/150‑1 (MQ=255)
                                                  tgAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCCGCACGATACTGCTTTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAgg                                                                                                 >  2:1315858/1‑150 (MQ=255)
                                                        gcTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGaaa                                                                                           >  1:977380/1‑150 (MQ=255)
                                                         cTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCTGGgcggcgtggc                                                                                                                                          <  2:600261/102‑5 (MQ=255)
                                                         cTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCTGGgcggcgtggc                                                                                                                                          >  1:600261/1‑98 (MQ=255)
                                                         cTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCa                                                                                                                          <  2:1438489/118‑1 (MQ=255)
                                                         cTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCa                                                                                                                          >  1:1438489/1‑118 (MQ=255)
                                                          tCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAAc                                                                                         >  2:1760835/1‑150 (MQ=255)
                                                           cATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAAcc                                                                                        <  1:917843/150‑1 (MQ=255)
                                                           cATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAAcc                                                                                        <  1:917819/150‑1 (MQ=255)
                                                                    ccgtgacgaACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCa                                                                               <  1:1905949/146‑1 (MQ=255)
                                                                    ccGCGACGAACGCCGCACGAAACAACGACCACCCCGCCGTGCTGTGCGGGACGCTCCCGGTGTCTACGCCACGCCGCCCCGGCCGCGAGTACCCGGCTAGGAACGCAAAGGGGACGGAGTACGCACCTGAGGCGGAAAACACCTGTACCa                                                                               <  2:515839/150‑1 (MQ=255)
                                                                        gacgaACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCg                                                                           <  1:1166467/150‑1 (MQ=255)
                                                                        gacgaACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCg                                                                           >  2:1166467/1‑150 (MQ=255)
                                                                           gaACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGAt                                                                        <  1:1536938/150‑1 (MQ=255)
                                                                           gaACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGAt                                                                        <  1:427828/150‑1 (MQ=255)
                                                                                gtccACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATcagca                                                                   >  1:65785/4‑150 (MQ=255)
                                                                                 cgcACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAgg                                                                                                 <  1:467691/119‑1 (MQ=255)
                                                                                 cgcACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAgg                                                                                                 >  2:467691/1‑119 (MQ=255)
                                                                                   cACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGc                                                                                                <  1:684987/118‑1 (MQ=255)
                                                                                   cACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGc                                                                                                >  2:684987/1‑118 (MQ=255)
                                                                                      cGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATcagcaacagca                                                             >  2:1134629/1‑150 (MQ=255)
                                                                                       gAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATcagcaacagcaa                                                            >  1:810736/1‑150 (MQ=255)
                                                                                        aacaacaacCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAg                                                           >  1:68482/1‑150 (MQ=255)
                                                                                           aacaacCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGcc                                                        >  2:1899627/1‑150 (MQ=255)
                                                                                             caacCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCgg                                                      >  2:396969/1‑150 (MQ=255)
                                                                                              aacCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGt                                                     >  1:1906811/1‑150 (MQ=255)
                                                                                              aacCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGt                                                     <  2:68482/150‑1 (MQ=255)
                                                                                                 cccGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGcc                                                  >  1:263243/1‑150 (MQ=255)
                                                                                                 cccGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGcc                                                  >  2:467437/1‑150 (MQ=255)
                                                                                                        ccgTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGcc                                           <  1:1101513/150‑1 (MQ=255)
                                                                                                        ccgTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGcc                                           <  1:467437/150‑1 (MQ=255)
                                                                                                                gCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTgc                                   <  1:213185/150‑1 (MQ=255)
                                                                                                                    gCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCAt                               >  1:1605473/1‑150 (MQ=255)
                                                                                                                     cACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACGCCAGTGGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATa                              >  1:1277430/1‑150 (MQ=255)
                                                                                                                     cACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATa                              <  2:263243/150‑1 (MQ=255)
                                                                                                                      aCGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAg                             <  1:1116339/150‑1 (MQ=255)
                                                                                                                          tGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCg                         <  1:1558111/150‑1 (MQ=255)
                                                                                                                          tGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCg                         <  1:902923/150‑1 (MQ=255)
                                                                                                                           gCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCgg                        <  2:1386925/150‑1 (MQ=255)
                                                                                                                                   aCAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCTGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGc                <  1:515798/150‑1 (MQ=255)
                                                                                                                                    cAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGcc               >  1:1073872/1‑150 (MQ=255)
                                                                                                                                     aaGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCg              >  1:11306/1‑150 (MQ=255)
                                                                                                                                      aGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAACCCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGGCGa             <  2:177039/150‑1 (MQ=255)
                                                                                                                                      aGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGa             <  1:1597839/150‑1 (MQ=255)
                                                                                                                                      aGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGa             >  2:752757/1‑150 (MQ=255)
                                                                                                                                      aGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGa             <  2:11306/150‑1 (MQ=255)
                                                                                                                                          gcgccgcCTGGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAGCAAGCCGATGCa         >  1:1608294/2‑150 (MQ=255)
                                                                                                                                          aCGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGc                                            <  1:1600326/115‑1 (MQ=255)
                                                                                                                                          aCGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGc                                            >  2:1600326/1‑115 (MQ=255)
                                                                                                                                             ccggcCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGATGCAGaa      <  2:916876/150‑1 (MQ=255)
                                                                                                                                             ccgccCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGATGCAGaa      >  1:1471585/1‑150 (MQ=255)
                                                                                                                                              cgccCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGATGCAGAAg     >  1:1028606/1‑150 (MQ=255)
                                                                                                                                                cccAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGATGCAGAAGGt   >  1:1478184/1‑150 (MQ=255)
                                                                                                                                                cccAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGATGCAGAAGGt   <  1:1899627/150‑1 (MQ=255)
                                                                                                                                                cccAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCACTTTCTGCGCATAGCCCGGCAACCAGCCGATGCAGAAGGt   <  1:1315858/150‑1 (MQ=255)
                                                                                                                                                 ccAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCGCAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGATGCAGAAGGTg  <  1:1280782/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGCCAGCAGGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGATGCAGAAGGTG  >  NC_002947/4165525‑4165819

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