Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,664,248 G→A 8.1% G134D (GGC→GAC)  nuoH → NADH‑quinone oxidoreductase subunit H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,664,2480GA 8.1% 212.2 / 5.7 74G134D (GGC→GAC) nuoHNADH‑quinone oxidoreductase subunit H
Reads supporting (aligned to +/- strand):  ref base G (30/38);  new base A (2/4);  total (32/42)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

CTGGCACCGGTAGTGGCCATGAGCGCCCTGCTGATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGCAGGTGGGTTCGTTCAACATGCGCGACATCGTTGAGTAC  >  NC_002947/4664100‑4664393
                                                                                                                                                      |                                                                                                                                                 
cTGGCACCGGTAGTGGCCATGAGCGCCCTGCTGATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGACATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCgg                                                                                                                                                   >  6:236663/1‑149 (MQ=255)
  ggCACCGGTAGTGGCCATGAGCGCCCTGCTGATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCt                                                                                                                                                 >  7:203684/1‑149 (MQ=255)
   gCACCGGTAGTGGCCATGAGCGCCCTGCTGATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTg                                                                                                                                                <  2:239352/149‑1 (MQ=255)
   gCACCGGGAGTGGCCATGAGCGCCCTGCTGATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTg                                                                                                                                                <  8:162299/149‑1 (MQ=255)
    cACCGGTAGTGGCCATGAGCGCCCTGCTGATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTgg                                                                                                                                               <  5:53465/149‑1 (MQ=255)
        ggTAGTGGCCATGAGCGCCCTGCTGATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGgtcgt                                                                                                                                           <  5:84111/149‑1 (MQ=255)
        ggTAGTGGCCATGAGCGCCCTGCTGATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGgtcgt                                                                                                                                           <  3:379564/149‑1 (MQ=255)
        ggTAGTGGCCATGAGCGCCCTGCTGATCGCCTGCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGgtcgt                                                                                                                                           <  2:305639/149‑1 (MQ=255)
         gTAGTGGCCATGAGCGCCCTGCTGATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGgtcgtc                                                                                                                                          <  5:329384/149‑1 (MQ=255)
               gCCATGAGCGCCCTGCTGATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGaacaa                                                                                                                                    <  3:270187/149‑1 (MQ=255)
               gCCATGAGCGCCCTGCTGATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGaac                                                                                                                                      >  8:19896/1‑147 (MQ=255)
               gCCATGAGCGCCCTGCTGATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGaac                                                                                                                                      <  7:19896/147‑1 (MQ=255)
                ccATGAGCGCCCTGCTGATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGaacaac                                                                                                                                   <  7:103799/149‑1 (MQ=255)
                ccATGAGCGCCCTGCTGATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGACGACCAGCCGGCGaacagc                                                                                                                                   <  1:71781/149‑3 (MQ=255)
                 cATGAGCGCCCTGCTGATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTTGGTCTACGCGGTGCTGTTCGACGACCAGCCGGCGAacagca                                                                                                                                  <  3:82218/149‑1 (MQ=255)
                  aTGAGCGCCCTGCTGATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAAcaacaa                                                                                                                                 <  7:31307/149‑1 (MQ=255)
                      gcgcCCTGCTGATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTAc                                                                                                                             >  3:221951/1‑149 (MQ=255)
                                gATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTgg                                                                                                                   >  6:381025/1‑149 (MQ=255)
                                gATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGACCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTgg                                                                                                                   >  7:209200/1‑149 (MQ=255)
                                 aTCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTCTACGCCCGCTGCTCGTCGAACAACAAGTACGCCCTGCTggg                                                                                                                  >  2:72122/1‑149 (MQ=255)
                                 aTCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACATGAACATCGGCCTGCTGTTCTTCTTAGACATGGACGGCCTGTCGTTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTggg                                                                                                                  >  7:307523/1‑149 (MQ=255)
                                    gCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCa                                                                                                                <  2:235628/148‑1 (MQ=255)
                                     ccTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGc                                                                                                              >  5:11894/1‑149 (MQ=255)
                                            ggTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTg                                                                                                       >  6:74006/1‑149 (MQ=255)
                                                aTCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTc                                                                                                   <  1:199400/149‑1 (MQ=255)
                                                 tCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGACGACCAGCCGGCGaaca                                                                                                                                     <  7:279237/114‑1 (MQ=255)
                                                 tCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGACGACCAGCCGGCGaaca                                                                                                                                     >  8:279237/1‑114 (MQ=255)
                                                         aCCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAGGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGa                                                                                           <  5:225108/148‑1 (MQ=255)
                                                          ccccGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGAcc                                                                                         >  5:54163/1‑149 (MQ=255)
                                                          ccccGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTGCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTCCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTGCCGTGCTTCGGCACAgc                                                                                           >  6:212784/1‑146 (MQ=255)
                                                           cccGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCg                                                                                        >  3:92812/1‑149 (MQ=255)
                                                           cccGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCg                                                                                        >  5:166392/1‑149 (MQ=255)
                                                                    gggCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCCTGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGGCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTAc                                                                                <  1:187181/148‑1 (MQ=255)
                                                                      gCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGaa                                                                             <  6:54163/149‑1 (MQ=255)
                                                                             cGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTcc                                                                      >  7:194529/1‑149 (MQ=255)
                                                                                ccTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTgg                                                                   <  1:3181/149‑1 (MQ=255)
                                                                                 cTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTgg                                                                   <  8:194529/148‑1 (MQ=255)
                                                                                  tGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGATAGTTCCCCGGCTGGTCGTCGAACAACAATTTCGCCCTGCTGGGCATGTTGCGTGCTTCGGCACAGACCGTGTCGTAACAACTTTTCCTggg                                                                  >  1:107833/1‑149 (MQ=255)
                                                                                       aTCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGgc                                                            >  1:142307/1‑149 (MQ=255)
                                                                                        tCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGgcg                                                           <  8:203684/149‑1 (MQ=255)
                                                                                          ggCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGgcgc                                                          <  7:303226/148‑1 (MQ=255)
                                                                                          ggCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCt                                                         <  8:90682/149‑1 (MQ=255)
                                                                                           gCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACTCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCt                                                         >  5:261755/1‑148 (MQ=255)
                                                                                           gCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCt                                                         >  3:293306/1‑148 (MQ=255)
                                                                                            ccTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGa                                                       <  6:166392/149‑1 (MQ=255)
                                                                                            ccTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGa                                                       >  8:56156/1‑149 (MQ=255)
                                                                                             ctgctgTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGAt                                                      >  3:291555/1‑149 (MQ=255)
                                                                                                ctgTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTATCGCGGTGCTGTTCGCCGACTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATgg                                                    >  1:257260/1‑149 (MQ=255)
                                                                                                ctgTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATgg                                                    <  5:236663/148‑1 (MQ=255)
                                                                                                 tgTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGc                                                  >  7:187035/1‑149 (MQ=255)
                                                                                                      ttcttcGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCgtggt                                             >  7:72358/1‑149 (MQ=255)
                                                                                                        cttcGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCgtggtg                                            <  6:261755/148‑1 (MQ=255)
                                                                                                         ttcGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCATGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCgtggtggt                                          >  6:101808/1‑149 (MQ=255)
                                                                                                           cGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGACTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCgtggtggt                                          <  2:257260/147‑1 (MQ=255)
                                                                                                            gCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGcc                                                                <  6:154551/124‑1 (MQ=255)
                                                                                                            gCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGcc                                                                >  5:154551/1‑124 (MQ=255)
                                                                                                            gCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGc                                        <  6:11894/148‑1 (MQ=255)
                                                                                                               aTGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGCAGGt                                    <  5:74006/149‑1 (MQ=255)
                                                                                                               aTGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGCAGGt                                    <  6:355307/149‑1 (MQ=255)
                                                                                                                 ggccggccTGTCGGGCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGCAGGTgg                                  <  4:221951/149‑1 (MQ=255)
                                                                                                                   ccggccTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGCAGGTGGgt                                <  8:209200/149‑1 (MQ=255)
                                                                                                                    cggccTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGCAGGTGGgtt                               >  3:367236/1‑149 (MQ=255)
                                                                                                                        cTGTCGGTCTACGCGGTGCTGTACGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGATCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGCAGGTGGgttcgt                            >  1:98484/1‑148 (MQ=255)
                                                                                                                         tGTAGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACTAAGTGTTCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGCAg                                      <  4:139106/137‑1 (MQ=255)
                                                                                                                          gTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGCAGGTGGgttcgttc                          >  2:163951/1‑148 (MQ=255)
                                                                                                                           tCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTTCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGCAGGTGGGTTCGTTCa                         <  5:101808/148‑1 (MQ=255)
                                                                                                                           tCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGCAGGTGGgttcgttc                          <  2:107833/147‑1 (MQ=255)
                                                                                                                                    cgcgGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGCAGGTGGGTTCGTTCAACATgcgc                 <  2:142307/148‑1 (MQ=255)
                                                                                                                                    cgcgGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGCAGGTGGGTTCGTTCAACATgcgc                 <  1:72122/148‑1 (MQ=255)
                                                                                                                                               ttCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGCAGGTGGGTTCGTTCAACATGCGCGACATCGtt        <  6:42871/147‑1 (MQ=255)
                                                                                                                                               ttCGCCGGCTGGTCGGCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGCAGGTGGGTTCGTTCAACATGCGCGACATCGtt        <  4:293306/147‑1 (MQ=255)
                                                                                                                                                  gCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGCAGGTGGGTTCGTTCAACATGCGCGACATCGTTGAGTa   >  5:69472/1‑149 (MQ=255)
                                                                                                                                                   ccGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGCAGGTGGGTTCGTTCAACATGCGCGACATCGTTGAGTAc  >  4:133655/1‑149 (MQ=255)
                                                                                                                                                    cGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGc                                                  >  8:147077/1‑100 (MQ=255)
                                                                                                                                                    cGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGc                                                  <  7:147077/100‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
CTGGCACCGGTAGTGGCCATGAGCGCCCTGCTGATCGCCTTCGTGGTCATCCCGATCACCCCGACCTGGGGCGTTGCCGACCTGAACATCGGCCTGCTGTTCTTCTTCGCCATGGCCGGCCTGTCGGTCTACGCGGTGCTGTTCGCCGGCTGGTCGTCGAACAACAAGTACGCCCTGCTGGGCAGCTTGCGTGCTTCGGCACAGACCGTGTCGTACGAAGTGTTCCTGGGCCTGGCGCTGATGGGCGTGGTGGTGCAGGTGGGTTCGTTCAACATGCGCGACATCGTTGAGTAC  >  NC_002947/4664100‑4664393

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: