Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,213,020 T→A 8.2% M15K (ATG→AAG)  mexC → multidrug RND transporter membrane fusion protein MexC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,213,0200TA 8.2% 210.8 / 6.3 71M15K (ATG→AAG) mexCmultidrug RND transporter membrane fusion protein MexC
Reads supporting (aligned to +/- strand):  ref base T (33/32);  new base A (4/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGAAGGGGTCATTAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTC  >  NC_002947/3212872‑3213165
                                                                                                                                                    |                                                                                                                                                 
tCTGAAGGGGTCATTAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATg                                                                                                                                                  >  1:327990/1‑150 (MQ=255)
  tGAAGGGGTCATTAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGc                                                                                                                                                <  2:1382899/150‑1 (MQ=255)
          tCATTAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCt                                                                                                                                        >  2:394531/1‑150 (MQ=255)
             ttAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCAATAAGGGTGGGCGCATTGGACCAGGCCATCGCCTTgg                                                                                                                                     >  2:36594/1‑150 (MQ=255)
                  agTCGTTTATGATTGTCATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCAAGGCGATCGCCATGGccaat                                                                                                                                <  1:412684/150‑4 (MQ=255)
                   ttCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGgt                                                                                                                               <  1:555039/149‑1 (MQ=255)
                    tCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGgtg                                                                                                                              >  1:675385/1‑150 (MQ=255)
                    acTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGgtg                                                                                                                              >  2:979645/2‑150 (MQ=255)
                          aTGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGcc                                                                                                                        >  2:721737/1‑150 (MQ=255)
                          aTGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGcc                                                                                                                        <  1:39695/150‑1 (MQ=255)
                                   aTGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcaga                                                                                                               >  1:20547/1‑150 (MQ=255)
                                     gagaCAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATg                                                                                                                                                  <  1:1247034/113‑1 (MQ=255)
                                     gagaCAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATg                                                                                                                                                  >  2:1247034/1‑113 (MQ=255)
                                      agaCAATATTGAGTCATGTGGGTTTGAAAAGTAGCATCCGCGCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAGTACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagca                                                                                                            <  2:449653/150‑1 (MQ=255)
                                      agaCAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagca                                                                                                            <  1:378239/150‑1 (MQ=255)
                                         cAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGa                                                                                                         <  1:973149/150‑1 (MQ=255)
                                         cAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGa                                                                                                         >  1:1282439/1‑150 (MQ=255)
                                          aaTATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTCTTCAATCAAGGTTCGGGTAGGTGAGGAGCAAGACGGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATGGCCTTGGCCGGGTGTGAGCCCGCAGCATAGCAGGAc                                                                                                        >  2:340839/1‑150 (MQ=255)
                                              ttGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTg                                                                                                    >  1:1353160/1‑150 (MQ=255)
                                                  cTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGc                                                                                                >  1:244025/1‑150 (MQ=255)
                                                   tCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGcc                                                                                               >  1:158290/1‑150 (MQ=255)
                                                     aTGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCaa                                                                                             >  2:211235/1‑150 (MQ=255)
                                                           gTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAg                                                                                       <  2:20547/150‑1 (MQ=255)
                                                           gTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAg                                                                                       >  1:75252/1‑150 (MQ=255)
                                                                   aGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCgg                                                                               <  1:809416/150‑1 (MQ=255)
                                                                             gCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGt                                                                     >  2:869368/1‑150 (MQ=255)
                                                                             gCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGt                                                                     <  1:869368/150‑1 (MQ=255)
                                                                                ccGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGAc                                                                  >  2:1362612/1‑150 (MQ=255)
                                                                                  gTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGa                                                                                                         >  2:300068/1‑109 (MQ=255)
                                                                                  gTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGa                                                                                                         <  1:300068/109‑1 (MQ=255)
                                                                                      tcaatcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGt                                                                     <  1:1408764/141‑1 (MQ=255)
                                                                                      tcaatcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGt                                                                     >  2:1408764/1‑141 (MQ=255)
                                                                                        aatcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcaga                                                                                                               >  2:211429/1‑97 (MQ=255)
                                                                                        aatcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcaga                                                                                                               <  1:211429/97‑1 (MQ=255)
                                                                                        aatcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTAcc                                                                                  >  1:1145006/1‑126 (MQ=255)
                                                                                        aatcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTAcc                                                                                  <  2:1145006/126‑1 (MQ=255)
                                                                                        aatcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGgcgc                                                          <  2:244025/150‑1 (MQ=255)
                                                                                          tcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCt                                                                                     >  1:260778/1‑121 (MQ=255)
                                                                                          tcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCt                                                                                     <  2:260778/121‑1 (MQ=255)
                                                                                            aaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGa                                                      >  2:1158469/1‑150 (MQ=255)
                                                                                                      tgtgtgATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGcc                                            <  2:75252/150‑1 (MQ=255)
                                                                                                      tgtgtgATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGcc                                            <  2:158290/150‑1 (MQ=255)
                                                                                                          tgATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagca                                                                                                            <  1:586307/82‑1 (MQ=255)
                                                                                                          tgATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagca                                                                                                            >  2:586307/1‑82 (MQ=255)
                                                                                                          tgATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTgg                                        <  2:722947/150‑1 (MQ=255)
                                                                                                                     aCAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTg                             >  1:141590/1‑150 (MQ=255)
                                                                                                                       aGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagca                                                                                                            >  2:70249/1‑69 (MQ=255)
                                                                                                                       aGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagca                                                                                                            <  1:70249/69‑1 (MQ=255)
                                                                                                                         gTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgccgg                         >  2:1154119/1‑150 (MQ=255)
                                                                                                                         gTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgccgg                         <  1:211235/150‑1 (MQ=255)
                                                                                                                         gTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgccgg                         <  1:394531/150‑1 (MQ=255)
                                                                                                                          ttCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTg                                                                             <  2:86948/97‑1 (MQ=255)
                                                                                                                          ttCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTg                                                                             >  1:86948/1‑97 (MQ=255)
                                                                                                                             gATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgccgggccg                     >  1:1091832/1‑150 (MQ=255)
                                                                                                                             gATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgccg                          >  1:700635/1‑145 (MQ=255)
                                                                                                                             gATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgccg                          <  2:700635/145‑1 (MQ=255)
                                                                                                                                   ggTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGGAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGa               <  1:1154119/150‑1 (MQ=255)
                                                                                                                                     tGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGc             >  1:516273/1‑150 (MQ=255)
                                                                                                                                     tGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTAGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGc             <  1:36594/150‑1 (MQ=255)
                                                                                                                                        gcgcATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGa          <  2:1353160/150‑1 (MQ=255)
                                                                                                                                          gcATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATg        >  1:504022/1‑150 (MQ=255)
                                                                                                                                            cccggCCAAGGCGATCGCCATGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGTGGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATgcg      >  2:274781/4‑150 (MQ=255)
                                                                                                                                               ggCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGt   <  1:979645/150‑1 (MQ=255)
                                                                                                                                               ggCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGt   <  1:1362612/150‑1 (MQ=255)
                                                                                                                                                gCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTggcggc                                    >  2:748098/1‑116 (MQ=255)
                                                                                                                                                gCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTggcggc                                    <  1:748098/116‑1 (MQ=255)
                                                                                                                                                gCCAAGGCGATCGCCATGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGa                                                                                        <  1:1230108/64‑1 (MQ=38)
                                                                                                                                                gCCAAGGCGATCGCCATGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGa                                                                                        >  2:1230108/1‑64 (MQ=38)
                                                                                                                                                gCCAAGGCGATCGCCATGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGCACTTCCGGGCCGGGTCGAGCCGATGCGCGTc  >  1:535888/1‑150 (MQ=255)
                                                                                                                                                   aTGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTg                                         >  2:453508/1‑108 (MQ=255)
                                                                                                                                                   aTGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTg                                         <  1:453508/108‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCTGAAGGGGTCATTAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTC  >  NC_002947/3212872‑3213165

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