Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,453,536 C→G 11.9% G1032G (GGC→GGG nicB → nicotinate dehydrogenase subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,453,5360CG11.9% 192.3 / 15.3 67G1032G (GGC→GGGnicBnicotinate dehydrogenase subunit B
Reads supporting (aligned to +/- strand):  ref base C (28/31);  new base G (5/3);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

CCTTGCTGGCGGCATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCGCTGAGCAACG  >  NC_002947/4453389‑4453678
                                                                                                                                                   |                                                                                                                                              
ccTTGCTGGCGGCATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGccc                                                                                                                                              <  2:453291/150‑1 (MQ=255)
 cTTGCTGGCGGCATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCATGGGGCCgg                                                                                                                                            >  2:2242620/1‑150 (MQ=255)
   tGCTGGCGGCATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCATGGGGCCggcc                                                                                                                                          <  1:2242620/150‑3 (MQ=255)
    gCTGGCGGCATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCggtg                                                                                                                                          <  1:874740/150‑1 (MQ=255)
    gCTGGCGGCATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGCAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCggtg                                                                                                                                          <  1:2369559/150‑1 (MQ=255)
      tGGCGGCATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGTATCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCggtggt                                                                                                                                        <  1:1016243/150‑1 (MQ=255)
        gcggcATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCa                                                                                                                                      <  2:1461094/150‑1 (MQ=255)
         cggcATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGCCCCATGGGCCCGGTGGTCAGt                                                                                                                                    <  1:1198978/150‑1 (MQ=255)
                     tGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGGCAc                                                                                                                         >  2:728572/1‑150 (MQ=255)
                        cTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCt                                                                                                                      <  2:820938/150‑1 (MQ=255)
                         tGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGccc                                                                                                                                              <  1:2037073/125‑1 (MQ=255)
                         tGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGccc                                                                                                                                              >  2:2037073/1‑125 (MQ=255)
                          gggAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGc                                                                                                                    >  1:896014/1‑150 (MQ=255)
                                ccccGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGTGTAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCaaaaa                                                                                                              <  1:744840/150‑1 (MQ=255)
                                acccGGCCATGAATGCCATGGGCAAGGCATCGACGCCGTGGAGCGAAGACGAGTGGTGCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTTGTCAGTGAGCTTGCCACCCTGCCaaaaa                                                                                                              >  2:279800/2‑150 (MQ=255)
                                 cccGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCAc                                                                                                                         <  2:1570791/138‑1 (MQ=255)
                                 cccGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCAc                                                                                                                         >  1:1570791/1‑138 (MQ=255)
                                 cccGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGGTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTCCCAAAAAg                                                                                                             >  2:608524/1‑150 (MQ=255)
                                   cGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGAGAGGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCg                                                                                                           <  1:724532/150‑1 (MQ=255)
                                     gCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGAc                                                                                                         <  1:1783243/150‑1 (MQ=255)
                                     gCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGAc                                                                                                         <  1:1688943/150‑1 (MQ=255)
                                        ttGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGAc                                                                                                         <  2:2104644/147‑1 (MQ=255)
                                        ttGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGAc                                                                                                         >  1:2104644/1‑147 (MQ=255)
                                           aaTGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCa                                                                                                                                      >  2:1319046/1‑115 (MQ=255)
                                           aaTGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCa                                                                                                                                      <  1:1319046/115‑1 (MQ=255)
                                           aaTGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTgcgc                                                                                                   <  2:896014/150‑1 (MQ=255)
                                            aTGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCggtg                                                                                                                                          >  1:1658997/1‑110 (MQ=255)
                                            aTGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCggtg                                                                                                                                          <  2:1658997/110‑1 (MQ=255)
                                                ccTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCAt                                                                                              >  1:2422583/1‑150 (MQ=255)
                                                          tCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTAc                                                                                    >  2:2243881/1‑150 (MQ=255)
                                                          tCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTAc                                                                                    <  1:906956/150‑1 (MQ=255)
                                                                          gCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAgg                                                                    <  2:1866305/150‑1 (MQ=255)
                                                                            gAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCTTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGc                                                                  >  2:631504/1‑150 (MQ=255)
                                                                              aGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGa                                                                >  2:2033549/1‑150 (MQ=255)
                                                                              aGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGa                                                                <  2:2009948/150‑1 (MQ=255)
                                                                              aGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCACCCACTACCTCAGTTCCCTGGAGGGCGa                                                                <  1:728572/150‑1 (MQ=255)
                                                                                  gAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGTGCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCg                                                            >  2:576899/1‑150 (MQ=255)
                                                                                  gAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCg                                                            >  1:1022205/1‑150 (MQ=255)
                                                                                    gttgttCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTTGAGGGCGAGCCGCa                                                          <  2:943302/150‑1 (MQ=255)
                                                                                        ttCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCCCCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCa                                                      >  2:385000/1‑150 (MQ=255)
                                                                                        ttCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAgg                                                                    >  2:245169/1‑136 (MQ=255)
                                                                                        ttCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAgg                                                                    <  1:245169/136‑1 (MQ=255)
                                                                                             cTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCGCTGGAGGGCGAGCCGCAGGCACTGgc                                                 >  1:2206521/1‑150 (MQ=255)
                                                                                             cTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGgc                                                 <  2:1022205/150‑1 (MQ=255)
                                                                                                 cTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGgccgcc                                             >  2:1947916/1‑150 (MQ=255)
                                                                                                  tCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCa                                            >  2:551027/1‑150 (MQ=255)
                                                                                                       aCCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCCCTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGcc                                       <  1:2134446/150‑1 (MQ=255)
                                                                                                        ccGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTccc                                                                          <  2:2326615/114‑1 (MQ=255)
                                                                                                        ccGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTccc                                                                          >  1:2326615/1‑114 (MQ=255)
                                                                                                         cGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTccc                                                                          <  2:1724407/113‑1 (MQ=255)
                                                                                                         cGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTccc                                                                          >  1:1724407/1‑113 (MQ=255)
                                                                                                         cGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCcgc                                     <  1:2243881/150‑1 (MQ=255)
                                                                                                         cGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCcgc                                     >  2:1708757/1‑150 (MQ=255)
                                                                                                            cTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGcc                                  <  2:1984921/150‑1 (MQ=255)
                                                                                                                tcCGAAAAACACGGCGTGGCGGCCGGCCCCATTGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAg                              >  2:784362/1‑150 (MQ=255)
                                                                                                                       aaCACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGaca                       >  1:2174861/1‑150 (MQ=255)
                                                                                                                        acacGGCGTGGCGGCCAGCCCCATGGGCCCGGTAGTCAGTGAGCTTGACACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGacac                      >  2:481226/1‑150 (MQ=255)
                                                                                                                         cacGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCATGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGAcacg                     >  2:1165125/1‑150 (MQ=255)
                                                                                                                          acGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGAcacgc                    >  2:1856791/1‑150 (MQ=255)
                                                                                                                            ggCGTGGCGGCCGGCCCCATGGGCCCGGTTGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGAcacgcac                  <  1:279800/150‑1 (MQ=255)
                                                                                                                                   cggccggcCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCGc           <  2:2174861/150‑1 (MQ=255)
                                                                                                                                    caccgggcCCATGGGGCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAg                                                         <  2:1643201/101‑1 (MQ=255)
                                                                                                                                    caccgggcCCATGGGGCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAg                                                         >  1:1643201/3‑103 (MQ=255)
                                                                                                                                    caccgggcCCATGGGGCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCGCt          >  1:126012/3‑150 (MQ=255)
                                                                                                                                        ggccCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACGGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCGCTGAGc      <  1:178334/150‑1 (MQ=255)
                                                                                                                                         gccCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGAc                                                                                                         >  1:2116318/1‑50 (MQ=255)
                                                                                                                                         gccCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGAc                                                                                                         <  2:2116318/50‑1 (MQ=255)
                                                                                                                                            ccATGGGGCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTg                           >  2:2081158/1‑125 (MQ=255)
                                                                                                                                            ccATGGGGCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTg                           <  1:2081158/125‑1 (MQ=255)
                                                                                                                                            ccATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCGCTGAGCAACg  >  2:1006066/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CCTTGCTGGCGGCATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCGCTGAGCAACG  >  NC_002947/4453389‑4453678

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