Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,119,244 G→C 7.1% E324Q (GAG→CAG)  PP_4505 → sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,119,2440GC 7.1% 335.7 / 13.0 108E324Q (GAG→CAG) PP_4505sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (50/50);  new base C (5/3);  total (55/53)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

Only 100 of 109 total aligned reads displayed.

CCTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATCAGC  >  NC_002947/5119099‑5119393
                                                                                                                                                 |                                                                                                                                                     
ccTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGgc                                                                                                                                                   <  1:930614/150‑1 (MQ=255)
ccTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGgc                                                                                                                                                   >  1:154649/1‑150 (MQ=255)
    gAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAg                                                                                                                                               <  1:1796133/150‑1 (MQ=255)
    gAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAg                                                                                                                                               <  2:793571/150‑1 (MQ=255)
        aCCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAg                                                                                                                                           >  2:2088021/1‑150 (MQ=255)
          cTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGt                                                                                                                                         <  1:2088021/150‑1 (MQ=255)
              tCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGa                                                                                                                                     <  2:1488216/150‑1 (MQ=255)
              tCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGa                                                                                                                                     <  1:424123/150‑1 (MQ=255)
              tCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGa                                                                                                                                     <  1:2112048/150‑1 (MQ=255)
                aGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGc                                                                                                                                   >  1:466520/1‑150 (MQ=255)
                   gAAATCCTCGCCCTGGCCCGGGTGGATTCCGAACAAGCCCATGCCGAGCCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTg                                                                                                                                >  1:1167905/1‑150 (MQ=255)
                         cTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTGTGAGCTTGCAGGGt                                                                                                                          >  1:1695882/1‑150 (MQ=255)
                          tCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGttt                                                                                                                                       <  1:998470/136‑1 (MQ=255)
                          tCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGtt                                                                                                                         >  2:1806033/1‑150 (MQ=255)
                          tCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCCGGTCTCGCCGGATCAGGACGTACGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGtt                                                                                                                         >  1:1233350/1‑150 (MQ=255)
                            gCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCCACATATAGTCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTgg                                                                                                                       <  1:1352904/150‑1 (MQ=255)
                              ccTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGcc                                                                                                                     >  2:885995/1‑150 (MQ=255)
                              ccTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGcc                                                                                                                     >  1:796764/1‑150 (MQ=255)
                                tGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCaa                                                                                                                   >  1:261657/1‑150 (MQ=255)
                                  gCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTCCAGGCGCacgtc                                                                                                                                            >  2:1771047/1‑119 (MQ=255)
                                   cccGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGc                                                                                                                >  1:354535/1‑150 (MQ=255)
                                    ccGGGTGGAGGCCGCACAAGCCCATTCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGCTGCGCTGCTGAGCGCCCCGGATCAGTACGTGCGCCTGGACGCGCAGCCAGGTTTAAGCTTGCAGGGTTGGCCAACGCt                                                                                                               <  2:250823/150‑1 (MQ=255)
                                     cGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTg                                                                                                              <  2:1244863/150‑1 (MQ=255)
                                      gggTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGa                                                                                                             >  1:949910/1‑150 (MQ=255)
                                                cGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTGCTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATc                                                                                                           <  2:248719/142‑1 (MQ=255)
                                                cGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATc                                                                                                           >  1:248719/1‑142 (MQ=255)
                                                 gAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGcc                                                                                                  <  1:1838867/150‑1 (MQ=255)
                                                 gAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGACCAACGCTGATCGAGCGTGcc                                                                                                  >  2:1030651/1‑150 (MQ=255)
                                                     aaGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTg                                                                                              >  1:1857482/1‑150 (MQ=255)
                                                       gCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGGTGGCCAACTCTAATCGAGCGTGCCGTTGAc                                                                                            >  2:305068/1‑150 (MQ=255)
                                                         ccATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGa                                                                                                             >  1:1880251/1‑131 (MQ=255)
                                                         ccATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGa                                                                                                             <  2:1880251/131‑1 (MQ=255)
                                                         ccATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACaa                                                                                          >  1:912784/1‑150 (MQ=255)
                                                         ccATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACaa                                                                                          >  2:1531343/1‑150 (MQ=255)
                                                                gAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACctgctg                                                                                   >  2:1381154/1‑150 (MQ=255)
                                                                gAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACctgctg                                                                                   >  2:2168142/1‑150 (MQ=255)
                                                                  aCCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTgcg                                                                                 >  2:1452840/1‑150 (MQ=255)
                                                                    acctCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCa                                                                               >  2:1329348/4‑150 (MQ=255)
                                                                       tCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATg                                                                            <  1:2168142/150‑1 (MQ=255)
                                                                         gACCTCAATGCCCTACTCGGCAGCGTGCTCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGcc                                                                          <  2:1154658/150‑1 (MQ=255)
                                                                         gACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGcc                                                                          <  1:1145129/150‑1 (MQ=255)
                                                                         gACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTCCAGGCGCac                                                                                                                                               >  1:570639/1‑80 (MQ=255)
                                                                         gACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTCCAGGCGCac                                                                                                                                               <  2:570639/81‑2 (MQ=255)
                                                                           ccTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCAGAAAAGGCGGTGCGCCTGGAGGCGCAGCCAGGGTTGAGCTTGCAGCGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCt                                                                        >  2:1533359/1‑150 (MQ=255)
                                                                            cTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTg                                                                       >  2:1769886/1‑150 (MQ=255)
                                                                            cTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTg                                                                       >  1:1640188/1‑150 (MQ=255)
                                                                             tCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTgc                                                                      >  1:1188041/1‑150 (MQ=255)
                                                                             tCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTgc                                                                      >  1:1613325/1‑150 (MQ=255)
                                                                             tCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTgc                                                                      >  1:1186265/1‑150 (MQ=255)
                                                                               aaTGCCCTACTCGGCAGCGTGCGCAAGCAGGCGCTGATTGGTGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCCAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTgcgc                                                                    <  2:1398220/150‑1 (MQ=255)
                                                                                aTGCCCTCCTCGGCAGCGTGCGCAAGGACGCCCTGATGAGCGCGCCGGATCAGGACGTGCGCTTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCt                                                                   <  2:469773/150‑1 (MQ=255)
                                                                                   cccTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACg                                                                                                                 <  1:2108448/101‑1 (MQ=255)
                                                                                    ccTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCaa                                                               >  1:956212/1‑150 (MQ=255)
                                                                                      tACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGtt                                                                                                                         <  1:1117390/90‑1 (MQ=255)
                                                                                      tACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGtt                                                                                                                         >  2:1117390/1‑90 (MQ=255)
                                                                                         tCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCcgg                                                          >  2:1021143/1‑150 (MQ=255)
                                                                                           ggCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCcggcg                                                        >  2:823245/1‑150 (MQ=255)
                                                                                              aGCGTGCTCAAGTACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGTCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGc                                                     <  2:466520/150‑1 (MQ=255)
                                                                                                    cgcAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCg                                               <  2:1640188/150‑1 (MQ=255)
                                                                                                     gcAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAgg                                                                                                                            <  1:929457/72‑1 (MQ=255)
                                                                                                     gcAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAgg                                                                                                                            >  2:929457/1‑72 (MQ=255)
                                                                                                     gcAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTACAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCgg                                              <  2:1283961/150‑1 (MQ=255)
                                                                                                       aaGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCTGGCCCGCCgggt                                            <  1:1769886/150‑3 (MQ=255)
                                                                                                         ggACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGc                                                                                                                      <  1:1718326/74‑1 (MQ=255)
                                                                                                         ggACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGc                                                                                                                      >  2:1718326/1‑74 (MQ=255)
                                                                                                         ggACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGa                                          <  1:823245/150‑1 (MQ=255)
                                                                                                         ggACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGa                                          <  2:2156933/150‑1 (MQ=255)
                                                                                                         ggACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGa                                          <  2:854450/150‑1 (MQ=255)
                                                                                                          gACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAg                                         <  2:1628368/150‑1 (MQ=255)
                                                                                                          gACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAg                                         <  1:1120210/150‑1 (MQ=255)
                                                                                                          gACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGATTCAACCCGGCGGGCCAGCCGGTCGAg                                         <  1:748772/150‑1 (MQ=255)
                                                                                                           aCGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCTGGCCCGCCgggttgaa                                        <  2:1095566/150‑7 (MQ=255)
                                                                                                           aCGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCTGGCCCGCCgggttgaa                                        <  2:1695882/150‑7 (MQ=255)
                                                                                                            cGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGt                                       <  2:956212/150‑1 (MQ=255)
                                                                                                            cGCCCTGCTGAGCGCCCCGACCCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTCAGCTTGCAGGGCTGCCCAACGCTCATCGAGCGTGATGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGt                                       <  1:1329348/150‑1 (MQ=255)
                                                                                                              cccTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCa                                     <  2:796764/150‑1 (MQ=255)
                                                                                                                ctgctgAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCAAAAACTGTGTGTTTGTGCCGTTGACAAGCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAgc                                   <  2:694787/150‑1 (MQ=255)
                                                                                                                 tgctgAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAAc                                                                                         <  1:543585/95‑1 (MQ=255)
                                                                                                                 tgctgAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAAc                                                                                         >  2:543585/1‑95 (MQ=255)
                                                                                                                    tgAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGcc                                                    <  1:1081156/129‑1 (MQ=255)
                                                                                                                    tgAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGcc                                                    >  2:1081156/1‑129 (MQ=255)
                                                                                                                    tgAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTCGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGCGCCAGCCGGTCGAGGTCAGcgccg                               >  2:915744/1‑150 (MQ=255)
                                                                                                                            ccGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAg                       >  1:1207526/1‑150 (MQ=255)
                                                                                                                               gaccAGGACGTGCGCCTGGCGGCGCAGGCATGTTTAAGCTTGCAGGGTTGGCCACCGCTGGTCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAg                    <  2:881732/147‑1 (MQ=255)
                                                                                                                                 tCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGa                  >  1:1929291/1‑150 (MQ=255)
                                                                                                                                 tCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGa                  >  2:1426840/1‑150 (MQ=255)
                                                                                                                                  cAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCAAGGTCAGCGCCGCGCGCGAGCAGGAc                 <  1:1531343/150‑1 (MQ=255)
                                                                                                                                  acctggctGCGCCTCCAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGAc                 >  2:639308/8‑150 (MQ=255)
                                                                                                                                    ctggctGCGCCTCCAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCg               >  2:2198816/6‑150 (MQ=255)
                                                                                                                                      aCGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCa             <  1:1426840/150‑1 (MQ=255)
                                                                                                                                        gTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCTGGCCCGCCgggt                                            <  1:1911263/117‑3 (MQ=255)
                                                                                                                                        gTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCTGGCCCGCCgggt                                            >  2:1911263/1‑115 (MQ=255)
                                                                                                                                         tGCGCCTCCAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAAcc                                                             <  1:1625635/99‑1 (MQ=255)
                                                                                                                                            gcCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCAGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGa       <  1:163679/150‑1 (MQ=255)
                                                                                                                                             ccTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGAt      <  2:2236202/150‑1 (MQ=255)
                                                                                                                                              cTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATc           <  1:1784201/144‑1 (MQ=255)
                                                                                                                                              cTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATc           >  2:1784201/1‑144 (MQ=255)
                                                                                                                                              cTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATc     >  1:516509/1‑150 (MQ=255)
                                                                                                                                              cTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATc     >  1:1094069/1‑150 (MQ=255)
                                                                                                                                                 gAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATCAgc  >  1:1244730/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CCTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATCAGC  >  NC_002947/5119099‑5119393

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