Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,712,993 T→A 10.2% intergenic (+34/‑106) adk → / → tsaB adenylate kinase/endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,712,9930TA10.2% 554.3 / 48.0 216intergenic (+34/‑106)adk/tsaBadenylate kinase/endopeptidase
Reads supporting (aligned to +/- strand):  ref base T (87/98);  new base A (14/8);  total (109/107)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-02

Only 100 of 221 total aligned reads displayed.

TGGATTTCTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCTGTTC  >  NC_002947/1712844‑1713142
                                                                                                                                                     |                                                                                                                                                     
tggATTTCTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCt                                                                                                                                                       >  2:4869343/1‑150 (MQ=255)
  gATTTCTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTg                                                                                                                                                     <  2:5988301/150‑1 (MQ=255)
    tttCTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCg                                                                                                                                                   <  1:7072926/150‑1 (MQ=255)
    tttCTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCATGCg                                                                                                                                                   >  1:8662707/1‑150 (MQ=255)
     ttCTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCgg                                                                                                                                                  <  2:3921552/150‑1 (MQ=255)
        tACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGcc                                                                                                                                               >  1:917123/1‑150 (MQ=255)
        tACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGcc                                                                                                                                               <  1:3877040/150‑1 (MQ=255)
         aCCAGAAGCTGTCGGCTGCCCATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCAAGCGGGCCg                                                                                                                                              <  2:3015754/150‑1 (MQ=255)
          ccAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCAAGCGGGCCGt                                                                                                                                             <  2:3166649/150‑1 (MQ=255)
           cAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCCGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGtt                                                                                                                                            <  2:6663225/150‑1 (MQ=255)
               aGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCg                                                                                                                                        >  2:6386615/1‑150 (MQ=255)
                  tGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTTTCGttt                                                                                                                                     >  1:8890742/1‑150 (MQ=255)
                  tGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGttt                                                                                                                                     <  1:5819594/150‑1 (MQ=255)
                  tGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGgtt                                                                                                                                     >  1:7174785/1‑150 (MQ=255)
                      ggCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAg                                                                                                                                 >  2:2605872/1‑150 (MQ=255)
                         tGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGcc                                                                                                                              >  2:1694491/1‑150 (MQ=255)
                         tGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCAAGCGGGCCGTTGt                                                                                                                                          >  2:3557240/1‑138 (MQ=255)
                          gCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCg                                                                                                                             <  2:5171714/150‑1 (MQ=255)
                           ccAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGc                                                                                                                            >  2:1434746/1‑150 (MQ=255)
                                ggcggcAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAgc                                                                                                                       <  1:307500/150‑1 (MQ=255)
                                  cggcAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAgcgc                                                                                                                     >  2:187504/1‑150 (MQ=255)
                                   ggcAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGcc                                                                                                                    >  2:2490769/1‑150 (MQ=255)
                                   ggcAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGcc                                                                                                                    <  1:7758201/150‑1 (MQ=255)
                                    gcAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGccc                                                                                                                   >  2:8432168/1‑150 (MQ=255)
                                     cAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGcccc                                                                                                                  <  2:8864423/150‑1 (MQ=255)
                                      aaGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGccccc                                                                                                                 >  1:1210470/1‑150 (MQ=255)
                                        gCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTc                                                                                                               <  2:7131495/150‑1 (MQ=255)
                                        gCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTc                                                                                                               <  2:8722533/150‑1 (MQ=255)
                                         ccGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCa                                                                                                              <  2:2063185/150‑1 (MQ=255)
                                          cGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTTAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAg                                                                                                             <  1:4620587/150‑1 (MQ=255)
                                            aaGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAATCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGa                                                                                                           <  1:1334372/150‑1 (MQ=255)
                                               tACAGCCACCTCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGacac                                                                                                        <  2:3691953/150‑1 (MQ=255)
                                               tACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGacac                                                                                                        <  2:1922070/150‑1 (MQ=255)
                                                aCAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGAcaca                                                                                                       >  1:7408731/1‑150 (MQ=255)
                                                aCAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCACAGGAcaca                                                                                                       >  1:8513921/1‑150 (MQ=255)
                                                         tCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCgctg                                                                                              <  2:5304264/150‑1 (MQ=255)
                                                            aaGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCgctgctg                                                                                           <  1:1489008/150‑1 (MQ=255)
                                                              ggTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCgctgctgct                                                                                         <  2:8845066/150‑1 (MQ=255)
                                                                  tCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCCCAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCtat                                                                                     >  1:8043950/1‑150 (MQ=255)
                                                                     gCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACt                                                                                  >  2:8381147/1‑150 (MQ=255)
                                                                     gCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACt                                                                                  >  2:7369049/1‑150 (MQ=255)
                                                                      cTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGGTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTg                                                                                 >  1:2057399/1‑150 (MQ=255)
                                                                         ggTGGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGGTAGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCGCTATACTGCCg                                                                              >  2:9045311/1‑150 (MQ=255)
                                                                           ttGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCGTGCGGGCCGTTGTAGTTTGTAGGCCGCACAGCGCCCCCGCAGGTAACAACGCCGCTGCTGCTCTATAATGCCGct                                                                            >  1:8840639/1‑150 (MQ=255)
                                                                            tGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGctc                                                                           <  2:7587018/150‑1 (MQ=255)
                                                                            tGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGctc                                                                           <  1:162895/150‑1 (MQ=255)
                                                                             gAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCGTGCGGGCCGTTGTCGTTTGTTGGCCGCACAGCGCCCCCGCAGGACCCAACCCCGCGGCTGCGCGATAAGGCCGcggt                                                                          >  1:3131295/1‑147 (MQ=255)
                                                                               aGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCttt                                                                        <  2:108420/150‑1 (MQ=255)
                                                                               aGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCttt                                                                        <  2:4442199/150‑1 (MQ=255)
                                                                               aGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCttt                                                                        <  2:7483694/150‑1 (MQ=255)
                                                                                gCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCtttt                                                                       <  2:3484719/150‑1 (MQ=255)
                                                                                   aTCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCa                                                                    <  2:900439/150‑1 (MQ=255)
                                                                                   aTCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCa                                                                    <  2:7793214/150‑1 (MQ=255)
                                                                                    tCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGc                                                                             >  2:2200430/1‑140 (MQ=255)
                                                                                    tCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCAt                                                                   >  2:7962469/1‑150 (MQ=255)
                                                                                    tCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGGTTGTAGGCCGCACAGCGCCCCCTCAGGACACCACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCAt                                                                   >  2:43352/1‑150 (MQ=255)
                                                                                      aCCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCCCCCCCGCACGACCCAACCCCCCAGCTGCTCTATACAGCCGCTCATTTTGAtta                                                                 >  1:7643907/1‑149 (MQ=255)
                                                                                       ccGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCg                                                                >  2:858196/1‑150 (MQ=255)
                                                                                                 ccTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATa                                                      <  2:4690211/150‑1 (MQ=255)
                                                                                                  cTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCAGGATAc                                                     >  1:8261279/1‑150 (MQ=255)
                                                                                                  cTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGGTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCAGGATAc                                                     >  1:8261217/1‑150 (MQ=255)
                                                                                                   tGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATAcc                                                    <  2:8026475/150‑1 (MQ=255)
                                                                                                    ggCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATAccc                                                   <  2:7214483/150‑1 (MQ=255)
                                                                                                     gCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCCCAGCGCCCCCACAGGACACAACCCCGCGGCTGCTCTATACAGCCGCTCTTTTTCATTAGCACCGGGGTACCCg                                                  >  1:1284648/1‑150 (MQ=255)
                                                                                                     gCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGGTGTCGGTTGGAGGCCGCCCCGCGCCCCCCCAGGACCCCACGCCGCGGCGGCTCGATACGGCCGCACTTTTTCATTCGCCCCAGGCGAcccc                                                  >  1:7342744/1‑149 (MQ=255)
                                                                                                         cccTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCg                                              >  2:6919237/1‑150 (MQ=255)
                                                                                                         accTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTTTCGTTTGTAGGCCGCACAGCCCCCCCAAACGACACAACGCCCCTGCTGCTCTATACAGACGCACATTTTAATTCGCACCAGGATACCCgatag                                              >  2:2299648/2‑146 (MQ=255)
                                                                                                           cTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGAt                                            <  1:5251089/150‑1 (MQ=255)
                                                                                                           cTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGAt                                            <  1:5185175/150‑1 (MQ=255)
                                                                                                           cTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGGTTGTATGCCGCACAGCGCCCCCTCAGGACCCAAACCCGCGGCTGGTCTATACTGCCGCTCGTTTTCAGTAGCCCCAGGATAccccgtcgct                                            >  1:7119109/1‑143 (MQ=255)
                                                                                                               gCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGacc                                        <  1:1054630/150‑1 (MQ=255)
                                                                                                                  gATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGaccacc                                     >  1:664198/1‑150 (MQ=255)
                                                                                                                   aTCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCGCCTGGATACCCGTTCGATGACCAccc                                    >  1:2895085/1‑150 (MQ=255)
                                                                                                                   aTCCGGCGCCTGACGCACCTCGACAACGACCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCt                                                           <  1:5425346/127‑1 (MQ=255)
                                                                                                                     ccGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCctg                                  >  1:6023873/1‑150 (MQ=255)
                                                                                                                       ggCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCctgct                                >  2:6481166/1‑150 (MQ=255)
                                                                                                                        gcgcCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTAGACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCctgctg                               <  2:5780536/150‑1 (MQ=255)
                                                                                                                          gcCTGACGCACCTCGACAACGGCCCGCGTGCGGGCCGGTGTCGGTTGTAGGCCGCACAGCGCCCCCCCAGGACCCCACCCCCCGGCTGCGCGATACAGCCGCGCTTTTTCCATCGCACCGGGATACCCCGGCGATGACCCCCCAGCGGGc                             >  1:7195863/1‑150 (MQ=255)
                                                                                                                             tGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCt                          <  1:5029316/150‑1 (MQ=255)
                                                                                                                              gACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCAGCGGGcccag                         >  2:7873705/1‑148 (MQ=255)
                                                                                                                               aCGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTgg                        <  1:4873953/150‑1 (MQ=255)
                                                                                                                                cGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGa                       <  2:4789634/150‑1 (MQ=255)
                                                                                                                                      tCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGc                 <  1:8815957/150‑1 (MQ=255)
                                                                                                                                        gACAACGGCCCGCAAGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGCTACCGCCa               >  2:8795128/1‑150 (MQ=255)
                                                                                                                                        gACAACGGCCCGCAAGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCGGGATACCGCCa               >  1:3672314/1‑150 (MQ=255)
                                                                                                                                         aCAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCAc              >  1:2118955/1‑150 (MQ=255)
                                                                                                                                          cAACGGCCCGCAAGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCAcc             >  2:421198/1‑150 (MQ=255)
                                                                                                                                           aaCGGCCCGCAAGCGGGCCGTTGGCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCGCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGCTGACCACCCTGCTGGCCCTGGATACCGCCACCg            >  2:6950387/1‑150 (MQ=255)
                                                                                                                                           aaCGGACCGCAAGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAAAGCCGCGGCTGCTCTATACTGCCGCTCTTTTTAATTCGCACCTGGATACCCGTTCGATGTCCACCCTGCTGGCCCGGGGTACCGCCAccc            >  1:1183658/1‑149 (MQ=255)
                                                                                                                                            aCGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGa           <  2:663842/150‑1 (MQ=255)
                                                                                                                                             cGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGaa          <  2:4048213/150‑1 (MQ=255)
                                                                                                                                             cGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGaa          <  2:5459971/150‑1 (MQ=255)
                                                                                                                                              ggCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCAATCGCACCTGGATACCCGTTCGATGACCACCCAGCTGGCCCGGGATACCGCCACCGAAg         >  1:3246797/1‑150 (MQ=255)
                                                                                                                                              ggCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTAATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCGGGATAACGCCACCGAAg         >  1:4285063/1‑150 (MQ=255)
                                                                                                                                                cccGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCACAGGACACAACGCCCCTGCTGCTCTAAACTGCCGCTCATTTTCATTCGCACCTGGATAACCCGTCGATGACCACCCTGCTGGCCCTGGATACCGCCCCCGAAGcc       >  1:178644/1‑150 (MQ=255)
                                                                                                                                                 ccGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCt      >  2:8323748/1‑150 (MQ=255)
                                                                                                                                                 ccGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCt      >  1:5962618/1‑150 (MQ=255)
                                                                                                                                                     ttGCGGGCCGTTGTCGTTTGTAGGCCGCATAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCGGGATACCGCCACCCAAGCCGGTTc  >  1:2691531/1‑150 (MQ=255)
                                                                                                                                                     ttGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTAGCACCTGGATACCCGTTAGATGTCCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCTGTTc  >  1:653556/1‑150 (MQ=255)
                                                                                                                                                     ttGCGGGCCGTTGTCGTTTGTAGGCCGCAAAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCAGCTGGCCCTGGATACCGCCACCCAAGCCTGTTc  >  1:6463838/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGATTTCTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCTGTTC  >  NC_002947/1712844‑1713142

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: