Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,123,192 T→G 5.8% intergenic (+27/‑32) PP_4504 → / → PP_4505 hypothetical protein/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,123,1920TG 5.8% 620.8 / 18.7 183intergenic (+27/‑32)PP_4504/PP_4505hypothetical protein/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (89/82);  new base G (5/6);  total (95/88)
Fisher's exact test for biased strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 185 total aligned reads displayed.

AAGTTGCTCACGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGACGCCTGGATCCTCA  >  NC_002947/5123043‑5123341
                                                                                                                                                     |                                                                                                                                                     
aaGTTGCTCACGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCAt                                                                                                                                                       <  1:392340/150‑1 (MQ=255)
  gTTGCTCACGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCa                                                                                                                                                     >  2:5112385/1‑150 (MQ=255)
     gCTCACGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACtt                                                                                                                                                  <  1:1240363/150‑2 (MQ=255)
      cTCACGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGc                                                                                                                                                 <  1:4939368/150‑1 (MQ=255)
       tCACGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGcg                                                                                                                                                >  1:2309810/1‑150 (MQ=255)
       tCACGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGcg                                                                                                                                                >  2:386002/1‑150 (MQ=255)
        cACGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGcgg                                                                                                                                               >  1:3349507/1‑150 (MQ=255)
             cTGACCAGGAGAAGAATTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAACTCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCGCGGCCCGCAGTGATGcgggccgg                                                                                                                                          <  2:3081501/150‑1 (MQ=255)
             cTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCAGTGATGcgggccgg                                                                                                                                          <  2:3897679/150‑1 (MQ=255)
             cTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCAGTGATGcgggccgg                                                                                                                                          <  2:3081444/150‑1 (MQ=255)
                aCCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCt                                                                                                                                       <  2:4940033/150‑1 (MQ=255)
                 ccAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCtt                                                                                                                                      <  1:5371617/150‑1 (MQ=255)
                    gcagAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTc                                                                                                                                   <  2:3621341/150‑1 (MQ=255)
                     cagAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTcc                                                                                                                                  >  1:1038250/1‑150 (MQ=255)
                       gaagaaGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTATCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCg                                                                                                                                <  1:4043415/150‑1 (MQ=255)
                          gaagTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGttt                                                                                                                             >  1:3141308/1‑150 (MQ=255)
                           aagTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTa                                                                                                                            <  1:3085795/150‑1 (MQ=255)
                           aagTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTa                                                                                                                            <  2:5451816/150‑1 (MQ=255)
                              tttGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTagga                                                                                                                         >  2:5482016/1‑150 (MQ=255)
                               ttGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTaggag                                                                                                                        <  1:4923710/150‑1 (MQ=255)
                                  aCGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTg                                                                                                                     <  2:796524/150‑1 (MQ=255)
                                    gAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCa                                                                                                                   >  2:4917172/1‑150 (MQ=255)
                                     aaCTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCAc                                                                                                                  <  2:2375216/150‑1 (MQ=255)
                                       cTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGCAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACtt                                                                                                                <  2:3809329/150‑1 (MQ=255)
                                       cTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACtt                                                                                                                <  2:2817356/150‑1 (MQ=255)
                                        tGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGc                                                                                                              >  1:1442104/1‑150 (MQ=255)
                                         gCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGc                                                                                                              <  1:458308/150‑1 (MQ=255)
                                         gCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGGTTAGGAGGTGCACTTGc                                                                                                              >  1:3946226/1‑150 (MQ=255)
                                           aaaaaGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGt                                                                                                            >  2:5189918/1‑150 (MQ=255)
                                           aaaaaGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGt                                                                                                            >  1:590081/1‑150 (MQ=255)
                                           aaaaaGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGt                                                                                                            >  2:5191402/1‑150 (MQ=255)
                                           aaaaaGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCAACACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGt                                                                                                            >  1:4076393/1‑150 (MQ=255)
                                            aaaaGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGtt                                                                                                           <  2:3475326/150‑1 (MQ=255)
                                            aaaaGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGtt                                                                                                           <  1:4824547/150‑1 (MQ=255)
                                                gAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACt                                                                                                       <  1:1158409/150‑1 (MQ=255)
                                                   cgcgCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGtt                                                                                                    <  1:2796887/150‑1 (MQ=255)
                                                    gcgcCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTc                                                                                                   <  2:844374/150‑1 (MQ=255)
                                                      gccgccAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTg                                                                                                 >  2:4557753/1‑150 (MQ=255)
                                                       ccgccAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTgg                                                                                                <  2:1599608/150‑1 (MQ=255)
                                                        cgccTAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGgc                                                                                               <  2:3204814/150‑1 (MQ=255)
                                                          ccAAGACAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGgcgc                                                                                             <  1:1136927/150‑1 (MQ=255)
                                                          ccAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCTCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGgcgc                                                                                             <  1:4533353/150‑1 (MQ=255)
                                                           cAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGTCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCa                                                                                            >  1:776932/1‑150 (MQ=255)
                                                           cAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCa                                                                                            <  1:981519/150‑1 (MQ=255)
                                                            aaGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCAt                                                                                           >  1:2433391/1‑150 (MQ=255)
                                                            aaGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGGTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCAt                                                                                           >  1:2893010/1‑150 (MQ=255)
                                                              gAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATcc                                                                                         >  1:2934313/1‑150 (MQ=255)
                                                                  gCCGAGTGGCCGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGc                                                                                     <  2:4299662/150‑1 (MQ=255)
                                                                  gCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGc                                                                                     <  2:4016877/150‑1 (MQ=255)
                                                                   ccGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGcc                                                                                    >  1:2283027/1‑150 (MQ=255)
                                                                    cGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCa                                                                                   <  1:218694/150‑1 (MQ=255)
                                                                        tGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGttt                                                                               <  2:5463214/150‑1 (MQ=255)
                                                                          gCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCAAGCCCGGCCCGCAGTGATGcgggccgggcat                                                                                                                                      >  2:4096486/1‑91 (MQ=37)
                                                                          gCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCAAGCCCGGCCCGCAGTGATGcgggccgggcat                                                                                                                                      <  1:4096486/93‑3 (MQ=37)
                                                                           cAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTctg                                                                            >  1:1899104/1‑150 (MQ=255)
                                                                           cAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTctg                                                                            <  1:3894250/150‑1 (MQ=255)
                                                                            aGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTctgg                                                                           >  1:877629/1‑150 (MQ=255)
                                                                            aGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTctgg                                                                           <  2:1414136/150‑1 (MQ=255)
                                                                            aGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTctgg                                                                           >  2:3396526/1‑150 (MQ=255)
                                                                             ggAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTctggc                                                                          >  1:2988855/1‑150 (MQ=255)
                                                                             ggAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGGTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTctggc                                                                          >  2:502560/1‑150 (MQ=255)
                                                                              gAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTctggct                                                                         >  2:1692487/1‑150 (MQ=255)
                                                                              gAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTctggct                                                                         >  2:2592561/1‑150 (MQ=255)
                                                                                 ttCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTGAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCtggctggc                                                                      >  2:4866098/1‑150 (MQ=255)
                                                                                   ccAGGCCTGGACAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCCGGCTTTGCCCGGGTAGGAGGTGCACTTGCGTTCACTGGTCTGGCGCATCCTGGCCAGTTTCGGGCTGGCCa                                                                    >  2:4250840/1‑150 (MQ=255)
                                                                                     aGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATc                                                                  >  2:939538/1‑150 (MQ=255)
                                                                                     aGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATc                                                                  >  1:5389573/1‑150 (MQ=255)
                                                                                         cTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCaggatgc                                                              <  1:2661612/150‑7 (MQ=255)
                                                                                         cTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCAccc                                                              >  1:2050408/1‑150 (MQ=255)
                                                                                          tGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCt                                                             >  2:183860/1‑150 (MQ=255)
                                                                                          aggAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGTGATTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCt                                                             <  2:5386791/149‑1 (MQ=255)
                                                                                            gAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTgg                                                           >  2:2410542/1‑150 (MQ=255)
                                                                                                 cTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCa                                                      >  1:1555670/1‑150 (MQ=255)
                                                                                                 cTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCa                                                      >  1:3112467/1‑150 (MQ=255)
                                                                                                 cTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCa                                                      >  2:4794118/1‑150 (MQ=255)
                                                                                                 cTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCa                                                      >  1:4759341/1‑150 (MQ=255)
                                                                                                  tGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAg                                                     <  1:5333985/150‑1 (MQ=255)
                                                                                                  tGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACGTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAg                                                     >  2:4613709/1‑150 (MQ=255)
                                                                                                           cgcCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGGTTAGTAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTcta                                            <  1:5112307/150‑3 (MQ=255)
                                                                                                              cAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATcc                                         >  2:4819619/1‑150 (MQ=255)
                                                                                                              cAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATcc                                         >  1:5584553/1‑150 (MQ=255)
                                                                                                                   cccAGTAAGCTGTCCGCCCCAGGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCctgctg                                    <  2:5583588/150‑1 (MQ=255)
                                                                                                                    ccAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTgg                                   <  2:1671701/150‑1 (MQ=255)
                                                                                                                        tAAGCTGTCCGCCCCATGCCCGGCCCGCATCACAGCCGGCCGGGCTTTACCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCGGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCGGGGGcc                               >  1:3490270/1‑149 (MQ=255)
                                                                                                                         aaGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGcac                              <  1:2386752/150‑1 (MQ=255)
                                                                                                                         aaGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGcac                              <  1:497807/150‑1 (MQ=255)
                                                                                                                             tGTCCGCCCCATGCCCGGCCCGCATCAATGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTAGCAGGCCTGTCGATCCAGCTGGGGCACAAGc                          >  1:4229471/1‑150 (MQ=255)
                                                                                                                               tCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTc                        <  1:473074/150‑1 (MQ=255)
                                                                                                                                 cGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCaa                      >  1:4442577/1‑150 (MQ=255)
                                                                                                                                  gCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAAc                     >  1:4719758/1‑150 (MQ=255)
                                                                                                                                    cccATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCAGTCGATCCTGCTGGGGCACATGCTCAACCa                   >  1:2608385/1‑150 (MQ=255)
                                                                                                                                     ccATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAg                  <  1:4018093/150‑1 (MQ=255)
                                                                                                                                      cATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAgg                 >  1:431201/1‑150 (MQ=255)
                                                                                                                                          cccggcccgCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGACGc             >  1:845691/1‑150 (MQ=255)
                                                                                                                                          cccggcccgCAGTGATGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGACGc             >  2:2106058/1‑150 (MQ=255)
                                                                                                                                           ccggcccgCAGTGATGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGACGcc            >  1:4154512/1‑150 (MQ=255)
                                                                                                                                           ccggcccgCAGTGATGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGACGcc            >  1:5460725/1‑150 (MQ=255)
                                                                                                                                                 cgCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGACGCCTGGATc      >  2:4955265/1‑150 (MQ=255)
                                                                                                                                                 cgCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGACGCCTGGATc      <  1:3433098/150‑1 (MQ=255)
                                                                                                                                                     tCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACGTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGACGCCTGGATCCTCa  <  2:3337798/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGTTGCTCACGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGACGCCTGGATCCTCA  >  NC_002947/5123043‑5123341

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