Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 546,010 G→T 6.3% intergenic (+47/‑26) rpoC → / → PP_mr08 DNA‑directed RNA polymerase subunit beta'/rpsL_psuedo ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947546,0100GT 6.3% 906.6 / 35.5 286intergenic (+47/‑26)rpoC/PP_mr08DNA‑directed RNA polymerase subunit beta'/rpsL_psuedo ncRNA
Reads supporting (aligned to +/- strand):  ref base G (103/164);  new base T (6/13);  total (109/177)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 287 total aligned reads displayed.

GCGAGCGCAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGACAAT  >  NC_002947/545863‑546159
                                                                                                                                                   |                                                                                                                                                     
gcgagcgCAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGaa                                                                                                                                                     <  1:1115261/150‑1 (MQ=255)
  gagcgCAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAAt                                                                                                                                                   >  2:4382073/1‑150 (MQ=255)
    gcgcAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGt                                                                                                                                                 <  1:4183881/150‑1 (MQ=255)
     cgcAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGtt                                                                                                                                                <  1:3546679/150‑1 (MQ=255)
     cgcAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGtt                                                                                                                                                <  1:1936816/150‑1 (MQ=255)
      gcAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTg                                                                                                                                               >  2:2221430/1‑150 (MQ=255)
      gcAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTg                                                                                                                                               <  1:2398791/150‑1 (MQ=255)
       cAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGc                                                                                                                                              <  2:1823544/150‑1 (MQ=255)
       cAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGc                                                                                                                                              <  1:3230768/150‑1 (MQ=255)
       cAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGc                                                                                                                                              >  1:2572860/1‑150 (MQ=255)
       cAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGc                                                                                                                                              <  1:1338777/150‑1 (MQ=255)
       cAAGCGTCGCCGTGATGCCCACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCCGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCAGCGGCGACAGGAAATGTTGc                                                                                                                                              >  2:195636/1‑150 (MQ=255)
        aaGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACCTGAAATGTTGcc                                                                                                                                             >  1:3273821/1‑150 (MQ=255)
        aaGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGcc                                                                                                                                             <  2:4003558/150‑1 (MQ=255)
        aaGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGcc                                                                                                                                             <  2:311286/150‑1 (MQ=255)
          gCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACTGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCAACATTTCATGTCGCCga                                                                                                                                           <  2:4259699/150‑1 (MQ=255)
          gCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCTATCATCGGCGACATGAAATGTTGCCga                                                                                                                                           <  2:5904006/150‑1 (MQ=255)
          gCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCAACATTTCATGTCGCCga                                                                                                                                           <  2:2108771/150‑1 (MQ=255)
           cGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCgat                                                                                                                                          <  1:2287409/150‑1 (MQ=255)
           cGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCgat                                                                                                                                          <  2:2270270/150‑1 (MQ=255)
            gTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCgatg                                                                                                                                         >  1:4915391/1‑150 (MQ=255)
               gCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATc                                                                                                                                      <  1:1610820/150‑1 (MQ=255)
                 cGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATcgg                                                                                                                                    <  1:5706866/150‑1 (MQ=255)
                     aTGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGCCCGATCATCGGCAACATTTCATGTCGCCTATGATcgcccga                                                                                                                                <  1:710211/150‑6 (MQ=37)
                        ccGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGa                                                                                                                             <  1:4808446/150‑1 (MQ=255)
                         cGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAg                                                                                                                            >  2:7219607/1‑150 (MQ=255)
                          gACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATcgggc                                                                                                                                  >  1:445915/1‑143 (MQ=255)
                               aCCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCgg                                                                                                                      <  2:831783/150‑1 (MQ=255)
                               aCCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCgg                                                                                                                      >  2:539580/1‑150 (MQ=255)
                                ccGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCggg                                                                                                                     <  2:2540660/150‑1 (MQ=255)
                                  gCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGc                                                                                                                   >  2:4821965/1‑150 (MQ=255)
                                     gCGTGAACGCGCCAGTGAAGTGGAAGCTGCACTGACCGACGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAACTGTGGCCGATGATCGGGCGGGGAGGGCCGGGGCCTc                                                                                                                <  1:3863436/150‑1 (MQ=255)
                                            aGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGa                                                                                                         <  1:2845548/150‑1 (MQ=255)
                                               gcCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGTTAATTAAGTACAGGGCATGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGTTGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTg                                                                                                      <  2:6318873/150‑1 (MQ=255)
                                                 cAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTggg                                                                                                    >  2:7278385/1‑150 (MQ=255)
                                                      aaGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCaa                                                                                               <  2:4625831/150‑1 (MQ=255)
                                                      aaGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCaa                                                                                               <  2:674591/150‑1 (MQ=255)
                                                        ggggAAGCTGCCCTGACCGAAGCGCTGAATTCCCGCGGTAATTAAGTACAGGGCAAGGCCCCGGCACGCCGCGGGCGATCATCGGCGACATGAACTGTGGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGa                                                                                             <  1:3729104/148‑1 (MQ=255)
                                                        gTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGa                                                                                             >  2:339261/1‑150 (MQ=255)
                                                          ggAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGTTATTTAAGTACCGGGCACGGCCCGGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATc                                                                                           <  1:4363424/150‑1 (MQ=255)
                                                          ggAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATc                                                                                           <  2:5555943/150‑1 (MQ=255)
                                                          ggAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATc                                                                                           >  2:6977703/1‑150 (MQ=255)
                                                           gAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATcc                                                                                          >  2:3215442/1‑150 (MQ=255)
                                                             aGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCctc                                                                                        <  2:5805190/150‑1 (MQ=255)
                                                                 gCACTGACCGAAGCGCTTAATTCCCGCGGTAATTAAGTCCAGGGCAAGGCCCCGGCAAGCCTCGGGCTATAATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTa                                                                                    <  1:5023266/150‑1 (MQ=255)
                                                                 gCACTGACCGAAGCGCTTAATACCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGAAATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTa                                                                                    <  1:6292849/150‑1 (MQ=255)
                                                                     tGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACt                                                                                >  2:1460395/1‑150 (MQ=255)
                                                                        ccGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACtttt                                                                             >  2:3623589/1‑150 (MQ=255)
                                                                        ccGAAGCGATGAATTCCAGCGGTAATTAAGTACTGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACtttt                                                                             <  1:4709490/150‑1 (MQ=255)
                                                                           aaGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTa                                                                          <  2:4220888/150‑1 (MQ=255)
                                                                                 tGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACtt                                                                               <  2:2880800/139‑1 (MQ=255)
                                                                                 tGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTa                                                                    >  2:115359/1‑150 (MQ=255)
                                                                                 tGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTa                                                                    >  2:2418706/1‑150 (MQ=255)
                                                                                 tGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCATa                                                                    >  1:1247953/1‑150 (MQ=255)
                                                                                   aaTTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTaaa                                                                  >  1:6531571/1‑150 (MQ=255)
                                                                                       ccAGCGGTAATTAAGTACCGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTg                                                              <  2:2622843/150‑1 (MQ=255)
                                                                                       ccAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTg                                                              <  2:6559641/150‑1 (MQ=255)
                                                                                        ctgCGGTAATTCAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTgg                                                             <  2:7152807/148‑1 (MQ=255)
                                                                                        cAGCGTTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTgg                                                             <  1:3441772/150‑1 (MQ=255)
                                                                                        cAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTgg                                                             <  2:858919/150‑1 (MQ=255)
                                                                                          gCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTg                                                           >  1:6952696/1‑150 (MQ=255)
                                                                                           cGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGa                                                          <  1:6709582/150‑1 (MQ=255)
                                                                                           cGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGa                                                          <  2:4299065/150‑1 (MQ=255)
                                                                                           cGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGa                                                          <  2:1221340/150‑1 (MQ=255)
                                                                                             gTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTTCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAgg                                                        >  1:5786444/1‑150 (MQ=255)
                                                                                             gTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAgg                                                        <  2:6171262/150‑1 (MQ=255)
                                                                                                      tACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGCCATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACATTTGTACCCTTAAATTTGGTGAGGCTCCGtctc                                               >  1:3095588/1‑150 (MQ=255)
                                                                                                       aCAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCg                                              <  2:5852740/150‑1 (MQ=255)
                                                                                                        cAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGc                                             >  1:5585145/1‑150 (MQ=255)
                                                                                                            gCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGt                                                            >  2:6016682/1‑131 (MQ=255)
                                                                                                            gCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGt                                                            <  1:6016682/131‑1 (MQ=255)
                                                                                                                ggCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAAttttt                                     <  1:7185773/150‑1 (MQ=255)
                                                                                                                ggCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATACTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAAttttt                                     >  1:2181644/1‑150 (MQ=255)
                                                                                                                ggCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTGGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAAttttt                                     <  2:2819359/150‑1 (MQ=255)
                                                                                                                  ccccGGCACGCCTCGGGCGATCTTCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTgg                                   <  1:806722/150‑1 (MQ=255)
                                                                                                                  ccccGGCAAGCCTCGGGCGATCATCGGCGACCGGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGGGCAAGATCCGCGTGAGACTTTTGTACCCTTAAATTTGGTGAGGCGCCGTCTCGCCAAGTGTTgg                                   >  1:6257477/1‑150 (MQ=255)
                                                                                                                  ccccGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTgt                                   <  1:6456668/150‑2 (MQ=255)
                                                                                                                      ggCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCttt                               >  1:3882026/1‑150 (MQ=255)
                                                                                                                         acgCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTGTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCtt                            <  2:3795909/148‑1 (MQ=255)
                                                                                                                          aGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTg                           <  2:4986222/150‑1 (MQ=255)
                                                                                                                           gCCTCGGGCGATCATCGGCGACATTAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAGTTTTTGGCTTTCTTGc                          >  2:1312962/1‑150 (MQ=255)
                                                                                                                             cTCGGGCGATCATCGGCGACATGACCTGGTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCaa                        <  2:1730012/150‑1 (MQ=255)
                                                                                                                             cTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCaa                        <  1:164319/150‑1 (MQ=255)
                                                                                                                             cTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTATTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCaa                        <  2:3075150/150‑1 (MQ=255)
                                                                                                                                gggCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGAc                     <  2:3703969/150‑1 (MQ=255)
                                                                                                                                 ggCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACa                    >  2:2695838/1‑150 (MQ=255)
                                                                                                                                   cGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAAt                  <  1:5644719/150‑1 (MQ=255)
                                                                                                                                    gATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATa                 <  1:539580/150‑1 (MQ=255)
                                                                                                                                    gATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATa                 >  1:4637369/1‑150 (MQ=255)
                                                                                                                                    gATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATa                 <  2:3875600/150‑1 (MQ=255)
                                                                                                                                    gATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATa                 <  2:5247182/150‑1 (MQ=255)
                                                                                                                                     atcatgGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTCTAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATag                <  1:1835986/150‑1 (MQ=255)
                                                                                                                                      tcatcGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCtt                            >  2:1602297/1‑137 (MQ=255)
                                                                                                                                      tcatcGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATaga               >  2:6693811/1‑150 (MQ=255)
                                                                                                                                      tcatcGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATaga               <  2:1349199/150‑1 (MQ=255)
                                                                                                                                         tcGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTc            <  2:7108081/150‑1 (MQ=255)
                                                                                                                                         tcGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTc            <  2:117851/150‑1 (MQ=255)
                                                                                                                                         tcGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTc            <  1:3646309/150‑1 (MQ=255)
                                                                                                                                           ggCGACATTAACTGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGc          <  2:1855461/150‑1 (MQ=255)
                                                                                                                                               aCATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGTTGAAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGATCAGTCTCGCCAATTTTTGTCTTTCTTGAAAGACAATATATTCTCAAGa      >  1:5878918/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCGAGCGCAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGACAAT  >  NC_002947/545863‑546159

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: