Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,119,243 G→C 7.4% L323L (CTG→CTC)  PP_4505 → sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,119,2430GC 7.4% 333.7 / 13.0 108L323L (CTG→CTC) 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 = 9.48e-01

Only 100 of 109 total aligned reads displayed.

GACCGCCTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATC  >  NC_002947/5119094‑5119390
                                                                                                                                                     |                                                                                                                                                   
gACCGCCTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTg                                                                                                                                                     <  2:927107/150‑1 (MQ=255)
     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                                                                                                                                  <  1:2112048/150‑1 (MQ=255)
                   tCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGa                                                                                                                                  <  2:1488216/150‑1 (MQ=255)
                   tCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGa                                                                                                                                  <  1:424123/150‑1 (MQ=255)
                     aGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGc                                                                                                                                >  1:466520/1‑150 (MQ=255)
                        gAAATCCTCGCCCTGGCCCGGGTGGATTCCGAACAAGCCCATGCCGAGCCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTg                                                                                                                             >  1:1167905/1‑150 (MQ=255)
                           aTCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAg                                                                                                                          <  2:1141570/150‑1 (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)
                                   ccTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGcc                                                                                                                  >  2:885995/1‑150 (MQ=255)
                                   ccTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGcc                                                                                                                  >  1:796764/1‑150 (MQ=255)
                                     tGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCaa                                                                                                                >  1:261657/1‑150 (MQ=255)
                                       gCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCTCCGGATCAGGACGTGCGCCTCCAGGCGCacgtc                                                                                                                                         <  1:1771047/123‑5 (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)
                                                   gCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAACGACGCCCTGCTGAGCGCGCCGGATCAGGACGTTCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGt                                                                                                  >  1:133592/1‑150 (MQ=255)
                                                     cGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTGCTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATc                                                                                                        <  2:248719/142‑1 (MQ=255)
                                                     cGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGc                                                                                                >  1:245935/1‑150 (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                                                                                       >  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)
                                                                           gTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAAt                                                                          >  1:881732/1‑150 (MQ=255)
                                                                            tCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATg                                                                         <  1:2168142/150‑1 (MQ=255)
                                                                              gACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGcc                                                                       <  1:1145129/150‑1 (MQ=255)
                                                                              gACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTCCAGGCGCac                                                                                                                                            <  2:570639/81‑2 (MQ=255)
                                                                              gACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTCCAGGCGCac                                                                                                                                            >  1:570639/1‑80 (MQ=255)
                                                                                ccTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCAGAAAAGGCGGTGCGCCTGGAGGCGCAGCCAGGGTTGAGCTTGCAGCGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCt                                                                     >  2:1533359/1‑150 (MQ=255)
                                                                                 cTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGCTCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTg                                                                    <  2:585059/150‑1 (MQ=255)
                                                                                 cTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTg                                                                    >  1:1640188/1‑150 (MQ=255)
                                                                                 cTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTg                                                                    >  2:1769886/1‑150 (MQ=255)
                                                                                  tCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTgc                                                                   >  1:1186265/1‑150 (MQ=255)
                                                                                  tCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTgc                                                                   >  1:1188041/1‑150 (MQ=255)
                                                                                  tCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTgc                                                                   >  1:1613325/1‑150 (MQ=255)
                                                                                    aaTGCCCTACTCGGCAGCGTGCGCAAGCAGGCGCTGATTGGTGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCCAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTgcgc                                                                 <  2:1398220/150‑1 (MQ=255)
                                                                                     aTGCCCTCCTCGGCAGCGTGCGCAAGGACGCCCTGATGAGCGCGCCGGATCAGGACGTGCGCTTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCt                                                                <  2:469773/150‑1 (MQ=255)
                                                                                        cccTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACg                                                                                                              >  2:2108448/1‑101 (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                                                                                                                   >  2:1718326/1‑74 (MQ=255)
                                                                                                              ggACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGc                                                                                                                   <  1:1718326/74‑1 (MQ=255)
                                                                                                              ggACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGa                                       <  1:823245/150‑1 (MQ=255)
                                                                                                              ggACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGa                                       <  2:854450/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                                                                                      >  2:543585/1‑95 (MQ=255)
                                                                                                                      tgctgAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAAc                                                                                      <  1:543585/95‑1 (MQ=255)
                                                                                                                         tgAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGcc                                                 >  2:1081156/1‑129 (MQ=255)
                                                                                                                         tgAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGcc                                                 <  1:1081156/129‑1 (MQ=255)
                                                                                                                         tgAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTCGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGCGCCAGCCGGTCGAGGTCAGcgccg                            >  2:915744/1‑150 (MQ=255)
                                                                                                                                 ccGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAg                    >  1:1207526/1‑150 (MQ=255)
                                                                                                                                    gaccAGGACGTGCGCCTGGCGGCGCAGGCATGTTTAAGCTTGCAGGGTTGGCCACCGCTGGTCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAg                 <  2:881732/147‑1 (MQ=255)
                                                                                                                                      tCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGa               >  2:1426840/1‑150 (MQ=255)
                                                                                                                                      tCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGa               >  1:1929291/1‑150 (MQ=255)
                                                                                                                                       cAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCAAGGTCAGCGCCGCGCGCGAGCAGGAc              <  1:1531343/150‑1 (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)
                                                                                                                                              tGCGCCTCCAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAAcc                                                          >  2:1625635/1‑99 (MQ=255)
                                                                                                                                                 gcCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCAGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGa    <  1:163679/150‑1 (MQ=255)
                                                                                                                                                  ccTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGAt   <  2:2236202/150‑1 (MQ=255)
                                                                                                                                                   cTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATc        >  2:1784201/1‑144 (MQ=255)
                                                                                                                                                   cTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATc  >  1:516509/1‑150 (MQ=255)
                                                                                                                                                   cTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATc  >  1:1094069/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GACCGCCTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATC  >  NC_002947/5119094‑5119390

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: ‑