Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,161,306 A→C 6.0% intergenic (+54/‑162) 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,3060AC 6.0% 1026.1 / 8.3 332intergenic (+54/‑162)PP_1013/yihSsensor histidine kinase/aldose‑ketose isomerase
Reads supporting (aligned to +/- strand):  ref base A (170/130);  new base C (3/26);  total (175/157)
Fisher's exact test for biased strand distribution p-value = 9.13e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 334 total aligned reads displayed.

ACATCGCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTGGCAAGGA  >  NC_002947/1161158‑1161455
                                                                                                                                                    |                                                                                                                                                     
aCATCGCCCTTAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGAc                                                                                                                                                      <  2:3624399/150‑1 (MQ=255)
aCATCGCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGAc                                                                                                                                                      <  2:3752914/150‑1 (MQ=255)
   tCGCCCATTGCCTTGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGTGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGa                                                                                                                                                   <  1:8266248/150‑1 (MQ=255)
   tCGCCCATAGCCTTGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGa                                                                                                                                                   <  2:6707065/150‑1 (MQ=255)
   tCGCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGa                                                                                                                                                   <  1:5775178/150‑1 (MQ=255)
   tCGCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGa                                                                                                                                                   >  2:4321640/1‑150 (MQ=255)
     gCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGt                                                                                                                                                 >  2:5030317/1‑150 (MQ=255)
     gCCAATAGCAATGGGGGCGTGGTGAGGTTGCTGAACCTGGGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGt                                                                                                                                                 <  1:1037068/150‑1 (MQ=39)
       ccATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCAACAACGACCGAGTCa                                                                                                                                               >  1:4192467/1‑150 (MQ=255)
        cATAGCCTTGGGGGCGTGGTGAGCTTGCTGAACCTGGGGGTGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAg                                                                                                                                              <  1:3157146/150‑1 (MQ=255)
        cATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAg                                                                                                                                              >  2:4004236/1‑150 (MQ=255)
          tAGCCATGGGGGCGCGGTGAGCTTGCTGACCCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGt                                                                                                                                            <  1:5453702/150‑1 (MQ=255)
            gCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGTGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGttt                                                                                                                                          <  1:1828816/150‑1 (MQ=255)
              cttGGGGGCGTGGTGAGCTTGCTGAACCTGGGGGTGGGGGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGc                                                                                                                                        <  1:7432376/148‑1 (MQ=255)
                tGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGcgt                                                                                                                                      >  1:6263693/1‑150 (MQ=255)
                 gggggCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGcgtc                                                                                                                                     <  2:3321296/150‑1 (MQ=255)
                 gggggCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGcgtc                                                                                                                                     <  2:3938636/150‑1 (MQ=255)
                    ggCGCGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTa                                                                                                                                  <  2:4221970/150‑1 (MQ=255)
                    ggCGAGGTGAGCTTGCTGAACCTGCGGGTGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTa                                                                                                                                  <  1:7668771/150‑1 (MQ=255)
                    ggCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTa                                                                                                                                  <  1:719988/150‑1 (MQ=255)
                    ggCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTa                                                                                                                                  <  1:8594430/150‑1 (MQ=255)
                    ggCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTa                                                                                                                                  <  2:2344690/150‑1 (MQ=255)
                    ggCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTa                                                                                                                                  <  2:5045723/150‑1 (MQ=255)
                    ggCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTa                                                                                                                                  <  2:7320335/150‑1 (MQ=255)
                    ggCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTa                                                                                                                                  <  1:4028356/150‑1 (MQ=255)
                       gAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTTCGGGGTAGGa                                                                                                                               >  2:2962751/1‑150 (MQ=255)
                       gAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGa                                                                                                                               <  2:5318362/150‑1 (MQ=255)
                       gAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGa                                                                                                                               <  1:4784677/150‑1 (MQ=255)
                       gAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGa                                                                                                                               <  2:6505865/150‑1 (MQ=255)
                       gAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGa                                                                                                                               <  2:1850476/150‑1 (MQ=255)
                       gAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCCCTCCTACAACGACCGCGTCAGGTTTGGGTCGTAGGa                                                                                                                               >  2:1496014/1‑150 (MQ=255)
                         ggTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGc                                                                                                                             <  1:2081625/150‑1 (MQ=255)
                           tGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCgg                                                                                                                           <  2:5734306/150‑1 (MQ=255)
                           tGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCgg                                                                                                                           <  2:2558480/150‑1 (MQ=255)
                           tGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCgg                                                                                                                           <  1:1298311/150‑1 (MQ=255)
                            gAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCggg                                                                                                                          <  2:7522603/150‑1 (MQ=255)
                            gAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCggg                                                                                                                          <  1:2516516/150‑1 (MQ=255)
                            gAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGCGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCggg                                                                                                                          <  2:5565259/150‑1 (MQ=255)
                                ttGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCGTCGCGGGTGAAACCCCGCCGAAAAACACCCAGTCAAGTTTGGGGGGTAGGTGGGGGttga                                                                                                                      >  1:248492/1‑148 (MQ=255)
                                 tGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCAc                                                                                                                     <  2:44981/150‑1 (MQ=255)
                                   cTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCAccc                                                                                                                   <  1:6840504/150‑1 (MQ=255)
                                      aaCCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGTAACCGCGCCAACAACGACCGAGTCAGGTTTGCGGAGTAGGAGCGGGTTGAcccccc                                                                                                                >  1:1298536/1‑147 (MQ=255)
                                        ccTGGGGGGGGGCGGCCTGGGGGTCACCCTGTAGCTGCCGGGGGGTAGGGATTGTGTTTTGGGGGGGCCTGTGCGGGCTTCTTCCGGGTTGAACCCGCGCCTACAACCACCGAGTCTGGTTTGCGTCTTTGGAGCGGGTTCACCCGCGaa                                                                                                              <  2:287382/150‑1 (MQ=255)
                                          tgggggCGGGGGGGCTGGGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGGGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACCACCGAGTCAGGTTTGCGTCGTTGGAGCGGGTTCACCCGCGAaga                                                                                                            <  2:1298585/147‑1 (MQ=39)
                                           gCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAagag                                                                                                           <  2:8818129/150‑1 (MQ=255)
                                              gttggggggggTGGGGGTGACCCTTTATCTGCGGGGGGTTTTGTTTTTTGTTTTTGGGGTGCCTGTGCGGGCTTCTTGGGGGGTTAACCCGCTCCTACAACGACCTTGACAGGTTTGCGTGGTTGGTGCGGGTTCACCCGCGAAGAGGcc                                                                                                        <  2:898576/147‑1 (MQ=255)
                                                 gggcgggcTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGc                                                                                                     >  2:7824136/1‑150 (MQ=255)
                                                   gcgggcTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGcac                                                                                                   >  2:8338096/1‑150 (MQ=255)
                                                    cgggcTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGGTGGAGCGGGTTCACCCGCGACGAGGGCCGcacc                                                                                                  >  2:6613587/1‑149 (MQ=255)
                                                         ggCGGGTGACCGTGGACGGCCCGCGGGGTAGGGATGGTTTTTTGGGGGGGCCTGGGGGGGCTTCTGCGCGGGGCACCCCGCTCCTACACCGCCCGAGCCTGGTTTGCGGGTTTGGCGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCAc                                                                                             <  2:6154192/149‑1 (MQ=255)
                                                             ggTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCaa                                                                                         >  2:1318232/1‑150 (MQ=255)
                                                                gACCCTGTATCTGCCGCGGGTTATGGATTGTGTTTTGGGGGGGCCTGGGCGGGCTTCTTCGCGGGTGAACCCGCTCCTCCAACGCCCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAAt                                                                                      <  2:7930982/150‑1 (MQ=37)
                                                                  cccTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCCCTCCTACAACCACCGCGGCAGGGTTGCGGCGTAGGGGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCcaaaagg                                                                                    >  1:5002405/1‑147 (MQ=37)
                                                                   ccTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGGGAACCCGCGCCGACAACGACCGAGTCAGGTTTGCGTCGTAGGGGCGGGTTCACCCCCGAAGAGGCCCGCACAGGCACCCCAAATGTc                                                                                   >  2:4846065/1‑150 (MQ=255)
                                                                    cTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACACCGACCGCGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCa                                                                                  <  2:5839193/150‑1 (MQ=255)
                                                                    cTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAAACCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAAAGTCa                                                                                  >  2:3390160/1‑150 (MQ=255)
                                                                       tATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTg                                                                               >  2:2697047/1‑150 (MQ=255)
                                                                       tATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTg                                                                               >  2:7140025/1‑150 (MQ=255)
                                                                                   ggTAGGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAAt                                                                   <  2:4976239/150‑1 (MQ=255)
                                                                                     tATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGGGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGAGCCGGGTCACCCCCGCAGAGGCCCGCACACGCACCCAAAATGTCACTGCCCGGGGACAATccc                                                                >  2:6678698/1‑150 (MQ=38)
                                                                                         gATTGAGTTTTGGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGACCCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCcgtc                                                             <  2:7233179/150‑1 (MQ=255)
                                                                                           ttGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAAACCGCTCCTACAAAGAACGGGTCCGGTTTGGGTGGGAGGTGCGGGTTCACCCGCGAAGAGGCCCGCCCAGGCACCCAAAATGTCACTGCCCGGTGGCAATACcgtcgg                                                           >  2:8115755/1‑149 (MQ=255)
                                                                                                ttttgGGGGGGCCTGGGCGGGCTCCTGCGCGGGTCACCCCGCCCCTACACCCCCGGTGTCTGGTTTGCGTGGTCGGCGCGGTTTCCCCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTc                                                      <  1:3452627/150‑1 (MQ=255)
                                                                                                tttggggggTGCGTGTGCGGCCTTCGGGGCGGGTGCACCCGCGCCTCCCCCGCCGGCGTCCGGTTTGCGTCGTAGGAGCGGGTTCCCCGGCGAAGAGGCCCGCACAGGCACCCATAATGTCACTGCCCGGGGACAATCCCGTCGTCCTTc                                                      <  2:6370356/145‑1 (MQ=255)
                                                                                                     ggggggCCTGGGGGGGGTTCTGCGGGGGTTAACCCGCCCCTCCACCGCCCGCGTCTGGTTGGCGTGGTTGGTGCGGGTTCCCCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGAc                                                 <  2:6083461/150‑1 (MQ=255)
                                                                                                           cgggggCGGGCTTCTGCGGGGTTGAACCCGCCCCTCCAACGCCCGAGTCCGTTTTGCTTCGTCGGTGGGGTTACACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGc                                           <  1:7850791/145‑1 (MQ=255)
                                                                                                             tgtgCGGGCTTCTTCGCGGGTGAACCCGCTCCTACCACGACCCAGTCAGGTTTGCGTCGGAGGAGCGGGTTCACCCGCGGAGAGGCCCCCCCAGGCACCCAAAATTTCACAGCCCCGGGACAAACCCGGCGTCCGGCGTGACATCCGccc                                         >  1:1448303/1‑149 (MQ=25)
                                                                                                             tgtgCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACCACCGAGTCAGGGTTGCGTCGTAGGAGCGGGGTCACCCCCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATACGCCa                                         >  2:1587516/1‑150 (MQ=38)
                                                                                                                  gggCTTCTTCGCGGGTGAAACCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCGTCGTGACATCCGCCATGAcc                                    >  2:3076610/1‑150 (MQ=255)
                                                                                                                      ttcttcGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCCCCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGc                                >  2:4076740/1‑150 (MQ=255)
                                                                                                                        cttcGCGGGTGAACCCGCCCCACCACCGCCGGTGTCTGGTTTGCGTCGTAGGCGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTg                              <  1:8028045/150‑1 (MQ=38)
                                                                                                                         ttcGCGGGTGAACCCGCTCCTACGACGCAAAACTGACTAGGTCGTTGTAGGAGCGGGTTCAACCGCGAAGAAGCCCCCAAAGGCACCCAAAATGTAACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTgg                             >  2:1925642/1‑150 (MQ=14)
                                                                                                                              gggTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCgaga                        >  1:5104865/1‑150 (MQ=255)
                                                                                                                                gTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAAt                      >  2:6898405/1‑150 (MQ=255)
                                                                                                                                   aaCCCGCTCCTACCACGCCCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCg                   <  2:4186376/150‑1 (MQ=255)
                                                                                                                                    aCCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGc                  >  2:6737139/1‑150 (MQ=255)
                                                                                                                                     cccGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCa                 >  2:7644200/1‑150 (MQ=255)
                                                                                                                                     cccGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCa                 >  2:2133048/1‑150 (MQ=255)
                                                                                                                                     cccGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCa                 >  2:8632059/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCaa                >  2:58036/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCaa                >  1:6103801/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCaa                >  2:4639986/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCaa                >  1:6777301/1‑150 (MQ=255)
                                                                                                                                       cGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCaaa               >  1:3859689/1‑150 (MQ=255)
                                                                                                                                       cGCACCAACACCGCCCGCGTCTGGTTGGGGGCGTCGGCGCGGGTCCCCCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCaaa               <  1:509117/150‑1 (MQ=255)
                                                                                                                                        gCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAAc              >  2:6423909/1‑150 (MQ=255)
                                                                                                                                        gCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAAc              >  1:6626560/1‑150 (MQ=255)
                                                                                                                                            cTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACt          >  1:2901742/1‑150 (MQ=255)
                                                                                                                                            cTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACt          >  2:5496289/1‑150 (MQ=255)
                                                                                                                                            cTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACt          >  1:5344455/1‑150 (MQ=255)
                                                                                                                                            cTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACt          >  1:704096/1‑150 (MQ=255)
                                                                                                                                             tACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTg         >  2:8423480/1‑150 (MQ=255)
                                                                                                                                             tACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAAGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTg         >  1:4972867/1‑150 (MQ=255)
                                                                                                                                              aCAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTgg        >  2:6337881/1‑150 (MQ=255)
                                                                                                                                               cAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTGGc       >  2:1089633/1‑150 (MQ=255)
                                                                                                                                                aaCGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTGGCa      >  2:6201115/1‑150 (MQ=255)
                                                                                                                                                aaCGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTGGCa      >  2:1652168/1‑150 (MQ=255)
                                                                                                                                                aaCGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTGGCa      >  2:5130218/1‑150 (MQ=255)
                                                                                                                                                aaCGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTGGCa      >  1:2322989/1‑150 (MQ=255)
                                                                                                                                                 aCGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTGGCaa     >  2:7960548/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACATCGCCCATAGCCATGGGGGCGAGGTGAGCTTGCTGAACCTGCGGGAGGGCGGGCTGCGGGTGACCCTGTATCTGCCGCGGGGTATGGATTGAGTTTTTGGGGTGCCTGTGCGGGCTTCTTCGCGGGTGAACCCGCTCCTACAACGACCGAGTCAGGTTTGCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGCACAGGCACCCAAAATGTCACTGCCCGGTGACAATCCCGTCGTCCTTCGTGACATCCGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTGGCAAGGA  >  NC_002947/1161158‑1161455

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: