Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 43,461:1 +A 13.1% intergenic (+27/+66) PP_0040 → / ← cadA‑I hypothetical protein/cadmium translocating P‑type ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294743,4611.A13.1% 281.4 / 42.4 84intergenic (+27/+66)PP_0040/cadA‑Ihypothetical protein/cadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base . (34/39);  new base A (6/5);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

CGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATCGCAACAG  >  NC_002947/43317‑43606
                                                                                                                                                 |                                                                                                                                                 
cGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCa                                                                                                                                              <  1:2316672/150‑1 (MQ=255)
cGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCa                                                                                                                                              <  1:1490006/150‑1 (MQ=255)
 ggTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCAc                                                                                                                                             >  2:2050190/1‑150 (MQ=255)
 ggTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCAc                                                                                                                                             <  1:1328058/150‑1 (MQ=255)
   tttCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCGATCAGCGCACGa                                                                                                                                           <  2:1165126/150‑1 (MQ=255)
    ttCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGAc                                                                                                                                          <  2:2297605/150‑1 (MQ=255)
            cGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGtac                                                                                                                                    >  1:1882651/1‑148 (MQ=255)
            cGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGtac                                                                                                                                    <  2:1882651/148‑1 (MQ=255)
             ggTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGtacttc                                                                                                                                 <  1:2271502/150‑1 (MQ=255)
                   aTGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGtacttctacttc                                                                                                                           >  1:1430198/1‑150 (MQ=255)
                    tGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGCCAAACACCGGCCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCa                                                                                                                          <  2:2301308/150‑1 (MQ=255)
                    tGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCa                                                                                                                          <  2:149357/150‑1 (MQ=255)
                    tGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCa                                                                                                                          <  2:156261/150‑1 (MQ=255)
                    tGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCa                                                                                                                          <  2:897205/150‑1 (MQ=255)
                    tGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCa                                                                                                                          <  2:2324668/150‑1 (MQ=255)
                     gAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGACATTCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAg                                                                                                                         <  1:1843508/150‑1 (MQ=255)
                         cTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCa                                                                                                                     >  2:452937/1‑150 (MQ=255)
                         cTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCa                                                                                                                     >  2:1756228/1‑150 (MQ=255)
                          tCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCa                                                                                                                          <  1:1526057/144‑1 (MQ=255)
                          tCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCa                                                                                                                          >  2:1526057/1‑144 (MQ=255)
                                 tGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTg                                                                                                             >  2:291138/1‑150 (MQ=255)
                                 tGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTg                                                                                                             >  2:420944/1‑150 (MQ=255)
                                         ccTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGAc                                                                                                                                          >  2:1320293/1‑113 (MQ=255)
                                         ccTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGAc                                                                                                                                          <  1:1320293/113‑1 (MQ=255)
                                            gCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCAt                                                                                                  <  1:452937/150‑1 (MQ=255)
                                             cGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGCGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATg                                                                                                 <  1:291138/150‑1 (MQ=255)
                                             cGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATg                                                                                                 <  1:2050190/150‑1 (MQ=255)
                                                 ttACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTg                                                                                             <  2:442333/150‑1 (MQ=255)
                                                   aCGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGaa                                                                                           <  2:2107046/150‑1 (MQ=255)
                                                      gCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCt                                                                                        >  2:761213/1‑150 (MQ=255)
                                                         tGGCTCACCTGAGTGCTGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAg                                                                                     <  2:1216615/150‑1 (MQ=255)
                                                            cTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGGCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGCTTGAAGCTCAGGGa                                                                                  >  1:1257733/1‑150 (MQ=255)
                                                               aCCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCa                                                                               >  1:950364/1‑150 (MQ=255)
                                                                ccTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAg                                                                              <  1:1465991/150‑1 (MQ=255)
                                                                 cTGAGTGCGGAAAGGCCCTTGAAAGCCTCACGCGCATCGGGAGCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAgg                                                                             >  1:2300230/1‑150 (MQ=255)
                                                                    aGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGc                                                                          <  1:401177/150‑1 (MQ=255)
                                                                         ggAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGGAATTCCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTg                                                                     >  1:155672/1‑150 (MQ=255)
                                                                         ggAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTg                                                                     >  1:472936/1‑150 (MQ=255)
                                                                         ggAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTg                                                                     >  2:2061085/1‑150 (MQ=255)
                                                                            aaGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCgg                                                                  <  1:495404/150‑1 (MQ=255)
                                                                            aaGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCgg                                                                  <  1:1838438/150‑1 (MQ=255)
                                                                             aGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCg                                                                                                                                       <  1:90573/80‑1 (MQ=255)
                                                                             aGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCg                                                                                                                                       >  2:90573/1‑80 (MQ=255)
                                                                              ggCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTa                                                                <  1:1756228/150‑1 (MQ=255)
                                                                              ggCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTa                                                                <  1:420944/150‑1 (MQ=255)
                                                                               gCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTa                                                                <  1:351081/149‑1 (MQ=255)
                                                                               gCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTa                                                                >  2:351081/1‑149 (MQ=255)
                                                                               gCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAg                                                               >  1:171073/1‑150 (MQ=255)
                                                                               gCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAg                                                               >  2:1533324/1‑150 (MQ=255)
                                                                                 ccTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGc                                                                                                            >  1:515265/1‑103 (MQ=255)
                                                                                 ccTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGc                                                                                                            <  2:515265/103‑1 (MQ=255)
                                                                                 ccTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGaa                                                                                           <  2:1383695/120‑1 (MQ=255)
                                                                                 ccTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGaa                                                                                           >  1:1383695/1‑120 (MQ=255)
                                                                                 ccTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGc                                                             >  1:213803/1‑150 (MQ=255)
                                                                                   ttGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCa                                                           >  1:2161428/1‑150 (MQ=255)
                                                                                        aGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAgg                                                      <  1:1533324/150‑1 (MQ=255)
                                                                                              aaGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAg                                                <  2:950364/150‑1 (MQ=255)
                                                                                                      cGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGaa                                        <  2:1430198/150‑1 (MQ=255)
                                                                                                        ggATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCAt                                                                                                  <  1:1988303/90‑1 (MQ=255)
                                                                                                        ggATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCAt                                                                                                  >  2:1988303/1‑90 (MQ=255)
                                                                                                         gATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAgac                                     >  2:914550/1‑150 (MQ=255)
                                                                                                                             gTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCa                                                                               <  1:2441430/88‑1 (MQ=255)
                                                                                                                             gTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCa                                                                               >  2:2441430/1‑88 (MQ=255)
                                                                                                                             gTACGGGTCGTGCGCGATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGaa                                                                                           <  2:451615/76‑1 (MQ=25)
                                                                                                                             gTACGGGTCGTGCGCGATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGaa                                                                                           >  1:451615/1‑76 (MQ=25)
                                                                                                                             gTACGGGTCGTGCGCGATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGt                                            <  2:1142071/123‑1 (MQ=255)
                                                                                                                             gTACGGGTCGTGCGCGATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGt                                            >  1:1142071/1‑123 (MQ=255)
                                                                                                                              tACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTc                >  1:719070/1‑150 (MQ=255)
                                                                                                                               aCGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCa               <  1:761213/150‑1 (MQ=255)
                                                                                                                                cGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCAt              <  1:313803/150‑1 (MQ=255)
                                                                                                                                 gggTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCGCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATg             >  2:948577/1‑150 (MQ=255)
                                                                                                                                 gggTCGTGCGCGATCAGCGCACGACCCGTACTTTTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATg             >  2:2056198/1‑150 (MQ=255)
                                                                                                                                   gTCGTGCGCGATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGaa                                        <  2:2247636/121‑1 (MQ=255)
                                                                                                                                   gTCGTGCGCGATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGaa                                        >  1:2247636/1‑121 (MQ=255)
                                                                                                                                   gTCGTGCGCGATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGa           >  2:1117396/1‑150 (MQ=255)
                                                                                                                                    tCGTGCGCGATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGaa                                                                                           <  1:1919117/69‑1 (MQ=25)
                                                                                                                                    tCGTGCGCGATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGaa                                                                                           >  2:1919117/1‑69 (MQ=25)
                                                                                                                                      gTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATCg        >  1:395280/1‑150 (MQ=255)
                                                                                                                                       tGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATCGc       >  2:314183/1‑150 (MQ=255)
                                                                                                                                            tGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCa                                                                               >  1:1060344/1‑73 (MQ=255)
                                                                                                                                            tGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATCGCAACAg  >  1:2435564/1‑150 (MQ=255)
                                                                                                                                            tGATCGCGCACCACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCa                                                                               <  2:1060344/73‑1 (MQ=255)
                                                                                                                                               tCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAgac                                     >  1:1572422/1‑112 (MQ=255)
                                                                                                                                               tCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAgac                                     <  2:1572422/112‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATCGCAACAG  >  NC_002947/43317‑43606

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