Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 48,872 T→G 8.4% D168A (GAC→GCC) ‡ czcA‑I ← cation efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294748,8720TG 8.4% 268.5 / 15.4 103D168A (GAC→GCC) ‡czcA‑Ication efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (55/39);  new base G (5/4);  total (60/43)
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.91e-01

Only 100 of 103 total aligned reads displayed.

AACGCTGCGACCAGGTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCTTCTGGCAGCTGT  >  NC_002947/48725‑49018
                                                                                                                                                   |                                                                                                                                                  
aaCGCTGCGACCAGGTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTcc                                                                                                                                                  >  2:2155912/1‑150 (MQ=255)
    cTGCGACCAGGTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGa                                                                                                                                              <  1:981594/150‑1 (MQ=255)
      gCGACCAGGTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATc                                                                                                                                            >  2:686655/1‑150 (MQ=255)
        gACCAGGTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCAc                                                                                                                                          >  2:438328/1‑150 (MQ=255)
         aCCAGGTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCAcg                                                                                                                                         >  1:819510/1‑150 (MQ=255)
            aGGTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGGACTGGATCatcaag                                                                                                                                      <  2:1190249/150‑6 (MQ=255)
             ggTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCa                                                                                                                                     >  1:1317725/1‑150 (MQ=255)
                      gtgtCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGgtcggtcg                                                                                                                            >  2:1476893/1‑150 (MQ=255)
                       tgtCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGtcggtcgg                                                                                                                           <  2:1829426/150‑1 (MQ=255)
                        gtCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTCAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCggg                                                                                                                          <  2:1749130/150‑1 (MQ=255)
                        gtCAACTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACTCCCGGGACGTTACGCAGCTGAGGCTTTATGATCCAGTCCTGGATCACGCGCAGGTCGGTCggg                                                                                                                          <  1:337050/150‑1 (MQ=255)
                              ttGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTAc                                                                                                                    <  2:1317725/150‑1 (MQ=255)
                                          aGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTcc                                                                                                        >  2:853129/1‑150 (MQ=255)
                                           ggCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCct                                                                                                       >  1:1278266/1‑150 (MQ=255)
                                                  tGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGAc                                                                                                <  2:97120/150‑1 (MQ=255)
                                                      tCCGGCGCAACCAGGAACTGCTTGGGATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACcgcg                                                                                            >  1:1530680/1‑150 (MQ=255)
                                                      tCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCg                                                                                                           <  2:2242153/135‑1 (MQ=255)
                                                      tCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCg                                                                                                           >  1:2242153/1‑135 (MQ=255)
                                                      tCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGATGAGGCTTGATGATCCAGTCCTg                                                                                                                                                >  2:896212/1‑98 (MQ=255)
                                                      tCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGATGAGGCTTGATGATCCAGTCCTg                                                                                                                                                <  1:896212/98‑1 (MQ=255)
                                                        cGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGcc                                                                                          >  2:962269/1‑150 (MQ=255)
                                                         ggCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCa                                                                                         >  1:1614424/1‑150 (MQ=255)
                                                          gcgcAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACCCCCGGGACCTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATGACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCAt                                                                                        >  2:1552209/1‑150 (MQ=255)
                                                          gcgcAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCggg                                                                                                                          <  2:1040456/116‑1 (MQ=255)
                                                          gcgcAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCggg                                                                                                                          >  1:1040456/1‑116 (MQ=255)
                                                          gcgcAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCAt                                                                                        >  2:324987/1‑150 (MQ=255)
                                                           cgcAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATc                                                                                       <  1:438328/150‑1 (MQ=255)
                                                            gcAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCt                                                                                      <  1:1476893/150‑1 (MQ=255)
                                                                  aGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGGACTGGATCa                                                                                                                                           >  1:2288136/1‑91 (MQ=255)
                                                                  aGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGGACTGGATCa                                                                                                                                           <  2:2288136/91‑1 (MQ=255)
                                                                    gAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTc                                                                              >  2:347471/1‑150 (MQ=255)
                                                                           ttGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTc                                                                       <  1:969848/150‑1 (MQ=255)
                                                                           ttGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGGACTGGATCatc                                                                                                                                         <  1:2033957/84‑3 (MQ=255)
                                                                           ttGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGGACTGGATCatc                                                                                                                                         >  2:2033957/1‑82 (MQ=255)
                                                                              gCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTAcc                                                                                                            >  1:79529/1‑110 (MQ=255)
                                                                              gCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTAcc                                                                                                            <  2:79529/110‑1 (MQ=255)
                                                                               cATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTGGGGGTGTACGGCATACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCcaca                                                                   >  2:2064487/1‑150 (MQ=255)
                                                                                 tAACCACCGATGGTATTGATCTAGGCCACACCCGGGACGTTACGAAGATGAGGTTTGATGATCCAGTGCTGGATCACGCGAAGGTCGGTCGGGGTGTACGGGGTACCGTCCTCGTTGACAGCGCCATTTTCGGCTTCAACAGTCGACAgg                                                                 >  2:2276879/1‑150 (MQ=255)
                                                                                   aCCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCggg                                                                                                                          <  1:1389240/91‑1 (MQ=255)
                                                                                   aCCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCggg                                                                                                                          >  2:1389240/1‑91 (MQ=255)
                                                                                    ccgccgATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTAc                                                                                                             >  2:625686/1‑103 (MQ=255)
                                                                                    ccgccgATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTAc                                                                                                             <  1:625686/103‑1 (MQ=255)
                                                                                     cgccgATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTCAAAAAGCCGGTTCTCTTTGACCGCGCCATCTTCGGCTTCAAAAGTCCACAGAAAGa                                                             >  1:1955376/1‑150 (MQ=21)
                                                                                     cgccgATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCcaca                                                                   >  1:907804/1‑144 (MQ=255)
                                                                                     cgccgATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCcaca                                                                   <  2:907804/144‑1 (MQ=255)
                                                                                      gccgATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGAt                                                            >  1:1515912/1‑150 (MQ=255)
                                                                                        cgATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAgg                                                                 >  1:1673963/1‑143 (MQ=255)
                                                                                        cgATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAgg                                                                 <  2:1673963/143‑1 (MQ=255)
                                                                                         gATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGAc                                                                                                >  2:281622/1‑111 (MQ=255)
                                                                                         gATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGAc                                                                                                <  1:281622/111‑1 (MQ=255)
                                                                                            tgTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGcc                                                      <  1:1843912/149‑1 (MQ=255)
                                                                                              ttctGATCTCGGCCACACCCGGGACGTTATGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGCGGTGTACGGCGTATCATACATAAACACCTTCTTACCGTTCGCGTCAACAGTCCACAGGAAGATCTCGCCGa                                                    <  1:1622844/147‑1 (MQ=12)
                                                                                                ttGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCggg                                                                                                                          >  2:1849054/1‑78 (MQ=255)
                                                                                                ttGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCggg                                                                                                                          <  1:1849054/78‑1 (MQ=255)
                                                                                                 tGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGACCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGc                                                 <  1:347471/150‑1 (MQ=255)
                                                                                                   aTCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAgagat                                               >  2:1622844/1‑146 (MQ=255)
                                                                                                   aTCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCg                                               <  1:587926/150‑1 (MQ=255)
                                                                                                   aTCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCg                                               >  1:747311/1‑150 (MQ=255)
                                                                                                    tctcGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCct                                                                                                       <  2:742848/93‑1 (MQ=255)
                                                                                                    tctcGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCct                                                                                                       >  1:742848/1‑93 (MQ=255)
                                                                                                    tctcGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCgg                                              >  2:2210754/1‑150 (MQ=255)
                                                                                                     ctcGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGt                                             <  2:468713/150‑1 (MQ=255)
                                                                                                      tcGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTa                                            >  1:1255510/1‑150 (MQ=255)
                                                                                                      tcGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTa                                            >  2:912494/1‑150 (MQ=255)
                                                                                                         gCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAg                                                                  >  2:2066635/1‑125 (MQ=255)
                                                                                                         gCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAg                                                                  <  1:2066635/125‑1 (MQ=255)
                                                                                                          ccACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGgt                                                                                                                        >  1:1022002/1‑70 (MQ=255)
                                                                                                          ccACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGgt                                                                                                                        <  2:1022002/70‑1 (MQ=255)
                                                                                                           cacaCCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGc                                                                                                                 >  2:1741481/1‑76 (MQ=255)
                                                                                                           cacaCCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGc                                                                                                                 <  1:1741481/76‑1 (MQ=255)
                                                                                                           cacaCCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATAc                                       <  1:2155912/150‑1 (MQ=255)
                                                                                                               cccGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGa                                   <  1:324987/150‑1 (MQ=255)
                                                                                                               cccGGGACGTTACGCAGCTGAGGCTTGATGATCCAGGCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCATACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGa                                   <  1:2064487/150‑1 (MQ=255)
                                                                                                                ccGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGAc                                  >  1:1094950/1‑150 (MQ=255)
                                                                                                                ccGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGAc                                  >  2:1097087/1‑150 (MQ=255)
                                                                                                                  gggACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGc                                                                                                                 >  1:2386591/1‑69 (MQ=255)
                                                                                                                  gggACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGc                                                                                                                 <  2:2386591/69‑1 (MQ=255)
                                                                                                                    gACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCAt                              >  1:236604/1‑150 (MQ=255)
                                                                                                                      cGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGa                            <  2:1530680/150‑1 (MQ=255)
                                                                                                                        ttACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACg                          >  1:90201/1‑150 (MQ=255)
                                                                                                                            gCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGAcc                                                                                               <  1:977377/77‑1 (MQ=255)
                                                                                                                            gCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCATCAAGATGCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCct                      >  2:933694/1‑150 (MQ=255)
                                                                                                                                cTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTac                  >  1:979095/1‑150 (MQ=255)
                                                                                                                                   aGGCTTGATGATCCAGTCCTGGATCACACGAAGGACGGTGGGGGTGTACGGCGTACCGTCCTCTTTTAGCACAGCATCTTCGACTTCTACAGTCCATAGGAAGATCTCGCCGAGGCCGGTAGAAACCGAACCCAAGACGGTCTGTACAcc               >  1:1976692/1‑150 (MQ=255)
                                                                                                                                     gCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGGGCCCGGTAGATACCGGACCCATGACGGCCTCTACACCtt             >  1:127003/1‑150 (MQ=255)
                                                                                                                                     gCTTGAGGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCtt             <  2:819510/150‑1 (MQ=255)
                                                                                                                                       ttGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCTTCt           >  1:1715055/1‑150 (MQ=255)
                                                                                                                                        ggatgctCCAGTCCTGGATCGTGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCTTCTg          <  2:1715055/149‑1 (MQ=255)
                                                                                                                                         gatgatCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCTTCTgg         <  1:853129/150‑1 (MQ=255)
                                                                                                                                         gatgatCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCTTCTgg         >  2:2277853/1‑150 (MQ=255)
                                                                                                                                          gtgatCCAGGACTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCTTCTGGc        >  2:1444634/2‑150 (MQ=255)
                                                                                                                                          gtgatCCAGGACTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCTTCTGGc        >  2:1526349/2‑150 (MQ=255)
                                                                                                                                          gtgatCCAGGACTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACACATGACGGCCTCTAAACctgctgtc        >  2:1978725/2‑148 (MQ=255)
                                                                                                                                            gatCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTCCGGCGTCCCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTCGCTACCGGACCACTGCCGGCCTCTATACCTGCGGGctg      <  2:1955376/150‑3 (MQ=255)
                                                                                                                                            gatCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACTGTCCACAg                                                                  >  1:532278/1‑90 (MQ=255)
                                                                                                                                            gatCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAg                                                                  <  2:532278/90‑1 (MQ=255)
                                                                                                                                             atCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCTTCTGGCAGc     >  2:966297/1‑150 (MQ=255)
                                                                                                                                              tCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCTTCTGGCAGCt    >  2:2117361/1‑150 (MQ=255)
                                                                                                                                               ccAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCTTCTGGCAGCTg   >  1:2243228/1‑150 (MQ=255)
                                                                                                                                                cAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCCCGCCGAGGCCGGTAGATACCGGACCCAATACACCCTCTACACCTTCTGGCAGCTGt  <  1:933694/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACGCTGCGACCAGGTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCTTCTGGCAGCTGT  >  NC_002947/48725‑49018

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