Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 30,651 C→T 6.2% A74V (GCG→GTG)  PP_0026 → cobalt/cadmium/zinc exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294730,6510CT 6.2% 186.0 / 2.3 64A74V (GCG→GTG) PP_0026cobalt/cadmium/zinc exporter
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base T (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

GGCTTGACTGGCAGCTTCATGATTGCAGAAGTGATCGGTGCGTGGATTACCGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCTTGTATGAGGCCTACCAGCGCTTCTTCATGCCGGCGGAGATCGCTACCGGAGCGATGA  >  NC_002947/30503‑30800
                                                                                                                                                    |                                                                                                                                                     
ggCTTGACTGGCAGCTTCATGATTGCAGAAGTGATCGGTGCGTGGATTACCGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCcggcg                                                                                                                                                      >  2:988233/1‑150 (MQ=255)
  cTTGACTGGCAGCTTCATGATTGCAGAAGTGATCGGTGCGTGGATTACCGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGa                                                                                                                                                    >  2:1865321/1‑150 (MQ=255)
        tGGCAGCTTCATGATTGCAGAAGTGATCGGTGCGTGGATTACCGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGaa                                                                                                                                              <  1:1576843/150‑1 (MQ=255)
        tGGCAGATTCATGATTGCAGAAGTGATCGGTGCGTGGATTACCGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGaa                                                                                                                                              >  1:797269/1‑150 (MQ=255)
          gCAGCTTCATGATTGCAGAAGTGATCGGTGCGTGGATTACCGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGaaaa                                                                                                                                            <  2:1649315/150‑1 (MQ=255)
           cAGCTTCATGATTGCAGAAGTGATCGGTGCGTGGATTACCGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAg                                                                                                                                           >  1:767362/1‑150 (MQ=255)
                    gATTGCAGAAGTGATCGGTGCGTGGATTACCGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCgg                                                                                                                                  >  1:963319/1‑150 (MQ=255)
                     aTTGCAGAAGTGATCGGTGCGTGGATTACCGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAGCCTTCGGc                                                                                                                                 >  1:1396678/1‑150 (MQ=255)
                          aGAAGTGATCGGTGCGTGGATTACCGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATgc                                                                                                                            >  2:1806565/1‑150 (MQ=255)
                            aaGTGATCGGTGCGTGGATTACCGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATgcgc                                                                                                                          >  2:1930914/1‑150 (MQ=255)
                                     gTGCGTGGATTACCGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGaaa                                                                                                                 >  1:562869/1‑150 (MQ=255)
                                           ggATTACCGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTgg                                                                                                           <  1:1939807/150‑1 (MQ=255)
                                           ggATTACCGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTgg                                                                                                           <  2:331184/150‑1 (MQ=255)
                                               tACCGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTc                                                                                                       <  1:277745/150‑1 (MQ=255)
                                                  cGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCt                                                                                                                                <  1:1439042/122‑1 (MQ=255)
                                                  cGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCt                                                                                                                                >  2:1439042/1‑122 (MQ=255)
                                                  cGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGAc                                                                                                    >  1:667296/1‑150 (MQ=255)
                                                  cGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGAc                                                                                                    <  1:271273/150‑1 (MQ=255)
                                                     cAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATATCCCTGATTGCACTTTAGAGAGTCAATCTGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGtt                                                                                                 <  2:1674592/150‑1 (MQ=255)
                                                         cTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAAc                                                                                             >  1:1081984/1‑150 (MQ=255)
                                                          tGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAAcc                                                                                                                                       >  1:1386485/1‑107 (MQ=255)
                                                          tGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAAcc                                                                                                                                       <  2:1386485/107‑1 (MQ=255)
                                                          tGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACg                                                                                            >  1:742138/1‑150 (MQ=255)
                                                             cgcTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCg                                                                                         <  2:767362/150‑1 (MQ=255)
                                                               ctgctgTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATGTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGtg                                                                                       <  2:840482/150‑1 (MQ=255)
                                                               ctgctgTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGtg                                                                                       <  2:667296/150‑1 (MQ=255)
                                                                tgctgTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTc                                                                                                       <  1:323608/133‑1 (MQ=255)
                                                                tgctgTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTc                                                                                                       >  2:323608/1‑133 (MQ=255)
                                                                         aCGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTg                                                                                                            >  2:1545815/1‑119 (MQ=255)
                                                                         aCGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTg                                                                                                            <  1:1545815/119‑1 (MQ=255)
                                                                           gCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGGTTCGCGGGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCtg                                                                           >  1:119388/1‑150 (MQ=255)
                                                                           gCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCtg                                                                           >  2:1482130/1‑150 (MQ=255)
                                                                             aTCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCtggt                                                                         >  1:801148/1‑150 (MQ=255)
                                                                                tcACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGTGGATCAGAAAAGAACCTTCGGCt                                                                                                                                <  1:1443870/92‑1 (MQ=255)
                                                                                tcACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGTGGATCAGAAAAGAACCTTCGGCt                                                                                                                                <  1:1444276/92‑1 (MQ=255)
                                                                                tcACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGTGGATCAGAAAAGAACCTTCGGCt                                                                                                                                >  2:1443870/1‑92 (MQ=255)
                                                                                tcACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGTGGATCAGAAAAGAACCTTCGGCt                                                                                                                                >  2:1444276/1‑92 (MQ=255)
                                                                                tcACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGcc                                                                                                                       >  1:175117/1‑101 (MQ=255)
                                                                                tcACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGcc                                                                                                                       <  2:175117/101‑1 (MQ=255)
                                                                                  acaTGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCcggcgg                                                                                                                                                     >  1:1804032/1‑69 (MQ=255)
                                                                                  acaTGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCcggcgg                                                                                                                                                     <  2:1804032/69‑1 (MQ=255)
                                                                                  acaTGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGa                                                                    >  2:858516/1‑150 (MQ=255)
                                                                                  acaTGTTCACAGATAATGCCGTCTTGGAGATGTCCCTGATTGCACTTCAGCTGGCCAAGCGGTCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGa                                                                    <  2:119388/150‑1 (MQ=255)
                                                                                   caTGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTgg                                                                                                           <  1:1374781/110‑1 (MQ=255)
                                                                                   caTGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTgg                                                                                                           >  2:1374781/1‑110 (MQ=255)
                                                                                        tCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACa                                                              >  1:1078754/1‑150 (MQ=255)
                                                                                           caGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAACTATTGTTTGCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCt                                                           <  1:1806565/150‑1 (MQ=255)
                                                                                            aGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAGCCTTCGGCTATGCCCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCtt                                                          <  2:1396678/150‑1 (MQ=255)
                                                                                             gATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCTTg                                                         <  1:263168/150‑1 (MQ=255)
                                                                                              aTACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGa                                                                                                     <  1:223440/105‑1 (MQ=255)
                                                                                              aTACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGa                                                                                                     >  2:223440/1‑105 (MQ=255)
                                                                                                aCTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCTTGTAt                                                      <  2:801148/150‑1 (MQ=255)
                                                                                                  tGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCt                                                           >  1:555602/1‑143 (MQ=255)
                                                                                                  tGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCt                                                           <  2:555602/143‑1 (MQ=255)
                                                                                                  tGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCTTGTATGa                                                    <  2:797269/150‑1 (MQ=255)
                                                                                                   gccgccTTGGCGATCTCACTGATTTGACTGCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCCTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCTTGTATGAg                                                   >  1:1014002/1‑150 (MQ=255)
                                                                                                    ccgccTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCTTGTATGAgg                                                  <  1:1865321/150‑1 (MQ=255)
                                                                                                       ccTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGAGCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCTTGTATGAGGCCt                                               >  1:1563692/1‑150 (MQ=255)
                                                                                                         ttGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCTTGTATGAGGCCTAc                                             <  2:1336155/150‑1 (MQ=255)
                                                                                                                    aCTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCTTGTATGAGGCCTACCAGCGcttctt                                  >  2:742265/1‑150 (MQ=255)
                                                                                                                     cTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCTTGTATGAGGCCTACCAGCGCttcttc                                 <  1:742265/150‑1 (MQ=255)
                                                                                                                         ttGCACTTCAGATAGCCAAGCGGCCGGCGGATGGTAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCTTGTATGAGGCCTACCAGCGCTTCTTCCTGc                             >  1:1910981/1‑150 (MQ=255)
                                                                                                                               ttCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCTTGTATGAGGCCTACCAGCGCTTCTTCATGCcggcgg                       <  1:21328/150‑1 (MQ=255)
                                                                                                                                tCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCTTGTATGAGGCCTACCAGCGCTTCTTCATGCCGGCGga                      <  2:1081984/150‑1 (MQ=255)
                                                                                                                                  aGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTTGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCTTGTATGAGGCCTACCAGCGCTTCTTCATGCCGGCGgaga                    <  1:436765/150‑1 (MQ=255)
                                                                                                                                   gATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTGTGCACGCCTGGAAATTCGGGCTTCGACGTTCAACGCCGTCCTGCTCTTCCTGGTTGCGATGTACATCTTGTATGAGGCCTACCAGCGCTTCTTCATGACGGCGGAGAt                   >  1:137104/1‑150 (MQ=255)
                                                                                                                                        ccAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGACTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCTTGTATGAGGCCTACCAGCGCTTCTTCATGCCGGCGGAGATCGCTa              >  2:1403826/1‑150 (MQ=255)
                                                                                                                                               gccggcGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCTTGTATGAGGCCTACAAGCGCTTCTTCATGCCGGCGGAGATCGCTACCGGAGc       <  2:963319/150‑1 (MQ=255)
                                                                                                                                                    cggATCAGAAAAGAACCTTCGTCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCTTGTATGAGGCCTACCAGCGCTTCTTCATGCCGGCGGAGATCGCTACCGGAGCgatga  >  2:1609038/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCTTGACTGGCAGCTTCATGATTGCAGAAGTGATCGGTGCGTGGATTACCGGCAGCCTGGCGCTGCTGTCGGACGCATCTCACATGTTCACAGATACTGCCGCCTTGGCGATCTCACTGATTGCACTTCAGATAGCCAAGCGGCCGGCGGATCAGAAAAGAACCTTCGGCTATGCGCGCCTGGAAATTCTGGCTTCGACGTTCAACGCCGTGCTGCTCTTCCTGGTGGCGATGTACATCTTGTATGAGGCCTACCAGCGCTTCTTCATGCCGGCGGAGATCGCTACCGGAGCGATGA  >  NC_002947/30503‑30800

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: ‑