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

AAGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATCGCAACAG  >  NC_002947/43313‑43606
                                                                                                                                                |                                                                                                                                                      
aaGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATcg                                                                                                                                                  >  1:1614757/1‑150 (MQ=255)
 aGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATcgc                                                                                                                                                 <  2:1187817/150‑1 (MQ=255)
 aGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGGGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGCAAAGTACGGGACGTCCGCTGATcgc                                                                                                                                                 >  1:587817/1‑150 (MQ=255)
   tCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCgcgc                                                                                                                                               >  1:952368/1‑150 (MQ=255)
     ggTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCGATCAGCGCAc                                                                                                                                             >  1:495147/1‑150 (MQ=255)
      gTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACg                                                                                                                                            >  2:452799/1‑150 (MQ=255)
      gTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACg                                                                                                                                            >  1:1437288/1‑150 (MQ=255)
       tttCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGa                                                                                                                                           <  1:883032/150‑1 (MQ=255)
        ttCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGAc                                                                                                                                          <  1:452799/150‑1 (MQ=255)
          cGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGAccc                                                                                                                                        <  1:1703408/150‑1 (MQ=255)
            tAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCCGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGt                                                                                                                                      <  2:422814/150‑1 (MQ=255)
                 ggTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGCAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTGGTGCGCTAATCGCGCACGACCCGtacttc                                                                                                                                 >  2:1219521/1‑150 (MQ=255)
                      cATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCGATCAGCGCACGACCCGtacttttctct                                                                                                                            <  2:495147/150‑7 (MQ=38)
                        tGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGa                                                                                                                                           >  1:650892/1‑133 (MQ=255)
                        tGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGa                                                                                                                                           <  2:650892/133‑1 (MQ=255)
                         gAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAg                                                                                                                         >  2:456381/1‑150 (MQ=255)
                         gAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGATGATAGCGCACGACCCGTACTTCTACTTCAg                                                                                                                         >  2:452074/1‑150 (MQ=255)
                          aGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGAc                                                                                                                                          >  2:142122/1‑132 (MQ=255)
                          aGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGAc                                                                                                                                          <  1:142122/132‑1 (MQ=255)
                          aGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAgg                                                                                                                        >  1:1504479/1‑150 (MQ=255)
                             cTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCa                                                                                                                     >  1:965459/1‑150 (MQ=255)
                             cTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCa                                                                                                                     >  1:1058991/1‑150 (MQ=255)
                               cAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCt                                                                                                                                                        >  1:852778/1‑114 (MQ=255)
                               cAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCt                                                                                                                                                        <  2:852778/114‑1 (MQ=255)
                                         ttGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTa                                                                                                         <  1:1467453/150‑1 (MQ=255)
                                         ttGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGAGCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTa                                                                                                         >  1:1590130/1‑150 (MQ=255)
                                              cTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGtacttc                                                                                                                                 >  2:603354/1‑121 (MQ=255)
                                              cTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGtacttc                                                                                                                                 <  1:603354/121‑1 (MQ=255)
                                              cTGCGGCTTACGGCATGGCTCACCTGAGTCAGGAAACGCCCGTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGGCATTTCCCTGCTTAGCTTc                                                                                                    >  1:1826726/1‑150 (MQ=255)
                                                gCGGCTTACGGCATGGCTCAGCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAGAGTACGGGTCGTGCGCTGATCGCGGACGACCCGTACTTCTACTTCACGTCATTTCCCTGCTTCGCTTCAt                                                                                                  <  2:1283629/150‑1 (MQ=255)
                                                gCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCAt                                                                                                  <  2:1058991/150‑1 (MQ=255)
                                                  ggCTTACGGCATGGCTCACCTGAGTGCGGAATGGCCCTTGAAAGGCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGTTCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGGTTAGCTTCATGa                                                                                                <  2:587817/150‑1 (MQ=255)
                                                  ggCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGa                                                                                                >  1:955584/1‑150 (MQ=255)
                                                   gCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGAt                                                                                               <  1:933987/150‑1 (MQ=255)
                                                    cTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGAtt                                                                                              >  2:1303563/1‑150 (MQ=255)
                                                    cTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGAtt                                                                                              <  2:1504479/150‑1 (MQ=255)
                                                       aCGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAACGTACGGGTCGTGCGCTGATCGCGCACGTCCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGaa                                                                                           <  2:1826726/150‑1 (MQ=255)
                                                        cGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAg                                                                                          <  2:952368/150‑1 (MQ=255)
                                                         ggCATGGCTCAGCTGAGTGCGAAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGTGCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGc                                                                                         >  2:24942/1‑150 (MQ=255)
                                                             tGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGCCTGACGCTCAg                                                                                     >  2:409295/1‑150 (MQ=255)
                                                             tGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAg                                                                                     >  1:226501/1‑150 (MQ=255)
                                                              ggCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAgg                                                                                    <  2:226501/150‑1 (MQ=255)
                                                                     cTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGAtt                                                                                              <  2:376362/133‑1 (MQ=255)
                                                                     cTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGAtt                                                                                              >  1:376362/1‑133 (MQ=255)
                                                                     cTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAgg                                                                             <  1:1715763/150‑1 (MQ=255)
                                                                            cGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTcc                                                                                                                >  1:129144/1‑108 (MQ=255)
                                                                            cGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTcc                                                                                                                <  2:129144/108‑1 (MQ=255)
                                                                            cGGAAAGGCCCTGGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCt                                                                      <  2:747221/150‑1 (MQ=255)
                                                                             ggAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTg                                                                     <  2:965459/150‑1 (MQ=255)
                                                                               aaaGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCg                                                                   >  2:1188499/1‑150 (MQ=255)
                                                                                 aGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGt                                                                 >  2:298969/1‑150 (MQ=255)
                                                                                  ggCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTa                                                                <  2:773449/150‑1 (MQ=255)
                                                                                   gCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAg                                                               <  2:1614757/150‑1 (MQ=255)
                                                                                       ttGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGATTCTGCGGTAGGCCa                                                           >  2:1076108/1‑150 (MQ=255)
                                                                                           aaGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAg                                                       >  2:1012053/1‑150 (MQ=255)
                                                                                           aaGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAg                                                       >  1:1706872/1‑150 (MQ=255)
                                                                                               cTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGcc                                                            >  2:1157588/1‑141 (MQ=255)
                                                                                               cTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGcc                                                            <  1:1157588/141‑1 (MQ=255)
                                                                                                  aaGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAg                                                <  2:1437288/150‑1 (MQ=255)
                                                                                                    gcgcATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTg                                                                                             <  2:969537/103‑1 (MQ=255)
                                                                                                    gcgcATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTg                                                                                             >  1:969537/1‑103 (MQ=255)
                                                                                                            ggATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAga                                      >  2:270310/1‑150 (MQ=255)
                                                                                                             gATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAgac                                     >  1:791339/1‑150 (MQ=255)
                                                                                                               tCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAgacga                                   <  1:697232/150‑1 (MQ=255)
                                                                                                                 gacgaGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGAcc                                 >  1:1191039/1‑150 (MQ=255)
                                                                                                                       tGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTccc                                                                                                               >  2:76007/1‑66 (MQ=255)
                                                                                                                       tGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTccc                                                                                                               <  1:76007/66‑1 (MQ=255)
                                                                                                                         agagaAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAg                                                               <  1:1323455/112‑1 (MQ=255)
                                                                                                                         agagaAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGATAg                                                               >  2:1323455/1‑112 (MQ=255)
                                                                                                                          gagaAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGgc                                                     >  2:377902/1‑121 (MQ=255)
                                                                                                                          gagaAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGgc                                                     <  1:377902/121‑1 (MQ=255)
                                                                                                                          gagaAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGtgtg                        <  2:1590130/150‑1 (MQ=255)
                                                                                                                           agaAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTgg                       >  2:1802610/1‑150 (MQ=255)
                                                                                                                            gaAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGa                      <  1:1802610/150‑1 (MQ=255)
                                                                                                                             aaaaGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAgg                                                                                    >  2:207061/1‑87 (MQ=255)
                                                                                                                             aaaaGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAgg                                                                                    <  1:207061/87‑1 (MQ=255)
                                                                                                                               aaGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTTATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGcc                   <  1:1806751/150‑1 (MQ=255)
                                                                                                                                 gTACGGGTCGTGCGCGATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCt                 >  2:349217/1‑150 (MQ=255)
                                                                                                                                  tACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAgtg                          <  1:219810/140‑1 (MQ=255)
                                                                                                                                  tACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAgtg                          >  2:219810/1‑140 (MQ=255)
                                                                                                                                    cGGGTCGTGCGCGATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGa                      <  1:1815854/142‑1 (MQ=255)
                                                                                                                                    cGGGTCGTGCGCGATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGa                      >  2:1815854/1‑142 (MQ=255)
                                                                                                                                     gggTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATg             <  1:1012053/150‑1 (MQ=255)
                                                                                                                                     gggTCGTGCGCGATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGc                                                                    <  2:414536/95‑1 (MQ=38)
                                                                                                                                     gggTCGTGCGCGATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGc                                                                    >  1:414536/1‑95 (MQ=38)
                                                                                                                                      ggTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATgg            >  1:1483241/1‑150 (MQ=255)
                                                                                                                                      ggTCGTGCGCGATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATgg            >  1:860888/1‑150 (MQ=255)
                                                                                                                                        tCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGAt          >  1:1863417/1‑150 (MQ=255)
                                                                                                                                           tGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATCGc       <  2:955584/150‑1 (MQ=255)
                                                                                                                                             cgcTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATCGCaa     >  1:634528/1‑150 (MQ=255)
                                                                                                                                              gcTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCAt              <  1:1938087/140‑1 (MQ=255)
                                                                                                                                              gcTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCAt              >  2:1938087/1‑140 (MQ=255)
                                                                                                                                                tGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATCGCAACAg  <  1:1188499/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
AAGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATCGCAACAG  >  NC_002947/43313‑43606

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: