Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,531,992 T→A 12.4% intergenic (+32/‑136) PP_4018 → / → topB acetyltransferase/DNA topoisomerase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,531,9920TA12.4% 387.2 / 37.9 148intergenic (+32/‑136)PP_4018/topBacetyltransferase/DNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base T (76/51);  new base A (9/9);  total (86/62)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 149 total aligned reads displayed.

GCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGTTT  >  NC_002947/4531847‑4532139
                                                                                                                                                 |                                                                                                                                                   
 cccAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCcg                                                                                                                                                <  2:3259042/150‑1 (MQ=255)
  ccAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTTCCcgc                                                                                                                                               >  1:3542196/1‑150 (MQ=255)
  ccAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCcgc                                                                                                                                               <  2:2405735/150‑1 (MQ=255)
  ccAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCcgc                                                                                                                                               <  1:2084485/150‑1 (MQ=255)
   cAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCTTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCGCCTGAGCCGGCCTCTTCGCGGGTAAACCcgcg                                                                                                                                              <  2:572145/150‑1 (MQ=39)
    aGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGa                                                                                                                                             >  2:2907510/1‑150 (MQ=255)
    aGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGa                                                                                                                                             >  2:2907467/1‑150 (MQ=255)
    aGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGa                                                                                                                                             <  2:1436286/150‑1 (MQ=255)
     ggTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGaa                                                                                                                                            <  1:3354844/150‑1 (MQ=255)
     ggTGCAAGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa                                                                                                                                            <  2:852102/150‑1 (MQ=255)
      gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAag                                                                                                                                           >  2:1763491/1‑150 (MQ=255)
       tGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAaga                                                                                                                                          <  2:1261433/150‑1 (MQ=255)
       tGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAaga                                                                                                                                          <  2:1356203/150‑1 (MQ=255)
       tGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACCGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                          <  1:3564885/150‑1 (MQ=255)
        gCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAagag                                                                                                                                         <  1:3465906/150‑1 (MQ=255)
         cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTGTACCCGCGAAGAgg                                                                                                                                        >  1:3922779/1‑150 (MQ=255)
          acgccacgccGTTCTATGGGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGc                                                                                                                                       >  1:268285/1‑150 (MQ=255)
           cgccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGcc                                                                                                                                      <  2:1588516/150‑1 (MQ=255)
            gccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCg                                                                                                                                     <  2:2981667/150‑1 (MQ=255)
             ccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTAGCGGGTTTAACCGCGAAGAGGCCgg                                                                                                                                    >  1:4193189/1‑150 (MQ=255)
              cacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGc                                                                                                                                   >  1:2934213/1‑150 (MQ=255)
                 gccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCa                                                                                                                                >  2:2841924/1‑150 (MQ=255)
                 gccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCa                                                                                                                                >  1:2534296/1‑150 (MQ=255)
                 gccGTTCTAAGAGCGGCTGGGCTTTCGCGTAGTCAGAGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCa                                                                                                                                >  2:3168281/1‑150 (MQ=255)
                        tATGAGCGGATGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAATATGGCATCCCGCATGTCGACATGGTGAGAGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGCGGGTAAACCCCAAAAGCCCGGCTCACGCGcta                                                                                                                         >  1:224710/1‑150 (MQ=255)
                             gCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACAt                                                                                                                    >  2:4088874/1‑150 (MQ=255)
                                   gggCTTTCGCGTTGTCAGCGAGGAATTCCTCGACGCTGGCATCCCGCTTGTCGCCCTGGGGCGCGCCTCGGGCGCGTGATCCATCACCTGCACTGGCCTCTTCGGGGGTTTCCCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGc                                                                                                              <  2:2511217/150‑1 (MQ=255)
                                       tttCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGAc                                                                                                          >  2:1798549/1‑150 (MQ=255)
                                       tttCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGAc                                                                                                          >  1:3255131/1‑150 (MQ=255)
                                       tttCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGAc                                                                                                          >  2:885330/1‑150 (MQ=255)
                                           gcgTAGTCAGTGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACaaaa                                                                                                      >  2:1758581/1‑150 (MQ=255)
                                           gcgTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGAGATACATCCCCGCTGACaaaa                                                                                                      >  2:2885122/1‑150 (MQ=255)
                                                 tCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTg                                                                                                >  2:4441987/1‑150 (MQ=255)
                                                  cAGCGTGTAATTCCTCGAAGGTGGCATCCCGAATGTCGACTTGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGCCCTCTTCGCGGTTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGt                                                                                               <  2:3542196/150‑1 (MQ=255)
                                                   tgcttgtaaTCCCTCGAAGCGGGCACCCCGCTTGTCGACTGGGGGCGCGACTCGCGCGCCGGACCCAACACCTCCAGGGGCCTCTGCGGGGTTTTACCCGCTAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTa                                                                                              <  2:2404666/143‑1 (MQ=255)
                                                      gAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCAGTAcat                                                                                           >  2:3828248/1‑150 (MQ=255)
                                                      gAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCATCGCGGGGTGACCCGCGAAGAGGCCGGCTCAGGCGATATACCTCCCCGCTGACAAAAAACCAGTAcat                                                                                           >  1:2253824/1‑150 (MQ=255)
                                                         gAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACAACCCCGCAGACAAAAAACCAGTAcaacca                                                                                        >  2:1627787/1‑150 (MQ=255)
                                                           agtCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccatc                                                                                      <  2:4121482/148‑1 (MQ=255)
                                                            ttCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTTCCCGCGAAGAGGCCGGCGCAGGCGATATACATCCCCCCTGACAAAAAACCTGTAcagacaacc                                                                                     >  1:3834538/1‑143 (MQ=255)
                                                              ccTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAg                                                                                   >  2:2508156/1‑150 (MQ=255)
                                                                tCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTAGCGGGTTTACCCGCGAACAGGCCGGCTCAGGCGATATACATCCCCGCTGACCAAAAAACTGTACATCCAACCAGAt                                                                                 >  1:110692/1‑150 (MQ=255)
                                                                  gAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTGGCGGGTTTAACCGCGAAGAGGGCGGCTCAGGCGGTTTTAATACCCGCTTACAAAAAAACTGTAAATACATACACATaa                                                                               >  1:3038543/1‑150 (MQ=255)
                                                                    aggtGGAATCCCGCATGTCGACTTGGTGGGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATaaaa                                                                             <  2:3340419/147‑1 (MQ=255)
                                                                    aGCTGGCATCCCGCTTGTCGACTTGGTGCGCGACTGGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGGGGGTTTACCCGCGAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATaaaa                                                                             <  2:660617/150‑1 (MQ=255)
                                                                       tGGCACCCCGATTGTCGAATTGGGGCGCTACTCGGGCGCGTTATCAATAACCTGCACTGGCCTATTCGGGGTTTTACCCGCTAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACtt                                                                          <  2:443916/150‑1 (MQ=255)
                                                                            tCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATACCCCCTGACAAAAAAACAGTAAAACCAACCAGATAAAACTTTCctc                                                                     >  1:1269853/1‑150 (MQ=255)
                                                                              ccGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTTCCCGCCAAGAGGCCGGCGCAGGCGAAATACATCCCCCCGGACAAAAAACACGTACAACCCTCACCATAAAAATTGCCTCTg                                                                   >  1:2772053/1‑150 (MQ=255)
                                                                                    gTCGACTGGGGGCGCGACTCGCGCGCCTGAACCACCACCTGCACGGGCCTCTGCGGGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGcccc                                                             <  1:2976116/150‑1 (MQ=255)
                                                                                         gtggttGCGCGCGTCGGGCGCCTGATCCAGCAGCTGCAGGGGCCTGTTGGCTGGTTGACCCGGGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTTCGCCCCCGCtt                                                        <  2:3993097/145‑1 (MQ=255)
                                                                                         cATGGGGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCCGCGCCCCCGCtt                                                        <  2:2601936/150‑1 (MQ=255)
                                                                                                  cgACTCGCGCGCCTGATCCTTCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCt                                               <  2:320131/150‑1 (MQ=255)
                                                                                                   gACTCGCGCGCGTTAACAAACACGTGCACTGGCCTCTTGGGGGTTTTCCCGGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTa                                              <  1:2901632/150‑1 (MQ=255)
                                                                                                   gACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTTCCCGCGAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTa                                              >  1:1489502/1‑150 (MQ=255)
                                                                                                   gACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACATGCCTCTGCGCCCCCGCTTGCGCATCCTa                                              >  2:1009742/1‑150 (MQ=255)
                                                                                                   gACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCCGGTTTACCCCCGAAGAGGCCGGCTCAGGCGATATAAAAACCCGCTGACAAAAAACCAGTAAATCCATACAGATAAAAAATTCCACTGCGCCCCCCCTTGCGCATccaa                                              >  1:2039560/1‑148 (MQ=255)
                                                                                                      tCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATgg                                           >  2:3189150/1‑150 (MQ=255)
                                                                                                      tCGCGCGCCTGATCCATCACCTGCACTGGCCTCGTCGCGGGTTTACCCGCGCAGAGGCCGGCGCAGGCGCTATACATCCCCGCGGACCAAAAACCGGGACATCCATCCCGCTAAAACTTGCCGCGGCGCCCCCGCGTGCGCAGCCGATgg                                           >  1:2199803/1‑150 (MQ=255)
                                                                                                          gcgcCTGATCCACCACCGGCACGGGCCTCTGCGGGGGTTTCCCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCcca                                       <  2:2052372/150‑1 (MQ=255)
                                                                                                               tGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCaccacc                                  >  1:1465497/1‑150 (MQ=255)
                                                                                                                gctCCATCCCCTGCACGGGCCTCTTCGCGGGTTTACCGGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCAccc                                 <  2:3267471/148‑1 (MQ=255)
                                                                                                                gATCCATCACCTGCACTGGCCTCTTCGCGGGTTTGCCCGCGCAGAGGCCGGCGCAGGCGATAAACAACCCCGCGGACCAAAAACCCGTACAACCAACCAGATAAAACGTGCCGCGGCGCCCCCGCGTGCGCCGCCTAAGGcccccccccc                                 >  1:2360708/1‑150 (MQ=255)
                                                                                                                  tCCATCACCTGCACGGGCCTCTTGGGGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCg                               <  2:657578/150‑1 (MQ=255)
                                                                                                                    cATCACCTGCAGTGGCATCTTCGCGGTTTTACCCGCGAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGAt                             <  2:2876658/150‑1 (MQ=255)
                                                                                                                      tccccTGCACTGGCCTCTGCGGGGGTTTCCCCGCGATGCGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGa                           <  2:4147262/147‑1 (MQ=255)
                                                                                                                       cACCTGCAGTGGCCTGTTGGGGTTTTTACCCGCTAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGaa                          <  2:1993407/150‑1 (MQ=255)
                                                                                                                       cACCTGCACTGGCCTCTTCGGGGTTTTACCCGCTAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGaa                          <  2:4045453/150‑1 (MQ=255)
                                                                                                                                gggCCGCTGCGGGGGTTTCCCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCg                 <  1:2248598/149‑1 (MQ=255)
                                                                                                                                  gCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAAACTGTACATCCATCCAGATAAAAATTGCCTCTGCGCCCCCGCTTGCGCATCCAATGGCCCACCACCCCGATGAAGCGTCCCCCGc               >  1:1925818/1‑150 (MQ=255)
                                                                                                                                   ccTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGACGCGTCCCCGGcc              >  2:1860675/1‑150 (MQ=255)
                                                                                                                                   ccTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGcc              >  1:3182922/1‑150 (MQ=255)
                                                                                                                                    ctctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCa             >  2:3970111/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt            >  2:1920236/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt            >  1:3479544/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt            >  2:1039496/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt            >  1:3174268/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGGTAAACCCGCGAAGAGGCCAGTGCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt            >  2:3816380/1‑150 (MQ=39)
                                                                                                                                       ttCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc          >  1:787995/1‑150 (MQ=255)
                                                                                                                                       ttCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc          >  1:1118527/1‑150 (MQ=255)
                                                                                                                                       ttCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc          >  1:787946/1‑150 (MQ=255)
                                                                                                                                       ttCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc          >  2:3438460/1‑150 (MQ=255)
                                                                                                                                       ttCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc          >  1:215217/1‑150 (MQ=255)
                                                                                                                                       ttCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGcc                                         >  2:573012/1‑119 (MQ=255)
                                                                                                                                       ttCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGcc                                         <  1:573012/119‑1 (MQ=255)
                                                                                                                                       ttCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc          >  2:4479340/1‑150 (MQ=255)
                                                                                                                                       ttCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc          >  1:1738174/1‑150 (MQ=255)
                                                                                                                                        tCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGaa                          <  1:2815640/133‑1 (MQ=255)
                                                                                                                                        tCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGaa                          >  2:2815640/1‑133 (MQ=255)
                                                                                                                                         cgcgGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgcgc        >  2:2236108/1‑150 (MQ=255)
                                                                                                                                          gcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGcc       >  2:4499773/1‑150 (MQ=255)
                                                                                                                                          gcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGcc       >  1:1419635/1‑150 (MQ=255)
                                                                                                                                           ggggTTTCCCCGCGACGCGGCCGGCTCAGGCGATATACCTCCCCGCTGACAAAAAACCTGTACAGCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCt      <  2:3495401/149‑1 (MQ=255)
                                                                                                                                           cgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCt      >  2:1676957/1‑150 (MQ=255)
                                                                                                                                           cgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCt      <  1:4119681/150‑1 (MQ=255)
                                                                                                                                            gggTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTg     >  1:3069954/1‑150 (MQ=255)
                                                                                                                                            gggTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTg     >  1:3030189/1‑150 (MQ=255)
                                                                                                                                            gggTTCCCCCGCGACGAGCCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTg     <  1:2418976/150‑1 (MQ=255)
                                                                                                                                             ggTTTACCCGCGAAGCGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGt    <  2:2651487/150‑1 (MQ=255)
                                                                                                                                             ggTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGt    >  1:4481222/1‑150 (MQ=255)
                                                                                                                                               tttCCCCGCCACGCGCCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGttt  <  1:960064/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGTTT  >  NC_002947/4531847‑4532139

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: