Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,063,001 C→G 9.4% E446D (GAG→GAC) ‡ PP_3581 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,063,0010CG 9.4% 161.9 / 7.3 53E446D (GAG→GAC) ‡PP_3581membrane protein
Reads supporting (aligned to +/- strand):  ref base C (22/26);  new base G (3/2);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

GCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGCTG  >  NC_002947/4062853‑4063147
                                                                                                                                                    |                                                                                                                                                  
gCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCgctgc                                                                                                                                                    >  3:25617/1‑149 (MQ=255)
 cccGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCgctgct                                                                                                                                                   >  4:69195/1‑149 (MQ=255)
 cccGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCgctgct                                                                                                                                                   >  8:192594/1‑149 (MQ=255)
  ccGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTc                                                                                                                                                  <  7:192594/149‑1 (MQ=255)
   cGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTc                                                                                                                                                  >  5:289724/1‑148 (MQ=255)
      tCTTTTCCGCCAGGCCGATGAAGCGGTCGGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAgc                                                                                                                                              <  7:142203/149‑1 (MQ=255)
        ttttCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAgcgc                                                                                                                                            >  7:105648/1‑149 (MQ=255)
              gCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGCCa                                                                                                                                          <  6:100974/145‑1 (MQ=255)
              gCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGCCa                                                                                                                                          >  5:100974/1‑145 (MQ=255)
              gCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGgcg                                                                                                                                      >  4:132142/1‑149 (MQ=255)
               ccAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGgcgc                                                                                                                                     <  6:324495/149‑1 (MQ=255)
                  ggCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAAGTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAg                                                                                                                                  <  2:253835/149‑1 (MQ=255)
                       aTGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTg                                                                                                                             <  3:132142/149‑1 (MQ=255)
                          aaGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTgca                                                                                                                           <  6:289724/148‑1 (MQ=255)
                             cGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGcagcag                                                                                                                       >  5:207289/1‑149 (MQ=255)
                                    gcagcagcaGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAAc                                                                                                                >  2:88095/1‑149 (MQ=255)
                                    gcagcagcaGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAAc                                                                                                                >  7:134949/1‑149 (MQ=255)
                                    cccgcagcagGGCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAAc                                                                                                                >  5:315798/4‑149 (MQ=255)
                                      agcagcagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAActg                                                                                                              <  8:134949/149‑1 (MQ=255)
                                         agcagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAActgctg                                                                                                           >  4:73956/1‑149 (MQ=255)
                                           cagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGgcg                                                                                                                                      <  5:107256/120‑1 (MQ=255)
                                           cagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGgcg                                                                                                                                      >  6:107256/1‑120 (MQ=255)
                                           cagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGgc                                                                                                         <  7:294709/149‑1 (MQ=255)
                                             ggCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGgcgc                                                                                                       <  6:207289/149‑1 (MQ=255)
                                              gCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGgcgcg                                                                                                      <  6:256898/149‑1 (MQ=255)
                                               cccTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGgcgcgc                                                                                                     <  8:105648/149‑1 (MQ=255)
                                                              tCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGCCa                                                                                                                                          >  4:177597/1‑97 (MQ=255)
                                                              tCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGCCa                                                                                                                                          <  3:177597/97‑1 (MQ=255)
                                                                ccAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGaca                                                                                    <  3:73956/149‑1 (MQ=255)
                                                                  aGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATc                                                                                        >  4:35072/1‑143 (MQ=255)
                                                                  aGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATc                                                                                        <  3:35072/143‑1 (MQ=255)
                                                                   gcgcAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACAt                                                                                 >  5:78323/1‑149 (MQ=255)
                                                                    cgcAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATc                                                                                <  3:69195/149‑1 (MQ=255)
                                                                      cagcaACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCga                                                                              <  6:315798/149‑1 (MQ=255)
                                                                       agcaACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCgaa                                                                             <  6:155704/149‑1 (MQ=255)
                                                                          aaCGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCgaag                                                                            >  5:25944/1‑147 (MQ=255)
                                                                           aCGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGaagaag                                                                         <  4:25617/149‑1 (MQ=255)
                                                                                 cGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGTAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGt                                                                     <  1:88095/147‑1 (MQ=255)
                                                                                         aTCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGGAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGGGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAAGCGTTCtt                                                           <  4:104761/149‑1 (MQ=255)
                                                                                             gTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCg                                                       >  8:47651/1‑149 (MQ=255)
                                                                                                        tGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCAt                                                                                                   >  4:313313/1‑94 (MQ=255)
                                                                                                        tGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCAt                                                                                                   <  3:313313/94‑1 (MQ=255)
                                                                                                            cGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGt                                                                                             <  7:151591/96‑1 (MQ=255)
                                                                                                            cGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGt                                                                                             >  8:151591/1‑96 (MQ=255)
                                                                                                            cGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATgg                                        <  7:47651/149‑1 (MQ=255)
                                                                                                             gTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGt                                       <  6:25944/149‑1 (MQ=255)
                                                                                                              tCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGt                                       >  1:1693/1‑148 (MQ=255)
                                                                                                                  aCTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATg                                  <  6:78323/149‑1 (MQ=255)
                                                                                                                            gTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCgg                        >  4:208604/1‑149 (MQ=255)
                                                                                                                                 gCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTa                   >  3:284662/1‑149 (MQ=255)
                                                                                                                                 gCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTa                   >  2:70482/1‑149 (MQ=255)
                                                                                                                                  cAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTa                   <  2:1693/148‑1 (MQ=255)
                                                                                                                                    aGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACg                 >  4:279752/1‑148 (MQ=255)
                                                                                                                                    aGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACg                 >  7:8829/1‑148 (MQ=255)
                                                                                                                                           cTGGCGCTGGAGCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCt         >  5:143619/1‑149 (MQ=255)
                                                                                                                                             ggCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCg       >  2:208899/1‑149 (MQ=255)
                                                                                                                                                gctgctCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGc    >  2:53587/1‑149 (MQ=255)
                                                                                                                                                 ctgctCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGCt   >  8:139936/1‑149 (MQ=255)
                                                                                                                                                    ctCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGCTg  <  3:208604/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGCTG  >  NC_002947/4062853‑4063147

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