Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 543,101 T→C 12.1% A446A (GCT→GCC rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947543,1010TC12.1% 954.1 / 91.3 347A446A (GCT→GCCrpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base T (167/137);  new base C (21/21);  total (189/158)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

Only 100 of 351 total aligned reads displayed.

GCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGG  >  NC_002947/542952‑543250
                                                                                                                                                     |                                                                                                                                                     
   aagaagATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTAtt                                                                                                                                                    >  2:8401352/1‑150 (MQ=255)
   aagaagATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTAtt                                                                                                                                                    >  2:8720096/1‑150 (MQ=255)
   aagaagATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTACTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTAtt                                                                                                                                                    >  1:3274512/1‑150 (MQ=255)
      aagaTGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAg                                                                                                                                                 >  1:2365866/1‑150 (MQ=255)
        gaTGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCt                                                                                                                                               >  2:6534426/1‑150 (MQ=255)
              cctgcgcGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCAccc                                                                                                                                         >  2:7500543/4‑150 (MQ=255)
               gAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCg                                                                                                                                        >  2:6997459/1‑150 (MQ=255)
               gAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTGTGAACCCGTACTGATCGAAGGTACGGCTATTCAGCTGCACCCg                                                                                                                                        >  2:5338609/1‑150 (MQ=255)
               gAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCCCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGATGGTAAGGCTATTCAGCTGCACCCg                                                                                                                                        >  2:1235249/1‑150 (MQ=255)
                aGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCACAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGc                                                                                                                                       <  2:2599108/150‑1 (MQ=255)
                  cgcgAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTg                                                                                                                                     <  1:4262033/150‑1 (MQ=255)
                  cgcgAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTg                                                                                                                                     <  1:8540229/150‑1 (MQ=255)
                        cTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCtgt                                                                                                                               <  1:7934676/150‑1 (MQ=255)
                           ccAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGtgct                                                                                                                            <  1:4507362/150‑1 (MQ=255)
                            cAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGtgctg                                                                                                                           <  2:283234/150‑1 (MQ=255)
                                     gggACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACg                                                                                                                  >  2:8012895/1‑150 (MQ=255)
                                           ttCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGATTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACt                                                                                                            <  1:7312812/150‑1 (MQ=255)
                                                gCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGAc                                                                                                       >  2:6397998/1‑150 (MQ=255)
                                                gCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGAc                                                                                                       <  1:505181/150‑1 (MQ=255)
                                                  tGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACgg                                                                                                     >  2:4823079/1‑150 (MQ=255)
                                                   gAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGt                                                                                                    <  2:5115798/150‑1 (MQ=255)
                                                   gAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGt                                                                                                    >  1:2159319/1‑150 (MQ=255)
                                                   gAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGt                                                                                                    >  1:6353955/1‑150 (MQ=255)
                                                    aaGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTg                                                                                                   >  1:1988114/1‑150 (MQ=255)
                                                     aGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGa                                                                                                  >  1:5672038/1‑150 (MQ=255)
                                                     aGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGa                                                                                                  >  1:7644041/1‑150 (MQ=255)
                                                      gTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGAc                                                                                                 <  1:3626249/150‑1 (MQ=255)
                                                      gTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGAc                                                                                                 <  1:1752356/150‑1 (MQ=255)
                                                        gATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCa                                                                                               <  2:7168046/150‑1 (MQ=255)
                                                           tCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCt                                                                                                                                               >  1:2041825/1‑99 (MQ=255)
                                                           tCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCt                                                                                                                                               <  2:2041825/99‑1 (MQ=255)
                                                               gAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGcc                                                                                        <  2:6245172/150‑1 (MQ=255)
                                                                aaCACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCg                                                                                       >  2:3056165/1‑150 (MQ=255)
                                                                aaCACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCg                                                                                       <  1:8381268/150‑1 (MQ=255)
                                                                    ccccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCa                                                                                   <  2:4021775/150‑1 (MQ=255)
                                                                     cccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCCGATGGCCGTTCAc                                                                                  >  2:7369168/1‑150 (MQ=255)
                                                                     cccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCAc                                                                                  >  2:1570882/1‑150 (MQ=255)
                                                                     cccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCAc                                                                                  >  2:7256355/1‑150 (MQ=255)
                                                                     cccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCAc                                                                                  >  2:1358884/1‑150 (MQ=255)
                                                                     cccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCAc                                                                                  >  2:8670476/1‑150 (MQ=255)
                                                                      ccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACg                                                                                 >  1:691268/1‑150 (MQ=255)
                                                                         tACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGc                                                                              >  2:5222082/1‑150 (MQ=255)
                                                                            tGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGc                                                                           >  2:3464222/1‑150 (MQ=255)
                                                                            tGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGc                                                                           >  1:1447257/1‑150 (MQ=255)
                                                                             gCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCt                                                                          <  1:1581114/150‑1 (MQ=255)
                                                                              cTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTg                                                                         >  1:6102802/1‑150 (MQ=255)
                                                                               tgaaCCGCGCGCCTACGCTGCACCGTCTTGGTATCCAGGCCTTCGAACCCACGCTGATCGAAGGCAAGGCTATCCAGCTGCATCCGCTCGTTTGTACGGCGTTCAACGCCGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCGCTgt                                                                        >  1:4447194/3‑149 (MQ=255)
                                                                               tCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGa                                                                        >  2:5723498/1‑150 (MQ=255)
                                                                                 aaCCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGcc                                                                                                                >  1:1015167/1‑108 (MQ=255)
                                                                                    cGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTg                                                                   <  2:1750915/150‑1 (MQ=255)
                                                                                    cGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTg                                                                   >  1:4693311/1‑150 (MQ=255)
                                                                                     gTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTgg                                                                  >  1:5916232/1‑150 (MQ=255)
                                                                                         aCCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGc                                                              >  2:8603734/1‑150 (MQ=255)
                                                                                          ccGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGcc                                                             >  2:7746309/1‑150 (MQ=255)
                                                                                           cGACCCTTCACCGTCTGGGTATCCAGGCGTTTGCACCCGTACGGATCGAAGGGAAGGCGATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTGCGACGGTGACCAGATGGCCGTGCACGTGCCGCGGACCCTGGAAGccc                                                            >  2:4533899/1‑150 (MQ=255)
                                                                                           cGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGccc                                                            >  1:7578155/1‑150 (MQ=255)
                                                                                            tacgctgcACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCa                                                           >  2:8787265/8‑150 (MQ=255)
                                                                                                 ttCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCg                                                      >  2:8546242/1‑150 (MQ=255)
                                                                                                 tgcACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCttg                                                      >  1:2908040/3‑148 (MQ=255)
                                                                                                  tCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGa                                                     >  1:7115275/1‑150 (MQ=255)
                                                                                                  tCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGa                                                     >  2:3726177/1‑150 (MQ=255)
                                                                                                   cACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGaa                                                    <  2:2365866/150‑1 (MQ=255)
                                                                                                   cACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGaa                                                    <  1:1339466/150‑1 (MQ=255)
                                                                                                    aCCGTCTTGGTATTAAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAg                                                   >  1:8743810/1‑150 (MQ=255)
                                                                                                    aCCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAg                                                   >  2:4323448/1‑150 (MQ=255)
                                                                                                       gTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAgcac                                                <  1:3671650/150‑3 (MQ=255)
                                                                                                       gTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAgcgc                                                <  1:1053012/150‑1 (MQ=255)
                                                                                                       gTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCAAGCTCGAAgcgc                                                <  1:5508269/150‑1 (MQ=255)
                                                                                                           tggTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTCAAGCTGAGCTTCAAgcgacaga                                            <  2:4202185/149‑6 (MQ=255)
                                                                                                           gggTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTgc                                            <  2:3924932/150‑1 (MQ=255)
                                                                                                            ggTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGcacgcgtg                                           >  2:3553668/1‑146 (MQ=255)
                                                                                                            ggTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTgcg                                           >  1:7229292/1‑150 (MQ=255)
                                                                                                            ggTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTgcg                                           >  2:6717058/1‑150 (MQ=255)
                                                                                                             gTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGcacgcgtgc                                          >  1:4118657/1‑145 (MQ=255)
                                                                                                                ttcAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCtga                                       >  2:7673250/3‑150 (MQ=255)
                                                                                                                tCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtga                                       >  2:7809504/1‑150 (MQ=255)
                                                                                                                    ggCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTgatgat                                   >  2:8187990/1‑150 (MQ=255)
                                                                                                                    ggCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTgatgat                                   >  2:7958057/1‑150 (MQ=255)
                                                                                                                    ggCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTgatgat                                   >  2:3342892/1‑150 (MQ=255)
                                                                                                                     gCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTACACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGatgatg                                  <  2:7069856/150‑1 (MQ=255)
                                                                                                                       gTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTc                                >  1:8942356/1‑150 (MQ=255)
                                                                                                                       gTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTc                                <  1:4017405/150‑1 (MQ=255)
                                                                                                                       gTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTc                                >  1:6238306/1‑150 (MQ=255)
                                                                                                                        tttGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCg                               <  2:7142128/150‑1 (MQ=255)
                                                                                                                               ccacgcTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACcaaca                        <  2:5086254/145‑1 (MQ=255)
                                                                                                                               ccGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACcaaca                        >  1:3165356/1‑150 (MQ=255)
                                                                                                                                 acgcTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCagcaac                      >  1:2669110/4‑150 (MQ=255)
                                                                                                                                    cTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATc                   <  1:8165478/150‑1 (MQ=255)
                                                                                                                                    cTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATc                   <  2:1022023/150‑1 (MQ=255)
                                                                                                                                       aTCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGt                                              <  2:7267715/120‑1 (MQ=255)
                                                                                                                                       aTCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTg                >  1:4080537/1‑150 (MQ=255)
                                                                                                                                        tCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACACCATCCTGt               <  1:6701418/150‑1 (MQ=255)
                                                                                                                                        tCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGt               <  1:5200165/150‑1 (MQ=255)
                                                                                                                                         cGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTc              <  1:4647927/150‑1 (MQ=255)
                                                                                                                                           aaGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgc            <  2:8960063/150‑1 (MQ=255)
                                                                                                                                            aGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgcc           <  1:958583/150‑1 (MQ=255)
                                                                                                                                             ggTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGc                                                         >  1:5984654/1‑103 (MQ=255)
                                                                                                                                             ggTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGc                                                         <  2:5984654/103‑1 (MQ=255)
                                                                                                                                             ggTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgcca          <  1:7978853/150‑1 (MQ=255)
                                                                                                                                                  ggCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACACCGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCaa     >  2:2214425/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGG  >  NC_002947/542952‑543250

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