Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 830,862 A→C 6.8% S46R (AGC→CGC)  PP_0715 → outer membrane efflux protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947830,8620AC 6.8% 180.7 / 4.8 59S46R (AGC→CGC) PP_0715outer membrane efflux protein
Reads supporting (aligned to +/- strand):  ref base A (31/24);  new base C (2/2);  total (33/26)
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.17e-01

CCTGAGCCATGCCCAACCGCTGCTATGGCCTGGGGCTGCTGATGGTGCTGGCCGGTTGCACTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGCGTATCGCTGAAGCCCGTGCCCAGTTGGCCATC  >  NC_002947/830719‑831008
                                                                                                                                               |                                                                                                                                                   
ccTGAGCCATGCCCAACCGCTGCTATGGCCTGGGGCTGCTGATGGTGCTGGCCGGTTGCACTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGcgcg                                                                                                                                                >  3:302383/1‑149 (MQ=255)
     gCCATGCCCAACCGCTGCTATGGCCTGGGGCTGCTGATGGTGCTGGCCGGTTGCACTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCgcggcggc                                                                                                                                           >  8:305452/1‑149 (MQ=255)
      ccATGCCCAACCGCTGCTATGGCCTGGGGCTGCTGATGGTGCTGGCCGGTTGCACTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCgcggcggcg                                                                                                                                          >  7:211862/1‑149 (MQ=255)
            ccAACCGCTGCTATGGCCTGGGGCTGCTGATGGTGCTGGCCGGTTGCACTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGCCGCAGCGCGGCGGCGCCCGAc                                                                                                                                    >  4:40936/1‑149 (MQ=255)
             cAACCGCTGCTATGGCCTGGGGCTGCTGATGGTGCTGGCCGGGTGCACTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGTAGCACGCCCATGATCGAGCAGGCAGGGCGCCGCCGCGCTGCGccctgcc                                                                                                                                   <  1:41489/149‑5 (MQ=255)
              aaCCGCTGCTATGGCCTGGGGCTGCTGATGGTGCTGGCCGGTTGCACTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCCGGGCGCAGCGCGGCGGCGCCCGCCCt                                                                                                                                  >  8:277791/1‑149 (MQ=255)
                ccGCTGCTATGGCCTGGGGCTGCTGATGGTGCTGGCCGGTTGCACTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGTCGGGCGCCGCCGCGCTGCGCCctgcctgc                                                                                                                                <  3:280423/149‑1 (MQ=38)
                   ctgctATGGCCTGGGGCTGCTGATGGTGCTGGCCGGTTGCACTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGt                                                                                                                              >  6:141659/1‑148 (MQ=255)
                        aTGGCCTGGGGCTGCTGATGGTGCTGGCCGGTTGCACTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑gtgg                                                                                                                       >  4:55510/1‑149 (MQ=255)
                             cTGGGGCTGCTGATGGTGCTGGCCGGTTGCACTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGc                                                                                                                  <  5:34683/149‑1 (MQ=255)
                                   ctgctgATGGTGCTGGCCGGTTGCACTCAGGTTGGTCCGGATTTCGAGAAGCCACAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGtt                                                                                                            >  5:72146/1‑149 (MQ=255)
                                   ctgctgATGGGGCTGGCCGGTTGCACTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGtt                                                                                                            <  5:148393/149‑1 (MQ=255)
                                      ctgATGGTGCTGGCCGGTTGCACTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGttt                                                                                                           <  1:101176/147‑1 (MQ=255)
                                      ctgATGGTGCTGGCCGGTTGCACTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGc                                                                                                         >  3:146170/1‑149 (MQ=255)
                                                           aCTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCt                                                                                    >  1:202912/1‑149 (MQ=255)
                                                            cTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTTGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTg                                                                                   <  8:43921/149‑1 (MQ=255)
                                                            cTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCt                                                                                    >  8:97260/1‑148 (MQ=255)
                                                            cTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCt                                                                                    >  1:104648/1‑148 (MQ=255)
                                                               aGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATa                                                                                <  7:131373/149‑1 (MQ=255)
                                                                  ttGGGCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGCCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGcc                                                                             <  3:55510/149‑1 (MQ=255)
                                                                       ccGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGc                                                                        >  4:351395/1‑149 (MQ=255)
                                                                         ggATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAgg                                                                         <  3:16425/146‑1 (MQ=255)
                                                                         ggATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAgg                                                                         >  4:16425/1‑146 (MQ=255)
                                                                          gATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCCGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTCATAGCCGAGGCCGa                                                                     <  7:305452/149‑1 (MQ=255)
                                                                           aTTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGAt                                                                    <  7:285173/149‑1 (MQ=255)
                                                                            tttCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATg                                                                   >  1:122345/1‑149 (MQ=255)
                                                                               cGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGcc                                                                 <  2:202912/148‑1 (MQ=255)
                                                                                gagaAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCcaa                                                               <  4:302383/149‑1 (MQ=255)
                                                                                 agaAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCcaac                                                              <  7:97260/149‑1 (MQ=255)
                                                                                  gaAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAg                                                                               >  4:233025/1‑131 (MQ=255)
                                                                                  gaAGCCCCAAGCGCAGTGGCTGGACGGTTGGAGCACGCCCAAGATCGAGCAGGCAGTGCGCAGCGCGGCGGCGCCCCACCTGCGTCA‑GTGGTGGGAGGTGTTTGCCGACTAGACCCTGGACGCGCTGATAGCCGAGGCCGATGCcaac                                                              >  7:39157/1‑148 (MQ=255)
                                                                                   aaGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCcaacaa                                                            >  1:761/1‑149 (MQ=255)
                                                                                   aaGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCGGATAGCCGAGGCCGATGCcaacaa                                                            <  8:211862/149‑1 (MQ=255)
                                                                                     gCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAacaaca                                                          >  8:188295/1‑149 (MQ=255)
                                                                                      ccccAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAacac                                                         >  2:125536/1‑149 (MQ=255)
                                                                                        ccAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCa                                                       >  6:198020/1‑149 (MQ=255)
                                                                                               ccGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAAc                                                     <  5:362062/144‑1 (MQ=255)
                                                                                               ccGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAAc                                                     >  6:362062/1‑144 (MQ=255)
                                                                                                     cTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGtggc                                          <  5:198020/149‑1 (MQ=255)
                                                                                                       ggACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGtggctg                                        >  3:235389/1‑149 (MQ=255)
                                                                                                       ggACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGtggctg                                        >  8:32400/1‑149 (MQ=255)
                                                                                                         aCGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGtggctggc                                      <  3:351395/149‑1 (MQ=255)
                                                                                                            gTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTCGTGGGAGGTGTTTGCCGACTCGACACTGGACCCGCTGATAGCCGAGGCCGATGCCAACAACACCAACATGAGCGTGGCTGGACTg                                   >  7:308404/1‑149 (MQ=255)
                                                                                                                        ccATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGCGTATCGCTg                         >  7:237497/1‑147 (MQ=255)
                                                                                                                              tgacgcaggtcgGGCGCCGCCGCGCTGCGCCCTGCCTGCGTCAGGTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCACCCTGCGCGTGGCTGGCCTGCGTATCGCTGAAGCCCg                  >  1:255542/5‑149 (MQ=37)
                                                                                                                               cGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGCGTATCGCTGaa                       >  6:314480/1‑142 (MQ=255)
                                                                                                                               cGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGCGTATCGCTGaa                       <  5:314480/142‑1 (MQ=255)
                                                                                                                                acgcaggtcgGGCGCCGCCGCGCTGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGCGTATCGCTGAAGCCCGTGc               >  2:287579/3‑149 (MQ=255)
                                                                                                                                   ccggctggGCGCAGCGCGGCGGCGCCCGACGTGCGTGA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTCGACGCGCTTATAGCCGACGCCGATGCCAACCACAGCAACCTCCGCGTGGCTGGCCTGCGTATCGCTGAAGCCCGTGCCCa            <  8:39157/147‑1 (MQ=255)
                                                                                                                                   caggcaggGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGCGTATCGCTGAAGCCCGTGCCCa            <  6:72146/149‑1 (MQ=255)
                                                                                                                                    aggcaggGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGc                                  <  2:63125/126‑1 (MQ=255)
                                                                                                                                    aggcaggGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGc                                  >  1:63125/1‑126 (MQ=255)
                                                                                                                                    aggcaggGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGCGTATCGCTGAAGCCCGTGCCCAg           >  4:9945/1‑149 (MQ=255)
                                                                                                                                      gcaggGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGCGTATCGCTGAAGCCCGTGCCCAGtt         <  4:235389/149‑1 (MQ=255)
                                                                                                                                            cgcAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGCCTCGCCCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGCGTATCGCTGAAGCCCGTGCCCAGTTGGCCAt   >  6:290819/1‑149 (MQ=255)
                                                                                                                                            cgcAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGCGTATCGCTGAAGCCCGTGCCCAGTTGGCCAt   <  2:104648/149‑1 (MQ=255)
                                                                                                                                             gcAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGCGTATCGCTGAAGCCCGTGCCCAGtt         <  5:70486/142‑1 (MQ=255)
                                                                                                                                             gcAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGCGTATCGCTGAAGCCCGTGCCCAGTTGGCCATc  <  2:122345/149‑1 (MQ=255)
                                                                                                                                             gcAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACACTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGCGTATCGCTGAAGCCCGTGCCCAGtt         >  6:70486/1‑142 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CCTGAGCCATGCCCAACCGCTGCTATGGCCTGGGGCTGCTGATGGTGCTGGCCGGTTGCACTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGCGTATCGCTGAAGCCCGTGCCCAGTTGGCCATC  >  NC_002947/830719‑831008

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑