Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,667,365 C→G 6.2% A377A (GCC→GCG nuoL → NADH‑quinone oxidoreductase subunit L

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,667,3650CG 6.2% 257.5 / 6.5 81A377A (GCC→GCGnuoLNADH‑quinone oxidoreductase subunit L
Reads supporting (aligned to +/- strand):  ref base C (35/40);  new base G (3/2);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

GGGACGCGGCGATCTTCCACCTGATGACCCACGCCTTCTTCAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTGGTAGGCGCGTTCATGA  >  NC_002947/4667223‑4667513
                                                                                                                                               |                                                                                                                                                    
gggACGCGGCGATCTTCCACCTGATGACCCACGCCTTCTTCAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGc                                                                                                                                                 >  8:81938/1‑148 (MQ=255)
gggACGCGGCGATCTTCCACCTGATGACCCACGCCTTCTTCAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCACAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGc                                                                                                                                                 >  8:219366/1‑148 (MQ=255)
      cggcgATCTTCCACCTGATGACCCACGCCTTCTTCAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTc                                                                                                                                          <  2:103858/149‑1 (MQ=255)
        gcgATCTTCCACCTGATGACCCACGCCTTCTTCAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCGTAGGCCAGcggcag                                                                                                                                        <  7:327817/149‑6 (MQ=255)
        gcgATCTTCCACCTGATGACCCACGCCTTCTTCAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGt                                                                                                                                        <  7:240283/149‑1 (MQ=255)
             cTTCCACCTGATGACCCACGCCTTCTTCAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCGTAGGCCAGcgg                                                                                                                                           >  7:280088/1‑139 (MQ=255)
             cTTCCACCTGATGACCCACGCCTTCTTCAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCGTAGGCCAGcgg                                                                                                                                           <  8:280088/141‑3 (MQ=255)
              ttCCACCTGATGACCCACGCCTTCTTCAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCgg                                                                                                                                  >  4:77961/1‑149 (MQ=255)
                ccACCTGATGACCCACGCCTTCTTCAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCggtg                                                                                                                                >  5:157267/1‑149 (MQ=255)
                  aCCTGATGACCCACGCCTTCTTCAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCggtggt                                                                                                                              <  2:18665/149‑1 (MQ=255)
                       atgcCCCCCGCCGTCCTCCACGCCCTGCTGTTCCTTGCCTCCGGTGCGGGCCCCGTTCCCTGACACCACGAGCAGAACATCTTCAAGATGCGCGGCCTGTGCAAGAAACAGCCGATCGCCTACGCCAGCGGCGTGGTCGGTGGTGCGGc                                                                                                                         <  8:152262/149‑1 (MQ=255)
                            ccACGCCTTCTTCAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTgg                                                                                                                    >  5:289591/1‑149 (MQ=255)
                             cACGCCTTCTTCAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTggc                                                                                                                   <  7:81938/149‑1 (MQ=255)
                               cGCCTTCTTCAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTggcgg                                                                                                                 <  1:88388/149‑1 (MQ=255)
                               cGCCTTCTTCAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTggcgg                                                                                                                 <  7:289077/149‑1 (MQ=255)
                                  cttcttCAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGccc                                                                                                              >  1:294026/1‑149 (MQ=255)
                                      ttcAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCggtgg                                                                                                                               >  3:220329/1‑128 (MQ=255)
                                      ttcAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGcc                                                                                                          >  3:350769/1‑149 (MQ=255)
                                      ttcAAGGCCCTGCTGTTCCTTGCCTCCGGGGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCggtgg                                                                                                                               <  4:220329/128‑1 (MQ=255)
                                        cAAGGCCCTGCTGTTCCTTGCCTCCGGTCCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGa                                                                                                        >  1:218596/1‑149 (MQ=255)
                                         aaGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGAt                                                                                                       >  3:268131/1‑149 (MQ=255)
                                          cggCCCTGCTGTTCCTTGCCTCCGGTGCGGGGATCGGTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATc                                                                                                      <  4:268131/148‑1 (MQ=255)
                                          aGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCCAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGGGCGGCACTGGCGGCCCTGCCGATc                                                                                                      <  7:246171/149‑1 (MQ=255)
                                          aGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATc                                                                                                      <  8:6468/149‑1 (MQ=255)
                                           ggCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCg                                                                                                     >  2:325638/1‑149 (MQ=255)
                                                 gctgTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCg                                                                                               >  6:376373/1‑149 (MQ=255)
                                                   tgTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTa                                                                                             <  7:309764/149‑1 (MQ=255)
                                                     ttCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGCTGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTa                                                                                             <  8:3325/147‑1 (MQ=255)
                                                     ttCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAg                                                                                            <  3:199377/148‑1 (MQ=255)
                                                            ccTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTa                                                                                     <  8:238621/148‑1 (MQ=255)
                                                            ccTCCGGTGCGGTGATCGTTGCCTGCCACAACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTAGGTGGTGCGGCACTCGCGGCCCTGCCGATCGTGACCGTAGGATTCTAc                                                                                    >  8:103429/1‑149 (MQ=255)
                                                               ccGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTcc                                                                                 <  1:375947/149‑1 (MQ=255)
                                                               ccGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTc                                                                                  >  1:330180/1‑148 (MQ=255)
                                                                 ggtgcggtgATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCaa                                                                               <  5:158111/149‑1 (MQ=255)
                                                                     cggtgATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCGTCGTGGGCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGa                                                                            <  2:349603/148‑1 (MQ=255)
                                                                       gtgATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGAc                                                                           >  2:345845/1‑147 (MQ=255)
                                                                              ttGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCctc                                                                  >  1:367317/1‑149 (MQ=255)
                                                                                 cctgccACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTgg                                                               >  6:185581/1‑149 (MQ=255)
                                                                                     ccaccaCGAGCATAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGaa                                                            <  7:219366/148‑1 (MQ=255)
                                                                                      caccacGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGc                                                          <  3:77961/149‑1 (MQ=255)
                                                                                        ccacGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCt                                                        >  7:63527/1‑149 (MQ=255)
                                                                                        ccacGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCt                                                        >  6:72521/1‑149 (MQ=255)
                                                                                          acGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCg                                                                                                                                         <  7:95973/66‑1 (MQ=255)
                                                                                          acGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCg                                                                                                                                         >  8:95973/1‑66 (MQ=255)
                                                                                          acGAGCAGAAAATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTc                                                      <  8:325864/149‑1 (MQ=255)
                                                                                            gAGCAGAACATCTTCAAGATGGGCGGCCTGTGGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGTTGCGGCACTGGCGCCCATGCCGATCGTGCCCGTAGGATTCTACTCCAAGGACGAGATCCTCTGGGaa                                                            >  1:238623/1‑142 (MQ=255)
                                                                                            gAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGGACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGaa                                                            <  2:238623/141‑1 (MQ=255)
                                                                                            gAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGc                                                    <  2:367317/149‑1 (MQ=255)
                                                                                                aGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGc                                                <  3:175117/149‑1 (MQ=255)
                                                                                                  aaCATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCgg                                              >  7:362446/1‑149 (MQ=255)
                                                                                                   aCATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCgg                                              <  7:103429/148‑1 (MQ=255)
                                                                                                       cTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAaca                                         >  2:299997/1‑149 (MQ=255)
                                                                                                        ttCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAacac                                        >  5:350409/1‑149 (MQ=255)
                                                                                                        ttCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAacac                                        >  4:318036/1‑149 (MQ=255)
                                                                                                         tCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACAcc                                       >  2:54675/1‑149 (MQ=255)
                                                                                                           aaGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCgg                                     >  5:277916/1‑149 (MQ=255)
                                                                                                            aGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGc                                    >  5:268024/1‑149 (MQ=255)
                                                                                                            aGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGc                                    <  1:325638/149‑1 (MQ=255)
                                                                                                             gATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGcc                                   >  4:290704/1‑149 (MQ=255)
                                                                                                                  gcggcCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCctgctg                              <  6:268024/149‑1 (MQ=255)
                                                                                                                          tgGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGcc                        <  6:157267/147‑1 (MQ=255)
                                                                                                                             gaagaaACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGcc                    >  6:325194/1‑149 (MQ=255)
                                                                                                                             gaagaaACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGc                     <  2:294026/148‑1 (MQ=255)
                                                                                                                                gaaACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCGCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCATGCTGTATGCCGGCCggg                 >  6:66439/1‑149 (MQ=255)
                                                                                                                                gaaACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTgg                 >  4:255311/1‑149 (MQ=255)
                                                                                                                                gaaAATGCCGCTTGCATACGCCAGATTCGTGCTCTGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTATGGTAAGCCTTCTCCAGCTGCAACACCGGCCTGCTGTATGCCGGCCTgg                 >  4:325894/1‑149 (MQ=255)
                                                                                                                                 aaaCTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTGGt                <  5:185581/149‑1 (MQ=255)
                                                                                                                                 aaaCTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTGGt                <  2:330180/149‑1 (MQ=255)
                                                                                                                                  aaCTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTGGt                <  8:63527/148‑1 (MQ=255)
                                                                                                                                   aCTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTGGTAg              <  8:362446/149‑1 (MQ=255)
                                                                                                                                   aCTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTGGTAg              <  1:54675/149‑1 (MQ=255)
                                                                                                                                    cTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTGGTAgg             <  6:289591/149‑1 (MQ=255)
                                                                                                                                     tGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAaca                                         >  5:171315/1‑120 (MQ=255)
                                                                                                                                     tGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAaca                                         <  6:171315/120‑1 (MQ=255)
                                                                                                                                      gaagCTGGCGTAGGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTGGTAGgcg           >  8:64903/4‑149 (MQ=255)
                                                                                                                                      gaagCTGGCGTAGGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTGGTAGgcg           >  5:20146/4‑149 (MQ=255)
                                                                                                                                       ccGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTGGTAGgcgc          >  1:76087/1‑149 (MQ=255)
                                                                                                                                        cGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTGGTAGGcgcg         <  1:299997/149‑1 (MQ=255)
                                                                                                                                         gCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTGGTAGGCGCGt        >  7:127953/1‑149 (MQ=255)
                                                                                                                                             gCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTGGTAGGCGCGTTc      <  7:334879/147‑1 (MQ=255)
                                                                                                                                             gCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTGGTAGGCGCGTTCAt    <  4:350769/149‑1 (MQ=255)
                                                                                                                                              ccTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTGGTAGGCGCGTTCAt    <  5:131799/148‑1 (MQ=255)
                                                                                                                                              ccTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTGGTAGGCGCGTTCAt    >  6:131799/1‑148 (MQ=255)
                                                                                                                                               cTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTGGTAGGCGCGTTCATGa  <  8:127953/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GGGACGCGGCGATCTTCCACCTGATGACCCACGCCTTCTTCAAGGCCCTGCTGTTCCTTGCCTCCGGTGCGGTGATCGTTGCCTGCCACCACGAGCAGAACATCTTCAAGATGGGCGGCCTGTGGAAGAAACTGCCGCTGGCCTACGCCAGCTTCGTGGTCGGTGGTGCGGCACTGGCGGCCCTGCCGATCGTGACCGTAGGCTTCTACTCCAAGGACGAGATCCTCTGGGAAGCCTTCGCCAGCGGCAACACCGGCCTGCTGTATGCCGGCCTGGTAGGCGCGTTCATGA  >  NC_002947/4667223‑4667513

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: