Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 228,897 A→G 6.2% S575G (AGT→GGT)  PP_0176 → putative efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947228,8970AG 6.2% 189.2 / 2.2 64S575G (AGT→GGT) PP_0176putative efflux transporter
Reads supporting (aligned to +/- strand):  ref base A (32/28);  new base G (2/2);  total (34/30)
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

GGCGGCCAAGCGCATGACCCGCTTCGCCTGGCGCGACATCGTCGAAATGACCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACTGTGGTCGAGGAAGTCGGTGCGCACTACGCCGCCTGCCTGCGTGCTGGC  >  NC_002947/228749‑229046
                                                                                                                                                    |                                                                                                                                                     
ggcggcCAAGCGCATGACCCGCTTCGCCTGGCGCGACATCGTCGAAATGACCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAg                                                                                                                                                      <  2:758231/150‑1 (MQ=255)
       aaGCGCATGACCCGCTTCGCCTGGCGCGACATCGTCGAAATGACCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGtgtg                                                                                                                                               <  1:1240651/150‑1 (MQ=255)
       aaGCGCATGACCCCCTTCGCCTGGCGCGACATCGGTCAAATGACCGATCCTGCTACATTTGCCGAACACCGCCAGCTTGTCTTACAGATGCTCGACCGCCTGATCCAGCACCTGCCACGTCTTTCGCAAACCGTCCAAGACAgacgtatc                                                                                                                                               <  2:3016/150‑8 (MQ=255)
                 cccGCTTCGCCTGGCGCGACATCGTCGAAATGACCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGACGGATCTCCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTg                                                                                                                                     <  1:1055827/150‑1 (MQ=255)
                    gCTTCGCCTGGCGCGACATCGTCGAAATGACCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGAcc                                                                                                                                  >  1:867547/1‑150 (MQ=255)
                    gCTTCGCCTGGCGCGACATCGTCGAAATGACCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGAcc                                                                                                                                  >  1:867528/1‑150 (MQ=255)
                           cTGGCGCGACATCGTCGAAATGACCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGtgtg                                                                                                                                               <  2:1273528/130‑1 (MQ=255)
                           cTGGCGCGACATCGTCGAAATGACCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGtgtg                                                                                                                                               >  1:1273528/1‑130 (MQ=255)
                              gcgcGACATCGTCGAAATGACCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCTGCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAgg                                                                                                                        <  1:925957/150‑1 (MQ=255)
                                gcgACATCGTGGAAATGACCGAGCCTGCGACATTGGCCGAACAGCGCCAGGTGGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGt                                                                                                                      <  1:1019252/150‑1 (MQ=255)
                                        gtcgAAATGACCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTg                                                                                                              <  1:883416/150‑1 (MQ=255)
                                         tcgAAATGACCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGc                                                                                                             >  2:1033065/1‑150 (MQ=255)
                                         tcgAAATGACCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGc                                                                                                             <  2:545144/150‑1 (MQ=255)
                                          cgAAATGACCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCt                                                                                                            <  2:325461/150‑1 (MQ=255)
                                           gAAATGACCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTc                                                                                                           <  2:867547/150‑1 (MQ=255)
                                           gAAATGACCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTc                                                                                                           <  2:867528/150‑1 (MQ=255)
                                             aaTGACCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGTACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGa                                                                                                         <  1:805863/150‑1 (MQ=255)
                                                 aCCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGAc                                                                                                        <  2:910145/147‑1 (MQ=255)
                                                 aCCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGAc                                                                                                        >  1:910145/1‑147 (MQ=255)
                                                         tGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTa                                                                                             >  2:723183/1‑150 (MQ=255)
                                                           cGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACa                                                                                           >  2:1084849/1‑150 (MQ=255)
                                                           cGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCCCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACa                                                                                           >  2:209681/1‑150 (MQ=255)
                                                             aCATTGGCCGAGCATGGCCAGGGTGTCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATg                                                                                         >  2:16939/1‑150 (MQ=255)
                                                                ttGGCCGAACACCGCCAGGTTGGCGTACAGATGCCCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCcg                                                                                      <  1:209681/150‑1 (MQ=255)
                                                                 tGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGAGCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCcgc                                                                                     <  1:391742/150‑1 (MQ=255)
                                                                    ccGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGCCAGTAGCGTGGCGCCAGCTGACCGGCGCTGAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCcgcgcg                                                                                  >  1:472833/1‑150 (MQ=255)
                                                                       aaCACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACGGTGgtgt                                                                                                                                                <  2:105671/85‑1 (MQ=255)
                                                                       aaCACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACGGTGgtgt                                                                                                                                                >  1:105671/1‑85 (MQ=255)
                                                                          acCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGtgtg                                                                                                                                               >  1:202789/1‑83 (MQ=255)
                                                                          acCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGtgtg                                                                                                                                               <  2:202789/83‑1 (MQ=255)
                                                                           ccGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCca                                                                           <  1:723183/150‑1 (MQ=255)
                                                                             gCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCcagc                                                                         >  2:1249542/1‑150 (MQ=255)
                                                                               cAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCcagcag                                                                       >  1:465723/1‑150 (MQ=255)
                                                                                aGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGGGGTGTTGCACTGCGTGACCTGCGCGTAGGGTTGAACCTTCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAgg                                                                      <  1:922633/150‑1 (MQ=255)
                                                                                  gTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGa                                                                                                         >  1:721629/1‑113 (MQ=255)
                                                                                  gTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGa                                                                                                         <  2:721629/113‑1 (MQ=255)
                                                                                             gATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGac                                                                                                                  >  1:927192/1‑92 (MQ=255)
                                                                                             gATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGaa                                                                                                                  <  2:927192/93‑1 (MQ=255)
                                                                                             gATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCt                                                         >  2:952411/1‑150 (MQ=255)
                                                                                                       cGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACtgtg                                               >  1:702624/1‑150 (MQ=255)
                                                                                                       cGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACtgtg                                               <  2:465723/150‑1 (MQ=255)
                                                                                                        gCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGGGACCTGCGCCTCAGGTTGAACCTGATCGACCTGCTTGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCACCTGAACACTGTgg                                              >  2:175132/1‑150 (MQ=255)
                                                                                                        gCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGGGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACTGAgg                                              >  1:352742/1‑150 (MQ=255)
                                                                                                         ccTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAg                                                                                                                         <  1:612376/74‑1 (MQ=255)
                                                                                                         ccTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAg                                                                                                                         >  2:612376/1‑74 (MQ=255)
                                                                                                            gATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACTGTGGTCGa                                          >  2:1004283/1‑150 (MQ=255)
                                                                                                            gATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACTGTGGTCGa                                          <  1:833581/150‑1 (MQ=255)
                                                                                                             aTGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACTGTGGTCGAg                                         <  1:1004283/150‑1 (MQ=255)
                                                                                                               gcagcaCCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTTCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACTGTGGTCGAGGa                                       >  2:431763/1‑150 (MQ=255)
                                                                                                                     ccTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACTGTGGTCGAGGAAGTCgg                                 >  1:199615/1‑150 (MQ=255)
                                                                                                                               cTGTCGCAAACCGGCCAAGACGGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACTGTGGTCGAGGAAGTCGGTGCGCACTAc                       >  1:109869/1‑150 (MQ=255)
                                                                                                                                tGTCGCAAACCGGCCAAGACAGTGTTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGa                                                               >  1:653955/1‑109 (MQ=255)
                                                                                                                                tGTCGCAAACCGGCCAAGACAGTGTTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGa                                                               <  2:653955/109‑1 (MQ=255)
                                                                                                                                tGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACTGTGGTCGAGGAAGTCGGTGCGCACTAcg                      >  2:148714/1‑150 (MQ=255)
                                                                                                                                 gTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACTGTGGTCGAGGAAGTCGGTGCGCACTAcgc                     >  2:658676/1‑150 (MQ=255)
                                                                                                                                          cGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCcagca                                                                        <  2:449333/90‑1 (MQ=255)
                                                                                                                                          cGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCcagca                                                                        >  1:449333/1‑90 (MQ=255)
                                                                                                                                          cGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACTGTCGTCGAGGAAGTCGGTGCGCACTACGCCgcctgcct            >  2:685960/1‑150 (MQ=255)
                                                                                                                                            gCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCTGCCAGCAGGCTCGCGAGCGCCTGAACACTGTGGTCGAGGGAGTCGGTGCGCACTACGCCGCctgcctgc          >  1:84861/1‑150 (MQ=255)
                                                                                                                                              cAAGACGGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACTGTGGTCGAGGAAGTCGGTGCGCACTACGCCGCCTGCCTGCGt        <  2:109869/150‑1 (MQ=255)
                                                                                                                                              cAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACTGTGGTCGAGGAAGTCGGTGCGCACTACGCCGCCTGCCTGCGt        >  1:946600/1‑150 (MQ=255)
                                                                                                                                                  aCAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACTGTGGTCGAGGAAGTCGGt                                >  2:1168684/1‑122 (MQ=255)
                                                                                                                                                  aCAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACTGTGGTCGAGGAAGTCGGt                                <  1:1168684/122‑1 (MQ=255)
                                                                                                                                                  aCAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACTGTGGTCGAGGAAGTCGGTGCGCACTACGCCGCCTGCCTGCGTGCTg    >  2:389486/1‑150 (MQ=255)
                                                                                                                                                    aGTGGTGTGGCACTGCGTGACATGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACTGTGGTCGAGGAAGTCGGTGCGCACTACGCCGCCTGCCTGCGTGCTGGc  >  1:343324/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCGGCCAAGCGCATGACCCGCTTCGCCTGGCGCGACATCGTCGAAATGACCGAGCCTGCGACATTGGCCGAACACCGCCAGGTTGGCGTACAGATGCTCGACCGCCTGATGCAGCACCTGCCACGTCTGTCGCAAACCGGCCAAGACAGTGGTGTGGCACTGCGTGACCTGCGCGTAGGGTTGAACCTGCTCGACCTGCTGGCGTACATGCCGCGCGCCGGCCAGCAGGCTCGCGAGCGCCTGAACACTGTGGTCGAGGAAGTCGGTGCGCACTACGCCGCCTGCCTGCGTGCTGGC  >  NC_002947/228749‑229046

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: