Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,161,299 T→A 6.0% intergenic (+47/‑169) PP_1013 → / → yihS sensor histidine kinase/aldose‑ketose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,161,2990TA 6.0% 888.8 / 5.2 308intergenic (+47/‑169)PP_1013/yihSsensor histidine kinase/aldose‑ketose isomerase
Reads supporting (aligned to +/- strand):  ref base T (119/155);  new base A (9/13);  total (137/171)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 323 total aligned reads displayed.

TGCGCGCAACATCGCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTG  >  NC_002947/1161150‑1161448
                                                                                                                                                      |                                                                                                                                                     
tGCGCGCAACATCGCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCt                                                                                                                                                       <  2:43988/150‑1 (MQ=255)
tGCGCGCAACATCGCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCt                                                                                                                                                       <  2:4277485/150‑1 (MQ=255)
tGCGCGAAACATCGCCCATAGCCATGGGGGCGTGGTGAGCTTGCTTAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCt                                                                                                                                                       <  2:3038783/150‑1 (MQ=255)
 gcgcgcAACATCGCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTTGGGTGCCTGTGCGGGCCTCTTCGCGGGTGAACCCGCTCCTa                                                                                                                                                      <  1:2088182/150‑1 (MQ=255)
   gcgcAACATCGCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACa                                                                                                                                                    >  2:8028045/1‑150 (MQ=255)
    cgcAACATCGCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACaa                                                                                                                                                   >  2:4784677/1‑150 (MQ=255)
    cgcAACATCGCCCATAGCAATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACaa                                                                                                                                                   <  2:4199947/150‑1 (MQ=255)
     gcAACATCGCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAAc                                                                                                                                                  >  1:8654843/1‑150 (MQ=255)
     gcAACATCGCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAAc                                                                                                                                                  <  1:1826274/150‑1 (MQ=255)
     gcAACATCGCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAAc                                                                                                                                                  <  1:2021602/150‑1 (MQ=255)
        aCATCGCCCTTAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGAc                                                                                                                                               <  2:3624399/150‑1 (MQ=255)
        aCATCGCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGAc                                                                                                                                               <  2:3752914/150‑1 (MQ=255)
          aTCGCCCTTAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCCTCTTCGCGGGTGAACCCGCTCCTacgacgcaaa                                                                                                                                             <  2:4530031/150‑5 (MQ=255)
           tCGCCCATAGCCTTGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGa                                                                                                                                            <  2:6707065/150‑1 (MQ=255)
           tCGCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGa                                                                                                                                            <  2:607101/150‑1 (MQ=255)
           tCGCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGa                                                                                                                                            <  1:4988813/150‑1 (MQ=255)
            cGCCCATCGCCTTGGGGGCGTGGTGAGCTTGCTGAACCTGGGGGTGGGGGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAg                                                                                                                                           <  1:785716/150‑1 (MQ=38)
             gCCAATAGCAATGGGGGCGTGGTGAGGTTGCTGAACCTGGGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGt                                                                                                                                          <  1:1037068/150‑1 (MQ=39)
                cATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAg                                                                                                                                       >  2:4004236/1‑150 (MQ=255)
                 attgcaaTGGGGGCGTGGTTAGGTTGCTTAACGTGGGGGTGGGGGGGCTGGGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAgg                                                                                                                                      <  1:3866195/144‑1 (MQ=37)
                    gCCATGGGGGCGAGGTGCGCTTGCGGAACCTGCGGGTGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGttt                                                                                                                                   <  2:4968273/150‑1 (MQ=255)
                    gCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGttt                                                                                                                                   <  1:5078156/150‑1 (MQ=255)
                      cATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGc                                                                                                                                 <  2:2053350/150‑1 (MQ=255)
                        tGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCACCAACAACGACCCAGTCAGGTTTGcgt                                                                                                                               >  1:3593998/1‑150 (MQ=255)
                        tGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAAACCCCACCAAAAACCACCCAGTCAGGTTTgggg                                                                                                                               >  1:6881451/1‑147 (MQ=255)
                         gggggCGAGTTGAGCTTGCTGAACCTGCGGGAGGGCGGGCGGGGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGcgtc                                                                                                                              <  1:7371749/150‑1 (MQ=255)
                         gggggCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGcgtc                                                                                                                              <  2:3321296/150‑1 (MQ=255)
                         gggggCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGcgtc                                                                                                                              <  2:3938636/150‑1 (MQ=255)
                           gggCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGcgtcgt                                                                                                                            <  2:1691365/150‑1 (MQ=255)
                            ggCGCGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTa                                                                                                                           <  2:4221970/150‑1 (MQ=255)
                            ggCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTa                                                                                                                           <  2:7320335/150‑1 (MQ=255)
                            ggCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTa                                                                                                                           <  2:5365462/150‑1 (MQ=255)
                             gCGAGGTGCGCTTGCTGAACCTGCGGGAGGGCGGGCGGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAg                                                                                                                          <  1:6759737/150‑1 (MQ=255)
                             gCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGTACCCCCGCCAACAACGACCCAGGGAGGGTTGcgtcgtgg                                                                                                                          >  1:6794159/1‑148 (MQ=255)
                               gAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGa                                                                                                                        <  2:6505865/150‑1 (MQ=255)
                               gAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGa                                                                                                                        <  1:2258711/150‑1 (MQ=255)
                               gAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGa                                                                                                                        <  2:1850476/150‑1 (MQ=255)
                               gAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGa                                                                                                                        <  1:4784677/150‑1 (MQ=255)
                               gAGGTGAGCTGGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGa                                                                                                                        <  1:8380410/150‑1 (MQ=255)
                                aGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAg                                                                                                                       <  2:7894737/150‑1 (MQ=255)
                                 ggggAGGTTGCTGACCCGGCGGGGGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGc                                                                                                                      <  1:7839786/147‑1 (MQ=39)
                                   tGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCgg                                                                                                                    <  1:1298311/150‑1 (MQ=255)
                                   tGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCgg                                                                                                                    <  2:2682969/150‑1 (MQ=255)
                                    gcgCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCggg                                                                                                                   <  2:5565201/148‑1 (MQ=255)
                                        ttGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCGTCGCGGGTGAAACCCCGCCGAAAAACACCCAGTCAAGTTTGGGGGGTAGGTGGGGGttga                                                                                                               >  1:248492/1‑148 (MQ=255)
                                        gtGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACACCGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCa                                                                                                               <  2:6333781/149‑1 (MQ=255)
                                         tGCTGAACCTGCGGGCGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCAc                                                                                                              <  1:5397133/150‑1 (MQ=255)
                                         tGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCAc                                                                                                              <  2:5623580/150‑1 (MQ=255)
                                              aaCGTGGGGGAGGGGGGGCTGCGGGTGACCCTGTATCTGCCGGGGGGTTTGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGCCCGATTCTGGTTTGCGTCGTAGGAGCGGGTTCACCcgcg                                                                                                         <  2:6125929/150‑1 (MQ=38)
                                              aaCCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGTAACCGCGCCAACAACGACCGAGTCAGGTTTGCGGAGTAGGAGCGGGTTGAcccccc                                                                                                         >  1:1298536/1‑147 (MQ=255)
                                                ccTGGGGGGGGGCGGCCTGGGGGTCACCCTGTAGCTGCCGGGGGGTAGGGATTGTGTTTTGGGGGGGCCTGTGCGGGCTTCTTCCGGGTTGAACCCGCGCCTACAACCACCGAGTCTGGTTTGCGTCTTTGGAGCGGGTTCACCCGCGaa                                                                                                       <  2:287382/150‑1 (MQ=255)
                                                  tGCGGTTGGGGGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGCGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                     <  2:1298586/150‑1 (MQ=255)
                                                   gCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAagag                                                                                                    <  2:8818129/150‑1 (MQ=255)
                                                      gttggggggggTGGGGGTGACCCTTTATCTGCGGGGGGTTTTGTTTTTTGTTTTTGGGGTGCCTGTGCGGGCTTCTTGGGGGGTTAACCCGCTCCTACAACGACCTTGACAGGTTTGCGTGGTTGGTGCGGGTTCACCCGCGAAGAGGcc                                                                                                 <  2:898576/147‑1 (MQ=255)
                                                           gcgggcTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGcac                                                                                            >  2:5753311/1‑150 (MQ=255)
                                                            cgggcTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGGTGGAGCGGGTTCACCCGCGACGAGGGCCGcacc                                                                                           >  2:6613587/1‑149 (MQ=255)
                                                                      gtcaCCCTTTATCTGCCGCGGGGTATGTATTGATTTTTTGGGGTGCCTGTGCGGGCTTATTCGCGGTTCAACCCGCACCAACACCCACCGAGCCAGGTTTGCGCCGTAGCAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCaaa                                                                                 <  2:4498265/147‑1 (MQ=37)
                                                                        gACCCTGTATCTGCCGCGGGTTATGGATTGTGTTTTGGGGGGGCCTGGGCGGGCTTCTTCGCGGGTGAACCCGCTCCTCCAACGCCCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAAt                                                                               <  2:7930982/150‑1 (MQ=37)
                                                                          ccgtGTATGGGCCGGGGGGTAGGGTTTGTGTTTTGGGGGTGCCGGGGGGGGCTGCTGCGCGGGTCACCCCGCACCAACACCGACCGAGTCTGGTGGGGGTGGTCGGCGGGGTTCCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGt                                                                             <  2:4137093/147‑1 (MQ=255)
                                                                          cccTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGGGCACCCCCCCCGACAACCACCGCGGGACGGTTGCGGCGGAGGAGCGGGGTGACCCCCGAAGAGGCCCCCCCACGCCCCCAAAAAGt                                                                             >  1:4204136/1‑150 (MQ=255)
                                                                            cTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACACCGACCGCGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCa                                                                           <  2:5839193/150‑1 (MQ=255)
                                                                            cTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAAACCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAAAGTCa                                                                           >  2:3390160/1‑150 (MQ=255)
                                                                            cTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAAACCGCGCCTACAACGACCGAGTCAGGTTTGCGGGGGAGGAGCGGGTTCAACCGCGAAGAGGCCCGCACAGGCCCCCAAAATGtaa                                                                           >  2:1215634/1‑148 (MQ=255)
                                                                             tGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCAc                                                                          >  2:8643266/1‑150 (MQ=255)
                                                                               tATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTg                                                                        >  2:7140025/1‑150 (MQ=255)
                                                                                             tATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCACGTTTTCGTCGTAGGAGCCGGTTCACCCCCCGACACGCCCCCACACGCACCCAAAAAGTCACAGCCCCGTGACaaaac                                                          >  2:6678766/1‑147 (MQ=25)
                                                                                                 gATTGAGTTTTGGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGACCCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCcgtc                                                      <  2:7233179/150‑1 (MQ=255)
                                                                                                        tggggggggggCCTGGGCGGGCGTCGGCGCGGGGGACCCCGCGCCTCCCCCGCCCGCGGCCGGTTTGCGTCGTCGGCGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTc                                               <  2:6370501/145‑1 (MQ=255)
                                                                                                         ttttGGGGGGCCGGGTGCGGGATTCTCGGGGGGGCACCCCGCCCATCCACCGACCGCGTCTGGTGGGGGTCGTTGGCGGGGGCCCCCCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTc                                               <  1:6246832/150‑1 (MQ=255)
                                                                                                               ggTGCCTGTGCGGGCTTCTTCGCGGGTTAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACAt                                        <  2:5035748/150‑1 (MQ=255)
                                                                                                                   ccTGTGGGGGATACTTCGCGGGTGAACCCGCCCCACCCACGCCCGAGTCTGGTTTGCGTCGCCGGCGGGGGATCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGc                                    <  1:1391078/150‑1 (MQ=255)
                                                                                                                           gggCTTCTTCGCGGGTGTACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGAcc                             >  2:7520285/1‑150 (MQ=255)
                                                                                                                           gggCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTAACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGAcc                             >  2:2522775/1‑150 (MQ=255)
                                                                                                                           gggATTCTGCGGGGGTCACCCCCCACCAACACCCCCCGTGTCTGTTTGCCGTCGTGGGCGGGGGTCCCCCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGAcc                             <  1:7531504/150‑1 (MQ=255)
                                                                                                                                 cttcGCGGGTGAACCCGCCCCACCACCGCCGGTGTCTGGTTTGCGTCGTAGGCGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTg                       <  1:8028045/150‑1 (MQ=38)
                                                                                                                                     gcgGGTGAACCCGCTCCGACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAAAACCGTCGTCCTTCGTGACATCCGCCAAGACCCCGCTGGCTa                   >  1:3072168/1‑150 (MQ=255)
                                                                                                                                     gcgGGGCAACCCGCTCCTACACCGCCCGCGTCAGGTTTGCGTCGTCGGCGCGGGTTCCCCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTa                   <  2:1466330/150‑1 (MQ=255)
                                                                                                                                       gggTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACCATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTCGa                 >  2:6185739/1‑150 (MQ=255)
                                                                                                                                       gggTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGa                 >  2:2341520/1‑150 (MQ=255)
                                                                                                                                         gTGACCCCGCTCCTACACCGCCCGTGTCTGTTTTGCGTCGTAGTAGCGGGTCCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAAt               <  1:6498658/150‑1 (MQ=255)
                                                                                                                                         gTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAAt               >  2:6898405/1‑150 (MQ=255)
                                                                                                                                          tGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAAtt              >  1:2596237/1‑150 (MQ=255)
                                                                                                                                           gAACCCGCTCCAACACCCCCCGAGTCTGGTTTGCGTCGTAGGAGGGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTc             <  2:8153315/150‑1 (MQ=39)
                                                                                                                                             aCCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGc           >  2:6737139/1‑150 (MQ=255)
                                                                                                                                              cccGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCa          >  2:7644200/1‑150 (MQ=255)
                                                                                                                                              cccGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCa          >  2:5003436/1‑150 (MQ=255)
                                                                                                                                               ccGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCaa         >  1:2888462/1‑150 (MQ=255)
                                                                                                                                               ccGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCaa         >  1:5766748/1‑150 (MQ=255)
                                                                                                                                               ccGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCaa         >  1:5101533/1‑150 (MQ=255)
                                                                                                                                               ccGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCaa         >  2:4639986/1‑150 (MQ=255)
                                                                                                                                               ccGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCaa         >  1:4043036/1‑150 (MQ=255)
                                                                                                                                                   tCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCAGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGa     >  1:6012710/1‑150 (MQ=255)
                                                                                                                                                     cTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACt   >  2:5106740/1‑150 (MQ=255)
                                                                                                                                                     cTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACt   >  1:3366405/1‑150 (MQ=255)
                                                                                                                                                     cTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACt   >  1:704804/1‑150 (MQ=255)
                                                                                                                                                      tACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTg  >  1:5265098/1‑150 (MQ=255)
                                                                                                                                                      tACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTg  >  2:3003112/1‑150 (MQ=255)
                                                                                                                                                      tACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTg  >  2:6429992/1‑150 (MQ=255)
                                                                                                                                                      tACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTg  >  2:8423480/1‑150 (MQ=255)
                                                                                                                                                      tACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTg  >  2:1497425/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TGCGCGCAACATCGCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTG  >  NC_002947/1161150‑1161448

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