Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 546,017 T→C 6.0% intergenic (+54/‑19) 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,0170TC 6.0% 900.3 / 27.9 286intergenic (+54/‑19)rpoC/PP_mr08DNA‑directed RNA polymerase subunit beta'/rpsL_psuedo ncRNA
Reads supporting (aligned to +/- strand):  ref base T (103/160);  new base C (5/12);  total (108/178)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

Only 100 of 297 total aligned reads displayed.

CGCAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGACAATCAGTGGA  >  NC_002947/545868‑546166
                                                                                                                                                     |                                                                                                                                                     
cgcAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGtt                                                                                                                                                       <  1:6438118/150‑1 (MQ=255)
 gcAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTg                                                                                                                                                      >  2:2221430/1‑150 (MQ=255)
 gcAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTg                                                                                                                                                      >  1:673542/1‑150 (MQ=255)
 gcAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTg                                                                                                                                                      <  1:2398791/150‑1 (MQ=255)
  cAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGc                                                                                                                                                     <  1:3323211/150‑1 (MQ=255)
  cAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGc                                                                                                                                                     <  1:3230763/150‑1 (MQ=255)
  cAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGc                                                                                                                                                     >  1:2572860/1‑150 (MQ=255)
  cAAGCGTCGCCGTGATGCCCACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCCGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCAGCGGCGACAGGAAATGTTGc                                                                                                                                                     >  2:195636/1‑150 (MQ=255)
   aaGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGcc                                                                                                                                                    <  2:4003558/150‑1 (MQ=255)
   aaGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGcc                                                                                                                                                    >  1:3273780/1‑150 (MQ=255)
   aaGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGcc                                                                                                                                                    <  2:311286/150‑1 (MQ=255)
   aaGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAAGGTTGcc                                                                                                                                                    >  1:1308705/1‑150 (MQ=255)
   aaGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACAGGAAATGTTGcc                                                                                                                                                    >  2:5994107/1‑150 (MQ=255)
    aGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCAACATTTCATGTCGCCg                                                                                                                                                   <  2:1169215/150‑1 (MQ=255)
     gCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACTGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCAACATTTCATGTCGCCga                                                                                                                                                  <  2:4259699/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)
       gTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCAACATTTCATGTCGCCgatg                                                                                                                                                <  2:1666878/150‑1 (MQ=255)
        tCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCgatga                                                                                                                                               >  1:5845820/1‑150 (MQ=255)
            cGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATcgg                                                                                                                                           <  2:2591383/150‑1 (MQ=255)
               gATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATcgggcg                                                                                                                                        <  1:4060898/150‑1 (MQ=255)
                  gCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCgggg                                                                                                                                     >  2:5016712/1‑150 (MQ=255)
                   ccGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGa                                                                                                                                    <  1:4808446/150‑1 (MQ=255)
                    cGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAg                                                                                                                                   >  2:7219607/1‑150 (MQ=255)
                     gACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAgg                                                                                                                                  >  1:7108081/1‑150 (MQ=255)
                     gACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAgg                                                                                                                                  >  2:3756997/1‑150 (MQ=255)
                          aCCGCTGCGTGTAAGCGCCAGTGACGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCACGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCgg                                                                                                                             <  2:6227488/150‑1 (MQ=255)
                          aCCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCgg                                                                                                                             <  2:3463345/150‑1 (MQ=255)
                          aCCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCgg                                                                                                                             <  2:831783/150‑1 (MQ=255)
                          aCCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCgg                                                                                                                             >  1:2002956/1‑150 (MQ=255)
                           ccGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCggg                                                                                                                            <  2:2540660/150‑1 (MQ=255)
                             gCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGc                                                                                                                          >  2:3679860/1‑150 (MQ=255)
                             gCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGc                                                                                                                          >  2:4821965/1‑150 (MQ=255)
                             gCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGc                                                                                                                          >  2:6709582/1‑150 (MQ=255)
                                gCGTGAACGCGCCAGTGAAGTGGAAGCTGCACTGACCGACGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAACTGTGGCCGATGATCGGGCGGGGAGGGCCGGGGCCTc                                                                                                                       <  1:3863436/150‑1 (MQ=255)
                                           ccAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTgg                                                                                                            <  2:5847441/150‑1 (MQ=255)
                                            cAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTggg                                                                                                           >  2:7278385/1‑150 (MQ=255)
                                             aGTCAAGTGGAAGCTGCACTGACCGAAGCGCTCAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGt                                                                                                          >  2:1092864/1‑150 (MQ=255)
                                                 aaGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCaa                                                                                                      <  2:674591/150‑1 (MQ=255)
                                                 aaGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCaa                                                                                                      <  2:4625831/150‑1 (MQ=255)
                                                   ggggAAGCTGCCCTGACCGAAGCGCTGAATTCCCGCGGTAATTAAGTACAGGGCAAGGCCCCGGCACGCCGCGGGCGATCATCGGCGACATGAACTGTGGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGa                                                                                                    <  1:3729104/148‑1 (MQ=255)
                                                    tGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGAt                                                                                                   >  1:5676215/1‑150 (MQ=255)
                                                     ggAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATc                                                                                                  >  2:6977703/1‑150 (MQ=255)
                                                      gAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATcc                                                                                                 >  2:3215399/1‑150 (MQ=255)
                                                          cTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCACGATCCTCtt                                                                                             >  1:2802159/1‑150 (MQ=255)
                                                            gCACTGACCGAAGCGCTTAATTCCCGCGGTAATTAAGTCCAGGGCAAGGCCCCGGCAAGCCTCGGGCTATAATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTa                                                                                           <  1:5023266/150‑1 (MQ=255)
                                                            gCACTGACCGAAGCGCTTAATACCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGAAATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTa                                                                                           <  1:6292849/150‑1 (MQ=255)
                                                            gCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTa                                                                                           >  2:6615497/1‑150 (MQ=255)
                                                                tGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACt                                                                                       >  2:1460395/1‑150 (MQ=255)
                                                                   ccGAAGCGATGAATTCCAGCGGTAATTAAGTACTGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACtttt                                                                                    <  1:4709490/150‑1 (MQ=255)
                                                                     gAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGt                                                                                  >  1:4322973/1‑150 (MQ=255)
                                                                      aaGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTa                                                                                 <  2:4220888/150‑1 (MQ=255)
                                                                            tGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACtt                                                                                      <  2:2880800/139‑1 (MQ=255)
                                                                            tGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTa                                                                           >  2:5436016/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)
                                                                                   ctgCGGTAATTCAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTgg                                                                    <  2:7152807/148‑1 (MQ=255)
                                                                                   cAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTgg                                                                    <  2:858919/150‑1 (MQ=255)
                                                                                   cAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTgg                                                                    <  1:3441773/150‑1 (MQ=255)
                                                                                   cAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTgg                                                                    <  2:387214/150‑1 (MQ=255)
                                                                                      cGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGa                                                                 <  2:1221340/150‑1 (MQ=255)
                                                                                      cGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGa                                                                 <  2:4299065/150‑1 (MQ=255)
                                                                                        gTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGCCATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAgg                                                               >  1:1872478/1‑150 (MQ=255)
                                                                                        gTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAgg                                                               <  2:6171200/150‑1 (MQ=255)
                                                                                         tacttAAGTACAGGGCAAGGCCCCGGCACGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCGTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGc                                                              >  2:1600972/4‑150 (MQ=255)
                                                                                           aTTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTc                                                            >  2:5889743/1‑150 (MQ=255)
                                                                                                 tACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGtctc                                                      >  1:1623817/1‑150 (MQ=255)
                                                                                                   cAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGc                                                    >  1:5585145/1‑150 (MQ=255)
                                                                                                       gCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGt                                                                   <  1:6016682/131‑1 (MQ=255)
                                                                                                           ggCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTGGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAAttttt                                            <  2:2819359/150‑1 (MQ=255)
                                                                                                            gCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTg                                           <  1:4382073/150‑1 (MQ=255)
                                                                                                             ccccGGCACGCCTCGGGCGATCTTCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTgg                                          <  1:806722/150‑1 (MQ=255)
                                                                                                             ccccGGCAAGCCTCGGGCGATCATCGGCGACCGGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGGGCAAGATCCGCGTGAGACTTTTGTACCCTTAAATTTGGTGAGGCGCCGTCTCGCCAAGTGTTgg                                          >  1:6257477/1‑150 (MQ=255)
                                                                                                                 ggCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCttt                                      >  1:3882026/1‑150 (MQ=255)
                                                                                                                      gCCTCGGGCGATCATCGGCGACATTAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAGTTTTTGGCTTTCTTGc                                 >  2:1312962/1‑150 (MQ=255)
                                                                                                                        cTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTATTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCaa                               <  2:3075150/150‑1 (MQ=255)
                                                                                                                         tCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTa                                                                                           >  1:5726779/1‑89 (MQ=255)
                                                                                                                         tCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCGTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAg                              >  1:7000072/1‑150 (MQ=255)
                                                                                                                            ggCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACa                           <  2:5322154/150‑1 (MQ=255)
                                                                                                                            ggCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACa                           >  2:2695838/1‑150 (MQ=255)
                                                                                                                            ggCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGAGCAAGATCCTCTATAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACa                           <  2:6354385/150‑1 (MQ=255)
                                                                                                                               gATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATa                        >  1:4637369/1‑150 (MQ=255)
                                                                                                                               gATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATa                        <  1:539580/150‑1 (MQ=255)
                                                                                                                               gATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATa                        <  2:1791929/150‑1 (MQ=255)
                                                                                                                                 tcatcGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATaga                      <  2:1349199/150‑1 (MQ=255)
                                                                                                                                    tcGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTc                   <  1:3646309/150‑1 (MQ=255)
                                                                                                                                    tcGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTc                   <  2:7108081/150‑1 (MQ=255)
                                                                                                                                      ggCGACATTAACTGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGc                 <  2:1855461/150‑1 (MQ=255)
                                                                                                                                       gCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCa                >  2:5605097/1‑150 (MQ=255)
                                                                                                                                         gACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCTATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAg              <  1:2949375/150‑1 (MQ=255)
                                                                                                                                          aCATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGTTGAAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGATCAGTCTCGCCAATTTTTGTCTTTCTTGAAAGACAATATATTCTCAAGa             >  1:5878918/1‑150 (MQ=255)
                                                                                                                                            aTGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGACa           >  1:2200294/1‑150 (MQ=255)
                                                                                                                                                 aTGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGACAATCAg      <  1:1743793/150‑1 (MQ=255)
                                                                                                                                                 aTGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGACAATCAg      <  2:3240908/150‑1 (MQ=255)
                                                                                                                                                  tGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGACAATCAGt     <  1:7219607/150‑1 (MQ=255)
                                                                                                                                                  tGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGACAATCAGt     <  1:1456255/150‑1 (MQ=255)
                                                                                                                                                  tGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGACAATCAGt     <  2:4965462/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGACAATCAGTGGA  >  NC_002947/545868‑546166

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: