Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,160,962 T→A 5.8% intergenic (+321/+547) PP_5739 → / ← PP_5402 hypothetical protein/NAD‑dependent acetyl‑CoA synthetase deacetylase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,160,9620TA 5.8% 1037.4 / 36.5 308intergenic (+321/+547)PP_5739/PP_5402hypothetical protein/NAD‑dependent acetyl‑CoA synthetase deacetylase subunit
Reads supporting (aligned to +/- strand):  ref base T (143/147);  new base A (12/6);  total (155/153)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

Only 100 of 321 total aligned reads displayed.

CAAGCCGCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCACGCCAGTCGTCCGG  >  NC_002947/6160814‑6161111
                                                                                                                                                    |                                                                                                                                                      
  aGCCGCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCaa                                                                                                                                                     <  1:5783995/150‑1 (MQ=255)
   gCCGCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAg                                                                                                                                                    <  1:2032850/150‑1 (MQ=255)
   gCCGCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAg                                                                                                                                                    <  1:9069361/150‑1 (MQ=255)
          gTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCTTGAGAGGGCGGCtt                                                                                                                                            <  1:496934/150‑1 (MQ=255)
           tGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCTTGAGAGGGCGGCttt                                                                                                                                           <  1:3423908/150‑1 (MQ=255)
           tGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCttt                                                                                                                                           <  2:3767126/150‑1 (MQ=255)
            gATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAATCCGCCCTCTCAAGGCGGCtttt                                                                                                                                          >  2:5929409/1‑150 (MQ=255)
             aTCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCttttt                                                                                                                                         >  1:2717412/1‑150 (MQ=255)
                cATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCtttttttt                                                                                                                                      <  1:4752275/150‑1 (MQ=255)
                cATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCtttttttt                                                                                                                                      <  2:7515340/150‑1 (MQ=255)
                cATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCtttttttt                                                                                                                                      <  2:1236052/150‑1 (MQ=255)
                cATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCtttttttg                                                                                                                                      <  2:6372003/150‑2 (MQ=255)
                  tttCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGc                                                                                                                                    >  1:1611471/1‑150 (MQ=255)
                     cGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGTTTTAGCTATAGTCCTCATCGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTATAATGCCGCCCTCTCAAGGCGGGTTTTTTTTGCCTc                                                                                                                                 >  2:4431297/1‑150 (MQ=255)
                     cGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGccaa                                                                                                                                 <  1:6396682/150‑3 (MQ=255)
                     cGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTc                                                                                                                                 <  1:3218568/150‑1 (MQ=255)
                      gCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCg                                                                                                                                >  2:7862525/1‑150 (MQ=255)
                         ttAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGAt                                                                                                                             <  2:2800788/150‑1 (MQ=255)
                           aGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTc                                                                                                                           <  1:6873898/150‑1 (MQ=255)
                               gCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTAt                                                                                                                       <  1:4610099/150‑1 (MQ=255)
                               gCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTAt                                                                                                                       <  1:4195531/150‑1 (MQ=255)
                                cTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATc                                                                                                                      <  2:3448416/150‑1 (MQ=255)
                                 tAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCg                                                                                                                     >  1:7051043/1‑150 (MQ=255)
                                 tAAATTTGTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCg                                                                                                                     <  2:4806297/150‑1 (MQ=255)
                                  aaaTTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCgg                                                                                                                    >  1:5983502/1‑150 (MQ=255)
                                           gTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACcgtcgt                                                                                                           >  2:3993060/1‑150 (MQ=255)
                                            tttGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGtg                                                                                                          >  1:5565541/1‑150 (MQ=255)
                                                gTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGtgctgc                                                                                                      >  1:4254000/1‑150 (MQ=255)
                                                   gCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAg                                                                                                   <  1:5915060/150‑1 (MQ=255)
                                                    cGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAgt                                                                                                  <  2:9156582/150‑1 (MQ=255)
                                                      tCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAgtgt                                                                                                >  1:6319551/1‑150 (MQ=255)
                                                        attaattaTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCg                                                                                              >  2:3083811/1‑150 (MQ=255)
                                                        attaattaTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCg                                                                                              >  2:8661231/1‑150 (MQ=255)
                                                         ttaattaTGTTTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGGCGACGAACCTCTCATTTCGTTGGCTATTTTGCCGCGCTCTCATGGGGGGTTTTTTTTGGCTGGGTTTGGTTTCGGTAGCGTGGTGCTGGGAgtgtggg                                                                                             >  2:8794061/1‑147 (MQ=255)
                                                         ttaattaTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGc                                                                                             <  1:994925/150‑1 (MQ=255)
                                                            attaTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGa                                                                                          <  1:7280673/150‑1 (MQ=255)
                                                             ttaTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGaa                                                                                         >  2:5280439/1‑150 (MQ=255)
                                                             ttaTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGaa                                                                                         <  2:1984717/150‑1 (MQ=255)
                                                             ttaTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGaa                                                                                         >  2:2520327/1‑150 (MQ=255)
                                                             ttaTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGaa                                                                                         >  2:7137949/1‑150 (MQ=255)
                                                             ttaTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGaa                                                                                         >  1:5058361/1‑150 (MQ=255)
                                                              taTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAAc                                                                                        >  1:232357/1‑150 (MQ=255)
                                                              taTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGACAGTGTCGCCGAAc                                                                                        >  2:8487431/1‑150 (MQ=255)
                                                                   tttAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGAt                                                                                   <  1:3766863/150‑1 (MQ=255)
                                                                    ttAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGAGTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATc                                                                                  >  2:6905528/1‑150 (MQ=255)
                                                                    ttAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCATCTCAAGGCGGCTTTTTTTTGCCTAGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATc                                                                                  >  1:5700223/1‑150 (MQ=255)
                                                                         tataGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACtgt                                                                             >  2:7568724/1‑150 (MQ=255)
                                                                          ataGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACtgtg                                                                            <  2:7387984/150‑1 (MQ=255)
                                                                           taGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGa                                                                           <  2:210505/150‑1 (MQ=255)
                                                                           taGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGa                                                                           <  1:2851441/150‑1 (MQ=255)
                                                                              tCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGa                                                                        <  2:8337374/150‑1 (MQ=255)
                                                                                cTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGAcc                                                                      >  1:8115082/1‑150 (MQ=255)
                                                                                cTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGAcc                                                                      >  1:5352961/1‑150 (MQ=255)
                                                                                cTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGAcc                                                                      >  2:297621/1‑150 (MQ=255)
                                                                                cTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGAcc                                                                      >  2:8081800/1‑150 (MQ=255)
                                                                                 tCAACGTTCGAGTGTTTCAGTACGTCGCCGACGAACCTCTCATTTCGTTGGGTAAAATGCCGCCCTCTCTAGGCGGGTTTTTTTTGGCTCGGATTCGTTTCGGTACCGTCGTGCTGGGAGTGTCGCGGGAGTTGATGTCTGTGTGGAcct                                                                     >  2:2262702/1‑149 (MQ=255)
                                                                                  cAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAg                                                                    <  1:4008319/150‑1 (MQ=255)
                                                                                  cAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAg                                                                    >  1:8072093/1‑150 (MQ=255)
                                                                                   aaCGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGa                                                                   <  2:4695778/150‑1 (MQ=255)
                                                                                   aaCGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGa                                                                   <  2:7273832/150‑1 (MQ=255)
                                                                                   aaCGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTATTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGa                                                                   <  2:486147/150‑1 (MQ=255)
                                                                                       ttCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTAt                                                               <  1:3783159/150‑1 (MQ=255)
                                                                                       ttCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTAt                                                               >  2:6749648/1‑150 (MQ=255)
                                                                                        tCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATc                                                              <  2:5395075/150‑1 (MQ=255)
                                                                                         cGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCa                                                             <  1:4524923/150‑1 (MQ=255)
                                                                                                 tCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTaca                                                     <  1:8521876/150‑1 (MQ=255)
                                                                                                  cAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTacac                                                    <  1:1743530/150‑1 (MQ=255)
                                                                                                       cGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAac                                               <  2:4588056/150‑1 (MQ=255)
                                                                                                         tCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAacca                                             >  1:8467623/1‑150 (MQ=255)
                                                                                                          cGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAaccac                                            <  1:564261/150‑1 (MQ=255)
                                                                                                             cgacgaACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTAGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCtg                                         <  1:1009017/150‑1 (MQ=255)
                                                                                                                cgaACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGa                                      <  2:1537158/150‑1 (MQ=255)
                                                                                                                cgaACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGa                                      <  1:4132179/150‑1 (MQ=255)
                                                                                                                 gaACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGAc                                     >  2:6100370/1‑150 (MQ=255)
                                                                                                                 gaACCTCTCAAGTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGAc                                     <  1:6807306/150‑1 (MQ=255)
                                                                                                                   aCCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTg                                   >  2:3326/1‑150 (MQ=255)
                                                                                                                   aCCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTg                                   >  1:3661524/1‑150 (MQ=255)
                                                                                                                      tctcAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTc                                >  1:7760066/1‑150 (MQ=255)
                                                                                                                       ctcAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCg                               >  1:2209159/1‑150 (MQ=255)
                                                                                                                        tcAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCgg                              >  1:1969166/1‑150 (MQ=255)
                                                                                                                         cAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGc                             <  1:68647/150‑1 (MQ=255)
                                                                                                                          aaTTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCGCGGATTCGGATCCGTACCGTCGTGCTGCCAGTGGCGCCGAACATTATCACAGTGGGGACCAGATTATCATTCGTACACACAAACACCCGTGACATTTCGGcc                            >  2:8169349/1‑150 (MQ=255)
                                                                                                                               gTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTg                       <  2:7226842/150‑1 (MQ=255)
                                                                                                                               gTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTg                       <  1:7862525/150‑1 (MQ=255)
                                                                                                                                 tGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGc                     <  1:7421507/150‑1 (MQ=255)
                                                                                                                                   aaaaaaaaGCCGCCTTGAGAGGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCaa                   >  1:7477035/4‑150 (MQ=255)
                                                                                                                                    cTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAAc                  >  2:6576928/1‑150 (MQ=255)
                                                                                                                                     tAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAAcc                 >  2:1722724/1‑150 (MQ=255)
                                                                                                                                      aaaaaGCCGCCTTGAGAGGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCtgt                                        >  2:1378795/1‑126 (MQ=255)
                                                                                                                                      aaaaaGCCGCCTTGAGAGGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCtgt                                        <  1:1378795/126‑1 (MQ=255)
                                                                                                                                      aaaaaGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCa                <  2:6975167/150‑1 (MQ=255)
                                                                                                                                      aaaaaGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCa                <  2:6674750/150‑1 (MQ=255)
                                                                                                                                        aaaGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCACg              >  1:3930918/1‑150 (MQ=255)
                                                                                                                                          aGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCACGcc            >  1:6136577/1‑150 (MQ=255)
                                                                                                                                           gccgccCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCACGCCa           <  1:1504518/150‑1 (MQ=255)
                                                                                                                                                 ctctCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGCACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCAGTGACAGTTCGGCCTGTTGGCAACCACGCCAgtcgtc     >  1:8724250/1‑150 (MQ=255)
                                                                                                                                                  tctcAAGGCGGCTTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCACGCCAgtcgtc     >  1:5320073/1‑150 (MQ=255)
                                                                                                                                                  tctcAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCACGCCAGTCGTcc    >  1:8904258/1‑150 (MQ=255)
                                                                                                                                                  tctcAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGGCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCACGCCAGTCGTcc    <  1:3053847/150‑1 (MQ=255)
                                                                                                                                                    tcAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCACGCCAGTCGTCCgg  >  1:1250909/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CAAGCCGCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCACGCCAGTCGTCCGG  >  NC_002947/6160814‑6161111

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