Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,354,147 C→T 5.3% intergenic (+72/+70) PP_2947 → / ← PP_2948 transcriptional regulator MvaT/GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,354,1470CT 5.3% 534.5 / 12.0 169intergenic (+72/+70)PP_2947/PP_2948transcriptional regulator MvaT/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (84/76);  new base T (5/4);  total (89/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 187 total aligned reads displayed.

GGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACGGAGGAAATCAC  >  NC_002947/3354000‑3354296
                                                                                                                                                   |                                                                                                                                                     
gggTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGcctc                                                                                                                                                     >  1:136484/1‑150 (MQ=255)
  gTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGcctcct                                                                                                                                                   <  2:1238905/150‑1 (MQ=255)
   tGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAGGCGACGCTGTCGCttcctc                                                                                                                                                  <  1:2956416/150‑1 (MQ=255)
   tGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCctcctc                                                                                                                                                  <  2:2177390/150‑1 (MQ=255)
   tGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCctcctc                                                                                                                                                  <  1:4967441/150‑1 (MQ=255)
   tGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCctcctc                                                                                                                                                  <  1:2752456/150‑1 (MQ=255)
       aaCCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGaaa                                                                                                                                              >  2:2814950/1‑150 (MQ=255)
       aaCCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTATTTTTCGAGGACGCGACAGCGTCGCCTCCTCGaaa                                                                                                                                              >  2:88253/1‑150 (MQ=255)
         ccACAAGCTGCTCAAGGACTGGCACGCCCAGTACGGGGCAGAAGAAGTGGAAGCCTGGGTGTCCTGAGTGGTGCTGCTGTCGCCCTCAGTGCGCCATCAGAAACGCCTTGTTCTTTTTCGAGGAGGCGACGCTGTCGCTTCCTCGaaaaa                                                                                                                                            <  2:149175/150‑3 (MQ=255)
          cacaAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAg                                                                                                                                           >  1:2458634/1‑150 (MQ=255)
          cacaAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAg                                                                                                                                           >  1:4135405/1‑150 (MQ=255)
           acaAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAgg                                                                                                                                          <  2:5211498/150‑1 (MQ=255)
           acaAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAgg                                                                                                                                          <  1:347294/150‑1 (MQ=255)
             aaGCTGCACAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCg                                                                                                                                       <  1:3718665/150‑1 (MQ=255)
              aGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCg                                                                                                                                       <  2:607598/150‑1 (MQ=255)
               gctgctCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGt                                                                                                                                      >  1:945692/1‑150 (MQ=255)
               gctgctCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGt                                                                                                                                      <  2:2286731/150‑1 (MQ=255)
                 tgctCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAAAGCCTTGTTCTTTTTCGAGGAAGAGACAGCGTCGCCTCCTCGAAACAGGGCGTTc                                                                                                                                    >  2:460155/1‑150 (MQ=255)
                    tCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATc                                                                                                                                 >  1:2781365/1‑150 (MQ=255)
                      aaGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAACCAGGGCGTTCATCgt                                                                                                                               >  1:4811770/1‑150 (MQ=255)
                       aGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgtg                                                                                                                              <  1:1804586/150‑1 (MQ=255)
                              gAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGc                                                                                                                       >  2:2929651/1‑150 (MQ=255)
                                    ccAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATg                                                                                                                 >  1:4613434/1‑150 (MQ=255)
                                    ccAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATg                                                                                                                 >  2:4803919/1‑150 (MQ=255)
                                     cAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGa                                                                                                                <  1:2313332/150‑1 (MQ=255)
                                      aGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAg                                                                                                               <  1:4183374/150‑1 (MQ=255)
                                           ggggCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGAc                                                                                                          >  2:2069417/1‑150 (MQ=255)
                                            gggCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACg                                                                                                         <  1:5421972/150‑1 (MQ=255)
                                             ggCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGt                                                                                                        >  1:467463/1‑150 (MQ=255)
                                             ggCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGt                                                                                                        <  1:4950826/150‑1 (MQ=255)
                                              gCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTg                                                                                                       <  1:1268414/150‑1 (MQ=255)
                                              gCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTg                                                                                                       >  1:1565391/1‑150 (MQ=255)
                                                agaagaAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCTTTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGcctc                                                                                                                                                     <  1:10119/102‑1 (MQ=255)
                                                agaagaAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGcctc                                                                                                                                                     >  2:10119/1‑102 (MQ=255)
                                                agaagaAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGt                                                                                                     <  2:5239709/150‑1 (MQ=255)
                                                     aagTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAAc                                                                                                <  2:2139636/150‑1 (MQ=255)
                                                      agTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACa                                                                                               <  1:4404962/150‑1 (MQ=255)
                                                        tGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAgg                                                                                             >  1:583752/1‑150 (MQ=255)
                                                          gtttgcTGTGTGTCCTGATTTGTGGTGATGTGGCCGTCTTTGCTGTTTGTGATACGGCTTGTTCTTTTTGGTGGATGGGTGTGGGTGGCCTCCTCGAAACTGGGCGTTGTTGGTTGGTGGGCGGAATGAGCTGAGGTGGTTTATCTGGaa                                                                                           <  2:3312595/145‑1 (MQ=255)
                                                          gAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAAc                                                                                                <  1:4901936/145‑1 (MQ=255)
                                                          gAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAAc                                                                                                >  2:4901936/1‑145 (MQ=255)
                                                          gAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGaa                                                                                           >  1:2795761/1‑150 (MQ=255)
                                                             tCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGtt                                                                                        >  2:4013799/1‑150 (MQ=255)
                                                               cTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTc                                                                                      <  2:4420757/150‑1 (MQ=255)
                                                                 gggTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGc                                                                                    <  2:945692/150‑1 (MQ=255)
                                                                         tGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAgc                                                                            >  2:767516/1‑150 (MQ=255)
                                                                         tGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAgc                                                                            >  2:2385718/1‑150 (MQ=255)
                                                                         tGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAgc                                                                            >  1:3611946/1‑150 (MQ=255)
                                                                         tGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAgc                                                                            >  2:3622207/1‑150 (MQ=255)
                                                                         tGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAgc                                                                            >  1:3004406/1‑150 (MQ=255)
                                                                         tGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAgc                                                                            >  2:44896/1‑150 (MQ=255)
                                                                         tGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAgc                                                                            >  1:1449205/1‑150 (MQ=255)
                                                                                tgCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAgt                                                                     >  2:4995064/1‑150 (MQ=255)
                                                                                tgCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAgt                                                                     >  2:2534879/1‑150 (MQ=255)
                                                                                 gCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAgtg                                                                    >  2:1659211/1‑150 (MQ=255)
                                                                                    gATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTcc                                                                 <  1:3323346/150‑1 (MQ=255)
                                                                                      tGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCAACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAATCAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGt                                                               >  2:4395273/1‑150 (MQ=255)
                                                                                        tCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGttt                                                             <  1:296073/150‑1 (MQ=255)
                                                                                            ccTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGaa                                                                                           <  2:5247561/116‑1 (MQ=255)
                                                                                             cTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATa                                                        >  1:3953381/1‑150 (MQ=255)
                                                                                              tCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATAc                                                       >  2:1465773/1‑150 (MQ=255)
                                                                                                aTTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACtt                                                     <  2:4091906/150‑1 (MQ=255)
                                                                                                     gcTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAggg                                                                                                                                         <  1:2167582/61‑1 (MQ=255)
                                                                                                       tATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGATGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGGAAAt                                              >  2:3841744/1‑150 (MQ=255)
                                                                                                       tATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAAt                                              >  1:840365/1‑150 (MQ=255)
                                                                                                       tATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAAt                                              >  1:654089/1‑150 (MQ=255)
                                                                                                       tATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAAt                                              >  1:4415753/1‑150 (MQ=255)
                                                                                                           tGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTg                                          >  2:4317343/1‑150 (MQ=255)
                                                                                                           tGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTg                                          >  1:1813136/1‑150 (MQ=255)
                                                                                                           tGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTg                                          <  1:2069417/150‑1 (MQ=255)
                                                                                                             aaaCGCCTTGTTCATTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGATTCAATACTTCGTAAATTATGGc                                        >  2:4652839/1‑150 (MQ=255)
                                                                                                                 cccTTGTACTATATCGACCAAGCGACAGCGTCGCCTCCACGAAACAGGCAGTTCATCGTCCGTGCGCCCAATGAGCTGACGTCGATAAACACGAAGTAACGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATATTGGCATCa                                    <  2:5075957/149‑1 (MQ=255)
                                                                                                                   cTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTaaa                                                                                                 >  1:54885/1‑87 (MQ=255)
                                                                                                                   cTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTaaa                                                                                                 <  2:54885/87‑1 (MQ=255)
                                                                                                                   cTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTg                                          <  1:1055292/142‑1 (MQ=255)
                                                                                                                   cTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTg                                          >  2:1055292/1‑142 (MQ=255)
                                                                                                                       ttCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGAt                              >  1:689135/1‑150 (MQ=255)
                                                                                                                         cTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTg                            <  2:2177381/150‑1 (MQ=255)
                                                                                                                           ttttCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTggg                          >  2:3073233/1‑150 (MQ=255)
                                                                                                                            tttCGAGGAGGCGACGCTGTCGCTTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGc                         >  2:4911115/1‑150 (MQ=255)
                                                                                                                              tCGAGGAGGCGACGCTGTCGCTTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCGCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGa                       >  1:3886845/1‑150 (MQ=255)
                                                                                                                                gaggagGCGACGCTGTCGCTTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCa                                    >  2:3985234/1‑135 (MQ=255)
                                                                                                                                gaggagGCGACGCTGTCGCTTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCa                                    <  1:3985234/135‑1 (MQ=255)
                                                                                                                                     aGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGaa                <  2:4191583/150‑1 (MQ=255)
                                                                                                                                      gCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGaaa               <  1:2218697/150‑1 (MQ=255)
                                                                                                                                      gCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGaaa               <  1:3681402/150‑1 (MQ=255)
                                                                                                                                      gCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGaaa               >  1:3850971/1‑150 (MQ=255)
                                                                                                                                        gACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACg             >  2:3486148/1‑150 (MQ=255)
                                                                                                                                          cAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACgga           <  2:1937852/150‑1 (MQ=255)
                                                                                                                                          cAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACgga           <  1:938658/150‑1 (MQ=255)
                                                                                                                                          cAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACgga           <  1:3508667/150‑1 (MQ=255)
                                                                                                                                          cAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACgga           <  1:920999/150‑1 (MQ=255)
                                                                                                                                          cAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACgga           <  2:736073/150‑1 (MQ=255)
                                                                                                                                          cAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACgga           <  2:330640/150‑1 (MQ=255)
                                                                                                                                           aGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTTTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACggag          <  2:3774374/150‑1 (MQ=255)
                                                                                                                                           aGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACggag          <  2:136505/150‑1 (MQ=255)
                                                                                                                                               tCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACGGAGGaaa      <  1:1066755/150‑1 (MQ=255)
                                                                                                                                                  cctcctCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACGGAGGAAATCa   >  1:5118296/1‑150 (MQ=255)
                                                                                                                                                   ctcctcGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACGGAGGAAATCAc  <  2:1449205/150‑1 (MQ=255)
                                                                                                                                                   ctcctcGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACGGAGGAAATCAc  <  2:4178801/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACGGAGGAAATCAC  >  NC_002947/3354000‑3354296

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