Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 48,872 T→G 6.1% 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 6.1% 209.9 / 3.7 75D168A (GAC→GCC) czcA‑Ication efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (36/34);  new base G (2/3);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CAACGCTGCGACCAGGTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCTTCTGGCAGCTGTTC  >  NC_002947/48724‑49020
                                                                                                                                                    |                                                                                                                                                    
cAACGCTGCGACCAGGTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTc                                                                                                                                                     >  1:1600733/1‑150 (MQ=255)
 aaCGCTGCGACCAGGTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTcc                                                                                                                                                    >  1:185294/1‑150 (MQ=255)
   cGCTGCGACCAGGTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTg                                                                                                                                                  <  2:1334535/150‑1 (MQ=255)
     cTGCGACCAGGTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGa                                                                                                                                                >  2:1562402/1‑150 (MQ=255)
     cTGCGACCAGGTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGa                                                                                                                                                <  2:480426/150‑1 (MQ=255)
      tGCGACCAGGTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGAt                                                                                                                                               >  2:972301/1‑150 (MQ=255)
         gACCAGGTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGGACTGGATCat                                                                                                                                            <  2:14796/150‑2 (MQ=255)
           ccAGGTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCAcgc                                                                                                                                          >  1:1416991/1‑150 (MQ=255)
              ggTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGTGACGTTACGCAGCTGAGGCTTGATGATCCAGGACTGGATCatcaagc                                                                                                                                       <  1:305606/150‑7 (MQ=255)
               gtaaTTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAg                                                                                                                                      >  2:1835689/4‑150 (MQ=255)
                              cTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTa                                                                                                                       >  1:1462975/1‑150 (MQ=255)
                                  taggtagCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGc                                                                                                                   <  1:84973/150‑1 (MQ=255)
                                  taggtagCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGc                                                                                                                   <  1:1662666/150‑1 (MQ=255)
                                   aggtagCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGAGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCg                                                                                                                  <  1:312861/150‑1 (MQ=255)
                                       agCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTAcc                                                                                                              <  2:1462975/150‑1 (MQ=255)
                                          cAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTgg                                                                                                                                                 <  1:2022/112‑1 (MQ=255)
                                          cAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTgg                                                                                                                                                 >  2:2022/1‑112 (MQ=255)
                                          cAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTACTAGGTCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTc                                                                                                           >  1:1476140/1‑150 (MQ=255)
                                                cTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCttt                                                                                                     <  2:185294/150‑1 (MQ=255)
                                                cTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCAACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACTCGCATGTCGGTCGGGGTGTACGGCGTACCGTCCTCttt                                                                                                     <  2:220965/150‑1 (MQ=255)
                                                 tttGGAACCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGGACCTAAGGCATGATGCCCAAGTGCTGGGTATCGAACAGGAACGTAGGAGCGTCCGTCGTACCATCCTCTTTg                                                                                                    >  2:365530/1‑150 (MQ=255)
                                                           gcgcAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCAt                                                                                          <  1:972301/150‑1 (MQ=255)
                                                                aCCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCg                                                                                     <  1:1835689/150‑1 (MQ=255)
                                                                 ccAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGGACTGGATCatcaa                                                                                                                                         <  2:445401/97‑5 (MQ=255)
                                                                 ccAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGGACTGGATCatcaa                                                                                                                                         >  1:445401/1‑93 (MQ=255)
                                                                   aGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCt                                                                                  >  1:1264097/1‑150 (MQ=255)
                                                                    ggAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCtt                                                                                 >  1:678425/1‑150 (MQ=255)
                                                                          gCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCgg                                                                                    >  1:1447327/1‑141 (MQ=255)
                                                                          gCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCgg                                                                                    <  2:1447327/141‑1 (MQ=255)
                                                                           cTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGt                                                                          <  1:365530/150‑1 (MQ=255)
                                                                              ggCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCca                                                                       <  1:707247/150‑1 (MQ=255)
                                                                                 aTAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACa                                                                            >  2:558670/1‑142 (MQ=255)
                                                                                 aTAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACa                                                                            <  1:558670/142‑1 (MQ=255)
                                                                                 aTAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCACCATCTTCGGCTTCAACAGTCCACAg                                                                    >  1:1366547/1‑150 (MQ=255)
                                                                                  tAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCcaca                                                                     >  1:1330171/1‑148 (MQ=255)
                                                                                  tAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCcaca                                                                     <  2:1330171/148‑1 (MQ=255)
                                                                                   aaCCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCggg                                                                                                                            <  2:224928/92‑1 (MQ=255)
                                                                                   aaCCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCggg                                                                                                                            >  1:224928/1‑92 (MQ=255)
                                                                                     ccgccgATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCcaca                                                                     >  1:552264/1‑145 (MQ=255)
                                                                                     ccgccgATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCcaca                                                                     <  2:552264/145‑1 (MQ=255)
                                                                                          gATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGAtctc                                                           >  1:64234/1‑150 (MQ=255)
                                                                                            tGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGc                                                         >  2:1264992/1‑150 (MQ=255)
                                                                                                aTTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCaa                                                                              <  2:741471/125‑1 (MQ=255)
                                                                                                aTTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCaa                                                                              >  1:741471/1‑125 (MQ=255)
                                                                                                    ccatcGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCg                                                 <  1:1754280/147‑1 (MQ=255)
                                                                                                     tctcGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTc                                                                         <  2:487060/125‑1 (MQ=255)
                                                                                                     tctcGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTc                                                                         >  1:487060/1‑125 (MQ=255)
                                                                                                     tctcGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCgg                                                >  2:1446495/1‑150 (MQ=255)
                                                                                                       tcGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTa                                              <  1:1446495/150‑1 (MQ=255)
                                                                                                       tcGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTa                                              >  1:54430/1‑150 (MQ=255)
                                                                                                       tcGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTa                                              >  1:212628/1‑150 (MQ=255)
                                                                                                            cacaCCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATAc                                         <  2:64234/150‑1 (MQ=255)
                                                                                                             acacCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCTCGCCATCTTCGGCTTCACCAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATAcc                                        >  1:147490/1‑150 (MQ=255)
                                                                                                                cccGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCCGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGGCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGa                                     <  2:1416991/150‑1 (MQ=255)
                                                                                                                 ccGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCttt                                                                                                     >  1:1094006/1‑85 (MQ=255)
                                                                                                                 ccGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCttt                                                                                                     <  2:1094006/85‑1 (MQ=255)
                                                                                                                  cGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGAcc                                   >  2:1614369/1‑150 (MQ=255)
                                                                                                                  cGGGACGATACGCAGCTGAGGCTTGATGATCCAGGCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACAGTCCTCTTTCACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCAGCTAGATACCGGAcc                                   >  2:621924/1‑150 (MQ=255)
                                                                                                                     gACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCAt                                <  2:212628/150‑1 (MQ=255)
                                                                                                                         ttACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTATCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACg                            >  1:764625/1‑150 (MQ=255)
                                                                                                                         ttACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTCCAACAGTCCACAGGAAGATCTCGCAGAGGCCGGTAGATACCGGACCCATGACg                            <  1:541430/150‑1 (MQ=255)
                                                                                                                           aCGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGc                          >  2:223727/1‑150 (MQ=255)
                                                                                                                             gCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCca                                                                       >  1:1652457/1‑103 (MQ=255)
                                                                                                                             gCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCca                                                                       <  2:1652457/103‑1 (MQ=255)
                                                                                                                             gCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCct                        <  1:1264992/150‑1 (MQ=255)
                                                                                                                                gCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTa                     <  2:54430/150‑1 (MQ=255)
                                                                                                                                gCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTa                     >  1:1030078/1‑150 (MQ=255)
                                                                                                                                gCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTa                     <  2:678425/150‑1 (MQ=255)
                                                                                                                                      gCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCtt               >  1:787849/1‑150 (MQ=255)
                                                                                                                                        ttGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCTTCt             <  1:1562402/150‑1 (MQ=255)
                                                                                                                                           atgatCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCtt                                                                                 >  2:166206/1‑79 (MQ=255)
                                                                                                                                           atgatCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCtt                                                                                 <  1:166206/79‑1 (MQ=255)
                                                                                                                                             gatCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCTTCTGGCAg        <  1:1352846/150‑1 (MQ=255)
                                                                                                                                               tCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCTTCTGGCAGCt      >  1:871814/1‑150 (MQ=255)
                                                                                                                                                   gTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCTTCTGGCAGCTGTTc  >  1:1576254/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAACGCTGCGACCAGGTCATTGAGTGTCAGCTTGTAGGTAGCCAGGCGCTTTGGATCCGGCGCAACCAGGAACTGCTTGGCATAACCGCCGATGGTATTGATCTCGGCCACACCCGGGACGTTACGCAGCTGAGGCTTGATGATCCAGTCCTGGATCACGCGCAGGTCGGTCGGGGTGTACGGCGTACCGTCCTCTTTGACCGCGCCATCTTCGGCTTCAACAGTCCACAGGAAGATCTCGCCGAGGCCGGTAGATACCGGACCCATGACGGCCTCTACACCTTCTGGCAGCTGTTC  >  NC_002947/48724‑49020

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: