Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,069,285 A→T 7.2% L108Q (CTG→CAG)  PP_5323 ← M23/M37 family peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,069,2850AT 7.2% 210.1 / 5.3 69L108Q (CTG→CAG) PP_5323M23/M37 family peptidase
Reads supporting (aligned to +/- strand):  ref base A (35/29);  new base T (3/2);  total (38/31)
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.42e-01

GCCCGCCAATCCACGGGAACGCATAGCGGAAGCCCTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTCGCGTTC  >  NC_002947/6069139‑6069422
                                                                                                                                                   |                                                                                                                                           
gCCCGCCAATCCACGGGAACGCATAGCGGAAGCCCTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGgcagc                                                                                                                                           <  8:377611/149‑1 (MQ=255)
            aCGGGAACGCATAGCGGAAGCCCTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCGCCGGCTGCTGCCGCC‑GCGCAcc                                                                                                                               <  3:94235/149‑2 (MQ=255)
                gAACGCATAGCGGAAGCCCTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGgcagca                                                                                                                                          >  4:33374/1‑134 (MQ=255)
                gAACGCATAGCGGAAGCCCTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGgcagca                                                                                                                                          <  3:33374/134‑1 (MQ=255)
                gAACGCATAGCGGAAGCCCTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATAc                                                                                                                           >  6:103861/1‑149 (MQ=255)
                    gCATAGCGGAAGCCCTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCgg                                                                                                                       >  2:389399/1‑149 (MQ=255)
                     cATAGCGGAAGCCCTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGGCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGGGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCggg                                                                                                                      <  1:303894/149‑1 (MQ=255)
                                  cTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCtt                                                                                                         >  3:122639/1‑149 (MQ=255)
                                  cTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCtt                                                                                                         <  5:103861/149‑1 (MQ=255)
                                       gCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCg                                                                                                    >  7:19526/1‑149 (MQ=255)
                                       gCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCg                                                                                                    >  1:30787/1‑149 (MQ=255)
                                           tttGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAACAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTtgt                                                                                                       >  3:394874/1‑142 (MQ=255)
                                           tttGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAACAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTtgt                                                                                                       <  4:394874/142‑1 (MQ=255)
                                              gCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGAc                                                                                               >  8:61516/1‑147 (MQ=255)
                                               cccGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCACCGCCTTGTGCGTTGACCAg                                                                                            <  4:298332/149‑1 (MQ=255)
                                                ccGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAgg                                                                                           <  5:252742/149‑1 (MQ=255)
                                                   ggTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCGCCTGAGCACAATGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGACGCCTTGTGGGTTGACAAGGTTg                                                                                        >  4:319614/1‑149 (MQ=255)
                                                   cgTCACTCAGGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTg                                                                                        <  5:101969/148‑1 (MQ=255)
                                                       aCCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCa                                                                                    <  2:353379/149‑1 (MQ=255)
                                                        cccATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGc                                                                                                                    <  4:292888/116‑1 (MQ=255)
                                                        cccATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGc                                                                                                                    >  3:292888/1‑116 (MQ=255)
                                                        cccATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAg                                                                                   >  4:50786/1‑149 (MQ=255)
                                                         ccATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGc                                                                                  >  7:249499/1‑149 (MQ=255)
                                                         ccATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGc                                                                                  >  8:115229/1‑149 (MQ=255)
                                                          cATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTg                                                                                                 >  1:142591/1‑133 (MQ=255)
                                                          cATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTg                                                                                                 <  2:142591/133‑1 (MQ=255)
                                                          cATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGcc                                                                                 <  3:50786/149‑1 (MQ=255)
                                                           agggCCTGCAGCAGGGTCGAACGGTGGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGTATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGc                                                                                  <  4:345115/145‑1 (MQ=255)
                                                           aTGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCg                                                                                >  3:95442/1‑149 (MQ=255)
                                                              gCCTGCAGGAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCaa                                                                             <  1:352739/149‑1 (MQ=255)
                                                               ccTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCaa                                                                             >  4:278573/1‑148 (MQ=255)
                                                                   cagcagGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAActcta                                                                        <  4:341575/149‑1 (MQ=255)
                                                                   cagcagGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAActct                                                                         <  3:2934/148‑1 (MQ=255)
                                                                      cagGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCTTTTTGTCCAGCCGCAActctacct                                                                     <  7:61516/149‑1 (MQ=255)
                                                                               aCGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGca                                                          >  4:91658/1‑149 (MQ=255)
                                                                                cGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGca                                                          >  8:46674/1‑148 (MQ=255)
                                                                                 ggTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTTCGTGGTGCG‑CGGCGGAAGCAGCCGGCGCAAGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCAAGc                                                                                  >  7:223526/1‑125 (MQ=255)
                                                                                 ggTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGACATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGcaca                                                        <  1:389399/149‑1 (MQ=255)
                                                                                     gTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCa                                                    <  4:122639/149‑1 (MQ=255)
                                                                                                gcCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGt                                         >  1:207790/1‑149 (MQ=255)
                                                                                                gcCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGt                                         <  3:278573/149‑1 (MQ=255)
                                                                                                 cccGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTCGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTg                                        >  2:122257/1‑149 (MQ=255)
                                                                                                        cgccCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGaaaa                                 <  2:207790/149‑1 (MQ=255)
                                                                                                        cgccCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGaaaa                                 <  2:30787/149‑1 (MQ=255)
                                                                                                          cccccTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGc                               >  5:2403/1‑149 (MQ=255)
                                                                                                             ccTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCa                            >  7:351396/1‑149 (MQ=255)
                                                                                                                gAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGaaaa                                 <  3:394815/141‑1 (MQ=255)
                                                                                                                gAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGaaaa                                 >  4:394815/1‑141 (MQ=255)
                                                                                                                gAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGc                         >  4:43514/1‑149 (MQ=255)
                                                                                                                gAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTACACCAGGTGATCCCCCCGCAACTCTACCTC‑T‑ACCCGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGc                         >  1:132255/1‑149 (MQ=255)
                                                                                                                  gCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCa                       >  5:300022/1‑149 (MQ=255)
                                                                                                                   cacaCTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCg                                     >  4:266953/1‑134 (MQ=255)
                                                                                                                   cacaCTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCg                                     <  3:266953/134‑1 (MQ=255)
                                                                                                                    acacTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCAcc                     >  2:84313/1‑149 (MQ=255)
                                                                                                                     cacTGAGCATCTGCGTGGTTGG‑GGGGGGCAGCAGGCGGCGGACGTTACGCGGGGCATGGCCGGCTTGTTGGTTGTCCAGGTTTTCCATCCGCAACTCTAACTC‑TCACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCAcc                     <  1:84313/149‑1 (MQ=255)
                                                                                                                     cacTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGAAGCTCGTTTCGTGCCGAAAAGCTCACGCCAcag                    >  6:270878/1‑147 (MQ=255)
                                                                                                                      acTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGt                   >  4:336235/1‑149 (MQ=255)
                                                                                                                          aGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCg                 >  5:111607/1‑147 (MQ=255)
                                                                                                                          aGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGaa               <  7:115229/149‑1 (MQ=255)
                                                                                                                           gCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACTCCACCGTCGAAg              <  6:2403/149‑1 (MQ=255)
                                                                                                                                tGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCAc                      <  1:364946/136‑1 (MQ=255)
                                                                                                                                tGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCAc                      >  2:364946/1‑136 (MQ=255)
                                                                                                                                 gtatcgtgCGCCGGCTGCTGCCGCC‑GCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAg              <  7:334783/138‑1 (MQ=255)
                                                                                                                                 gtatcgtgCGCCGGCTGCTGCCGCC‑GCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAg              >  8:334783/6‑143 (MQ=255)
                                                                                                                                 gCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTc        >  2:382976/1‑149 (MQ=255)
                                                                                                                                 gCGTGGTGCG‑CGGCGGCAACAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTc        >  1:392607/1‑149 (MQ=255)
                                                                                                                                  cGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTc        <  4:95442/148‑1 (MQ=255)
                                                                                                                                  cGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTc        <  3:91658/148‑1 (MQ=255)
                                                                                                                                   gtggtgCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACACCACCGTCGAAGGTTTTcgc      >  1:215042/1‑149 (MQ=255)
                                                                                                                                      gtgCGCCGGCTGCTGCCGCC‑GCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTCGCGtt   >  6:28729/1‑149 (MQ=255)
                                                                                                                                       tgcgccggctgCTGCCGCC‑GCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTCGCGTTc  >  1:394801/6‑149 (MQ=255)
                                                                                                                                       tgCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGGCGTCCAGCCGCAACTCTACCTCAC‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTCGCGtt   <  1:122257/149‑1 (MQ=255)
                                                                                                                                                   agcagcCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGca                                                          >  3:159160/1‑82 (MQ=255)
                                                                                                                                                   agcagcCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGca                                                          <  4:159160/82‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
GCCCGCCAATCCACGGGAACGCATAGCGGAAGCCCTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTCGCGTTC  >  NC_002947/6069139‑6069422

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