Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,509,940 T→C 6.3% Q41R (CAG→CGG)  serS ← serine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,509,9400TC 6.3% 186.3 / 2.4 63Q41R (CAG→CGG) serSserine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (29/30);  new base C (2/2);  total (31/32)
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.81e-01

CAGCTCAGCCTTGCCGGCAGCCAGTTCGTTGGCCATGCGCTCGACGTCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCACTTCCTGAAGTTGGCCGCGTAACAGTTTGGAATCGAGCATATCGTTCTCTCGTTATATG  >  NC_002947/4509795‑4510080
                                                                                                                                                  |                                                                                                                                             
cAGCTCAGCCTTGCCGGCAGCCAGTTCGTTGGCCATGCGCTCGACGTCAGCGCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTg                                                                                                                                              >  1:237249/1‑148 (MQ=255)
   cgcAGCCTTGCCGGCAGCCAGTTCGTTGGCCATGCGCTCGACGGCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCg                                                                                                                                         <  8:134985/147‑1 (MQ=255)
      aGCCTTGCCGGCAGCCAGTTCGTTGGCCATGCGCTCGACGTCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGcc                                                                                                                                       <  2:237249/148‑1 (MQ=255)
              cGGCAGCCAGTTCGTTGGCCATGCGCTCGACGTCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGgcgc                                                                                                                              >  7:163551/1‑149 (MQ=255)
              cGGCAGCCAGTTCGTTGGCCATGCGCTCGACGTCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTCGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGgcgc                                                                                                                              <  5:42390/149‑1 (MQ=255)
                     cAGTTCGTTGGCCATGCGCTCGACGTCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCa                                                                                                                       <  5:95337/149‑1 (MQ=255)
                      aGTTCGTTGGCCATGCGCTCGACGTCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAg                                                                                                                      <  6:320307/149‑1 (MQ=255)
                      aGTTCGTTGGCCATGCGCTCGACGTCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAg                                                                                                                      >  5:320307/1‑149 (MQ=255)
                       gTTCGTTGGCCATGCGCTCGACGTCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGt                                                                                                                     >  1:16708/1‑149 (MQ=255)
                           gTTGGCCATGCGCTCGACGTCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGAtt                                                                                                                 <  4:97841/149‑1 (MQ=255)
                           gTTGGCCATGCGCTCGACGTCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGAtt                                                                                                                 <  3:115171/149‑1 (MQ=255)
                            ttGGCCATGCGCTCGACGTCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCAGACCCGCACCgaacaa                                                                                                                                    <  8:282869/129‑6 (MQ=255)
                            ttGGCCATGCGCTCGACGTCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCAGACCCGCACCgaacaa                                                                                                                                    >  7:282869/1‑124 (MQ=255)
                                  aTGCGCTCGACGTCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATgcg                                                                                                          >  2:266554/1‑149 (MQ=255)
                                       cTCGACGTCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGa                                                                                                     <  7:63991/149‑1 (MQ=255)
                                       cTCGACGTCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGa                                                                                                     <  4:275192/149‑1 (MQ=255)
                                              tCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCAGACCCGCACCga                                                                                                                                        >  1:286457/1‑106 (MQ=255)
                                              tCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCAGACCCGCACCga                                                                                                                                        <  2:286457/107‑2 (MQ=255)
                                                aGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGc                                                                                           <  8:153153/149‑1 (MQ=255)
                                                  ccATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGTACCGCCTTGCGGCGCTCTTCCAGTGATTCGAt                                                                                                             >  1:186360/1‑130 (MQ=255)
                                                  ccATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGTACCGCCTTGCGGCGCTCTTCCAGTGATTCGAt                                                                                                             <  2:186360/130‑1 (MQ=255)
                                                  ccATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTg                                                                                         <  3:159317/149‑1 (MQ=255)
                                                          cGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAg                                                                                             >  3:298584/1‑138 (MQ=255)
                                                          cGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAg                                                                                             <  4:298584/138‑1 (MQ=255)
                                                            gggCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACgg                                                                                <  1:266554/149‑1 (MQ=255)
                                                             ggCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACgg                                                                                >  6:172854/1‑148 (MQ=255)
                                                                   gTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGc                                                                           >  7:197236/1‑147 (MQ=255)
                                                                      ttCGNCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGc                                                                      >  5:32137/1‑149 (MQ=255)
                                                                      ttCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTc                                                                                                                         >  4:247578/1‑99 (MQ=255)
                                                                      ttCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTc                                                                                                                         <  3:247578/99‑1 (MQ=255)
                                                                      ttCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTg                                                                                                                    >  3:6758/1‑104 (MQ=255)
                                                                      ttCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTg                                                                                                                    <  4:6758/104‑1 (MQ=255)
                                                                      ttCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGc                                                                      >  7:123678/1‑149 (MQ=255)
                                                                      ttCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGc                                                                      >  1:15134/1‑149 (MQ=255)
                                                                           ccttggccttggCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGct                                                                                                                             <  1:246322/90‑1 (MQ=255)
                                                                           ccttggccttggCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGct                                                                                                                             >  2:246322/1‑90 (MQ=255)
                                                                           ccttggccttggCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGAtt                                                                                                                 >  5:268955/1‑102 (MQ=255)
                                                                           ccttggccttggCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGAtt                                                                                                                 <  6:268955/102‑1 (MQ=255)
                                                                           ccttggccttggCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTcc                                                                 <  8:123678/149‑1 (MQ=255)
                                                                           ccttggccttggCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTcc                                                                 <  2:15134/149‑1 (MQ=255)
                                                                             ttggccttggccTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGa                                                                              <  5:72777/134‑1 (MQ=255)
                                                                             ttggccttggccTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGa                                                                              >  6:72777/1‑134 (MQ=255)
                                                                               ggccttggcctGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGc                                                               <  8:197236/147‑1 (MQ=255)
                                                                                     ggcctggccGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCCCTTCCt                                                       >  1:262255/1‑149 (MQ=255)
                                                                                           gccGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAg                                                                             <  3:15310/121‑1 (MQ=255)
                                                                                           gccGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAg                                                                             >  4:15310/1‑121 (MQ=255)
                                                                                               aTGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCACTTCCTGAAGTTGGcc                                             >  4:34984/1‑149 (MQ=255)
                                                                                                     ttGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCACTTCCTGAAGTTGGCCGCGTaa                                       <  3:34984/149‑1 (MQ=255)
                                                                                                     ttGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCACTTCCTGAAGTTGGCCGCGTaa                                       >  6:92578/1‑149 (MQ=255)
                                                                                                       ggAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCACTTCCTGAAGTTGGCCGCGTaa                                       >  2:200110/1‑147 (MQ=255)
                                                                                                           cGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCACTTCCTGAAGTTGGCCGCGTAACAGttt                                 <  6:32137/149‑1 (MQ=255)
                                                                                                           cGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCACTTCCTGAAGTTGGCCGCGTAACAGttt                                 <  8:163551/149‑1 (MQ=255)
                                                                                                              gCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCACTTCCTGAAGTTGGCCGCGTAACAGTTTgg                               >  4:39081/1‑148 (MQ=255)
                                                                                                                     gCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCACTTCCTGAAGTTGGCCGCGTAACAGTTTGGAATCGAg                        >  5:22828/1‑148 (MQ=255)
                                                                                                                     gCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCACTTCCTGAAGTTGGCCGCGTAACAGTTTGGAATCGAGc                       <  1:200110/149‑1 (MQ=255)
                                                                                                                       tCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCACTTCCTGAAGTTGGCCGCGTAACAGTTTGGAATCGAGCa                      >  7:24265/1‑148 (MQ=255)
                                                                                                                          gCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCACTTCCTGAAGTTGGCCGCGTAACAGTTTGGAATCGAGCATATc                  <  5:172854/149‑1 (MQ=255)
                                                                                                                            cTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGCCCAGGTGTCCGCCACTTCCTGAAGTTGGCAGCGTAACAGTTTGGAATCGAGCATATCttt               >  2:284588/1‑149 (MQ=255)
                                                                                                                             tGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCACTTCCTGAAGTTGGCCGCGTAACAGTTTGGAATCGAGCATATCGt                >  3:107969/1‑148 (MQ=255)
                                                                                                                               cAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCACTTCCTGAAGTTGGCCGCGTAACAGTTTGGAATCGAGCATATCGTtct             <  8:24265/149‑1 (MQ=255)
                                                                                                                                     ttCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCACTTCCTGAAGTTGGCCGCGTAACAGttt                                 <  4:120152/123‑1 (MQ=255)
                                                                                                                                     ttCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCACTTCCTGAAGTTGGCCGCGTAACAGttt                                 >  3:120152/1‑123 (MQ=255)
                                                                                                                                             gggTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATGCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCACTTCCTGAAGTTGGCCGCGTAACAGTTTGGAATCGAGCATATCGTTCTCTCGTTATATg  >  1:78984/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CAGCTCAGCCTTGCCGGCAGCCAGTTCGTTGGCCATGCGCTCGACGTCAGCCATCAGCGGGGCGATGTCTTCGCCCTTGGCCTTGGCCTGGCCGATGGACTTGGAACGGGCGTTACGCTCAGCCTGCAGTTGTTCGGTGCGGGTCTGCACCGCCTTGCGGCGCTCTTCCAGTGATTCGATGCGCGCGACATCCAGGCTGAAGCCACGGGAGGCCAGGCGGTCCGCCACTTCCTGAAGTTGGCCGCGTAACAGTTTGGAATCGAGCATATCGTTCTCTCGTTATATG  >  NC_002947/4509795‑4510080

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: