Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 830,866 C→G 7.1% A47G (GCG→GGG) ‡ PP_0715 → outer membrane efflux protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947830,8660CG 7.1% 177.7 / 3.2 56A47G (GCG→GGG) ‡PP_0715outer membrane efflux protein
Reads supporting (aligned to +/- strand):  ref base C (28/24);  new base G (2/2);  total (30/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 = 8.39e-01

CCTGAGCCATGCCCAACCGCTGCTATGGCCTGGGGCTGCTGATGGTGCTGGCCGGTTGCACTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGCGTATCGCTGAAGCCCGTGCCCAGTTGGCCATCGTGCA  >  NC_002947/830719‑831013
                                                                                                                                                   |                                                                                                                                                    
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                                                                                         >  1:104648/1‑148 (MQ=255)
                                                            cTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCt                                                                                         >  8:97260/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                                                                              >  4:16425/1‑146 (MQ=255)
                                                                         ggATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAgg                                                                              <  3:16425/146‑1 (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                                                          >  6:362062/1‑144 (MQ=255)
                                                                                               ccGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAAc                                                          <  5:362062/144‑1 (MQ=255)
                                                                                                     cTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGtggc                                               <  5:198020/149‑1 (MQ=255)
                                                                                                       ggACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGtggctg                                             >  8:32400/1‑149 (MQ=255)
                                                                                                       ggACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGtggctg                                             >  3:235389/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                                       >  1:63125/1‑126 (MQ=255)
                                                                                                                                    aggcaggGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGc                                       <  2:63125/126‑1 (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)
                                                                                                                                                  gcggcggcgCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGCGTATCGCTGAAGCCCGTGCCCAGTTGGCCATCGTGCa  <  7:32400/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCTGAGCCATGCCCAACCGCTGCTATGGCCTGGGGCTGCTGATGGTGCTGGCCGGTTGCACTCAGGTTGGTCCGGATTTCGAGAAGCCCCAAGCGCCGTGGCTGGACGGTTGGAGCACGCCCATGATCGAGCAGGCAGGGCGCAGCGCGGCGGCGCCCGACCTGCGTCA‑GTGGTGGGCGGTGTTTGCCGACTCGACCCTGGACGCGCTGATAGCCGAGGCCGATGCCAACAACACCAACCTGCGCGTGGCTGGCCTGCGTATCGCTGAAGCCCGTGCCCAGTTGGCCATCGTGCA  >  NC_002947/830719‑831013

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