Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,069,288 A→C 7.2% L107R (CTG→CGG)  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,2880AC 7.2% 214.7 / 7.7 69L107R (CTG→CGG) PP_5323M23/M37 family peptidase
Reads supporting (aligned to +/- strand):  ref base A (35/29);  new base C (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 = 9.80e-01

ACGGGAACGCATAGCGGAAGCCCTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTCGCGTTCAGCAGTACCTGG  >  NC_002947/6069151‑6069434
                                                                                                                                          |                                                                                                                                                    
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                                                                                                                     <  5:103861/149‑1 (MQ=255)
                      cTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCtt                                                                                                                     >  3:122639/1‑149 (MQ=255)
                           gCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCg                                                                                                                >  1:30787/1‑149 (MQ=255)
                           gCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCg                                                                                                                >  7:19526/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                                                                                                             <  2:142591/133‑1 (MQ=255)
                                              cATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTg                                                                                                             >  1:142591/1‑133 (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:30787/149‑1 (MQ=255)
                                                                                            cgccCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGaaaa                                             <  2:207790/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                          >  8:334783/6‑143 (MQ=255)
                                                                                                                     gtatcgtgCGCCGGCTGCTGCCGCC‑GCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAg                          <  7:334783/138‑1 (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                    <  3:91658/148‑1 (MQ=255)
                                                                                                                      cGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTc                    <  4:95442/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)
                                                                                                                                        gcagcCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTCGCGTTCAGCAGTACCTgg  <  6:300022/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
ACGGGAACGCATAGCGGAAGCCCTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTCGCGTTCAGCAGTACCTGG  >  NC_002947/6069151‑6069434

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