Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 43,457 Δ1 bp 13.3% intergenic (+23/+70) 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,4570T.13.3% 242.3 / 31.8 83intergenic (+23/+70)PP_0040/cadA‑Ihypothetical protein/cadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base T (33/39);  new base . (6/5);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CAATCAAGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATCGCAACAG  >  NC_002947/43308‑43606
                                                                                                                                                     |                                                                                                                                                      
cAATCAAGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTAAGAGAAAAGTACGGGTCGTGCGCt                                                                                                                                                        <  2:903263/150‑1 (MQ=255)
 aaTCAAGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTg                                                                                                                                                       >  1:156261/1‑150 (MQ=255)
 aaTCAAGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTg                                                                                                                                                       >  1:149357/1‑150 (MQ=255)
  aTCAAGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGa                                                                                                                                                      <  2:1750766/150‑1 (MQ=255)
   tCAAGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCCGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGAt                                                                                                                                                     <  2:1854680/150‑1 (MQ=255)
         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                                                                                                                                    <  2:1882651/148‑1 (MQ=255)
                     cGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGtac                                                                                                                                    >  1:1882651/1‑148 (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:2324668/150‑1 (MQ=255)
                             tGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCa                                                                                                                          <  2:897205/150‑1 (MQ=255)
                             tGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCa                                                                                                                          <  2:149357/150‑1 (MQ=255)
                             tGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCa                                                                                                                          <  2:156261/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                                                                                                                          >  2:1526057/1‑144 (MQ=255)
                                   tCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCa                                                                                                                          <  1:1526057/144‑1 (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                                                                                                                                       >  2:90573/1‑80 (MQ=255)
                                                                                      aGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCg                                                                                                                                       <  1:90573/80‑1 (MQ=255)
                                                                                       ggCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTa                                                                <  1:420944/150‑1 (MQ=255)
                                                                                       ggCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTa                                                                <  1:1756228/150‑1 (MQ=255)
                                                                                        gCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTa                                                                <  1:351081/149‑1 (MQ=255)
                                                                                        gCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTa                                                                >  2:351081/1‑149 (MQ=255)
                                                                                        gCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAg                                                               >  2:1533324/1‑150 (MQ=255)
                                                                                        gCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAg                                                               >  1:171073/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                                                                                                  >  2:1988303/1‑90 (MQ=255)
                                                                                                                 ggATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCAt                                                                                                  <  1:1988303/90‑1 (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)
                                                                                                                                                     |                                                                                                                                                      
CAATCAAGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATCGCAACAG  >  NC_002947/43308‑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: