Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,588,798 G→T 6.6% T721T (ACG→ACT PP_2266 → DNA‑directed RNA polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,588,7980GT 6.6% 186.4 / 3.8 61T721T (ACG→ACTPP_2266DNA‑directed RNA polymerase
Reads supporting (aligned to +/- strand):  ref base G (31/26);  new base T (2/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.63e-01

GTGAAAGCTGGTGAGGCCATGAAGTGGCTCCAGGACGCTGCAACCCTGGCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCGATGCGGCCCACCTCATGCTCACCGTTGTGCGTGGCGTGGA  >  NC_002947/2588655‑2588944
                                                                                                                                               |                                                                                                                                                    
gTGAAAGCTGGTGAGGCCATGAAGTGGCTCCAGGACGCTGCAACCCTGGCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCAt                                                                                                                                                 >  4:186074/1‑149 (MQ=255)
gTGAAAGCTGGTGAGGCCATGAAGTGGCTCCAGGACGCTGCAACCCTGGCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCAt                                                                                                                                                 >  6:249727/1‑149 (MQ=255)
  gAAAGCTGGTGAGGCCATGAAGTGGCTCCAGGACGCTGCACCCCTGGCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATc                                                                                                                                                >  5:231358/1‑148 (MQ=255)
  gAAAGCTGGTGAGGCCATGAAGTGGCTCCAGGACGCTGCAACCCTGGCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATc                                                                                                                                                >  7:249307/1‑148 (MQ=255)
       cTGGTGAGGCCATGAAGTGGCTCCAGGACGCTGCAACCCTGGCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACggg                                                                                                                                         >  5:203699/1‑149 (MQ=255)
         ggTGAGGCCATGAAGTGGCTCCAGGACGCTGCAACCCTGGCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGc                                                                                                                                       >  5:30578/1‑149 (MQ=255)
         ggTGAGGCCATGAACTGGCTCCAGGACGCTGCAACCCTGGCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGc                                                                                                                                       >  8:253201/1‑149 (MQ=255)
          gTGAGGCCATGAAGTGGCTCCAGGACGCTGCAACCCTGGCCGCAAAAGAAGAGCTGCAGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTTATGCACGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGccc                                                                                                                                      >  5:44186/1‑149 (MQ=255)
           tGAGGCCATGAAGTGGCTCCAGGACGCTGCAACCCTGGCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACCCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGccc                                                                                                                                     <  6:117980/149‑1 (MQ=255)
             aGGCCATGAAGTGGCTCCAGGACGCTGCAACCCTGGCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCgt                                                                                                                                   <  6:226083/149‑1 (MQ=255)
                 cATGAAGTGGCTCCAGGACGCTGCAACCCTGGCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTgc                                                                                                                               <  6:248168/149‑1 (MQ=255)
                                  aCGCTGCAACCCTGGCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGaa                                                                                                             <  3:186074/149‑1 (MQ=255)
                                    gCTGCAACCCTGGCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAaga                                                                                                           <  6:30578/149‑1 (MQ=255)
                                      tGCAACCCTGGCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACTGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAagag                                                                                                          >  2:38899/1‑148 (MQ=255)
                                      tGCAACCCTGGCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACTGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAagag                                                                                                          >  2:92136/1‑148 (MQ=255)
                                        cAACCCTGGCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCaa                                                                                                       <  7:253201/149‑1 (MQ=255)
                                                gCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCt                                                                                               >  4:337447/1‑149 (MQ=255)
                                                gCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCt                                                                                               >  1:169528/1‑149 (MQ=255)
                                                   gCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGa                                                                                            <  5:223010/149‑1 (MQ=255)
                                                      aaaGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACTGCCATCAACGGGGCCCGTGTGCGCCTTCTGATGAACGAAGAGCAAGACGCCCTGGACaa                                                                                         <  1:38899/149‑1 (MQ=255)
                                                       aagaagAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACTGCCATCAACCGGGGCCCGTGTGCTCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACaa                                                                                         <  1:92136/149‑1 (MQ=255)
                                                            gagCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAgcgc                                                                                     >  1:140206/1‑147 (MQ=255)
                                                             agCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAg                                                                                  <  3:296792/149‑1 (MQ=255)
                                                                  ccGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTc                                                                             >  4:88420/1‑149 (MQ=255)
                                                                        tgCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAgg                                                                        <  5:330888/148‑1 (MQ=255)
                                                                           tGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCGGTCGCAGGGGAt                                                                    >  4:242754/1‑149 (MQ=255)
                                                                             gACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGAt                                                                    >  8:180408/1‑147 (MQ=255)
                                                                                 aCTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGccc                                                              <  3:337447/149‑1 (MQ=255)
                                                                                 aCTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGccc                                                              <  2:140206/149‑1 (MQ=255)
                                                                                        tGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCGACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTc                                                       >  5:258307/1‑149 (MQ=255)
                                                                                          ggCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCTGTCGCAGGGGATCAGCCCCAACTTCg                                                      <  3:242754/148‑1 (MQ=255)
                                                                                          ggCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGt                                                     <  5:249727/149‑1 (MQ=255)
                                                                                           gCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCaa                                                                                   >  5:313478/1‑118 (MQ=255)
                                                                                           gCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCaa                                                                                   <  6:313478/118‑1 (MQ=255)
                                                                                           gCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTc                                                    <  5:173827/149‑1 (MQ=255)
                                                                                           gCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTc                                                    >  6:187631/1‑149 (MQ=255)
                                                                                             ttCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCa                                                  >  5:239623/1‑149 (MQ=255)
                                                                                               cccGGTGATGCAGGCGTTCAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGAGTGCGCCACCTGATGAACGAagatc                                                                                                         >  8:352720/1‑90 (MQ=255)
                                                                                                ccGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAg                                                                                        >  4:87451/1‑108 (MQ=255)
                                                                                                ccGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAg                                                                                        <  3:87451/108‑1 (MQ=255)
                                                                                                      aTGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCg                                          <  5:84619/148‑1 (MQ=255)
                                                                                                        gCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCGATg                                       >  2:16825/1‑149 (MQ=255)
                                                                                                         cAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCGATg                                       <  8:249307/148‑1 (MQ=255)
                                                                                                         cAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCGATGc                                      <  6:203699/149‑1 (MQ=255)
                                                                                                          aGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCGATGCg                                     >  3:353239/1‑149 (MQ=255)
                                                                                                           ggCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCGATGCgg                                    <  1:16825/149‑1 (MQ=255)
                                                                                                            gCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCt                                             <  8:187817/139‑1 (MQ=255)
                                                                                                            gCGTAAAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCt                                             >  7:187817/1‑139 (MQ=255)
                                                                                                                aCAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTACATGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCGATGCGGCCCa                                >  1:105837/1‑149 (MQ=255)
                                                                                                                 catcgccGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCGATGCGGCCCAc                               <  6:258307/143‑1 (MQ=255)
                                                                                                                    cTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCGATGCGGCCCACCTCa                           <  7:180408/149‑1 (MQ=255)
                                                                                                                     tCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGGCTACAGATAAACGAAGAGCAACACGCAGTGGACAAGCGAAAGCAGGCGCAGGGGATCAGCCCACAGGTCGGCAATATCAGCGGGGCGGCCCCCCTc                            >  2:329475/1‑147 (MQ=255)
                                                                                                                      ccGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCGATGCGGCCCACCTCATg                         >  7:108776/1‑149 (MQ=255)
                                                                                                                       cGTGTTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGCGCGCCTCCTGATGTACGAAGAGCAAGACGCCCTTGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCGATGCGGCCCACCTCATg                         <  2:313240/148‑1 (MQ=255)
                                                                                                                          ggTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCGATGCGGCCCACCTCATGCTCa                     >  2:28360/1‑149 (MQ=255)
                                                                                                                                   ccgcGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAGGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCGATGCGGCCCACCTCATGCTCACCGTTGTgc            >  7:364/1‑149 (MQ=255)
                                                                                                                                           aaaCGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGATCAAGACGCCCTGGACAAGCGAAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCGATGCGGCCCACCTCATGCTCACCGTTGTgcgtggcgtg    <  6:44186/149‑1 (MQ=255)
                                                                                                                                           aaaCGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCGATGCGGCCCACCTCATGCTCACCGTTGTgcgtggcgtg    >  7:48283/1‑149 (MQ=255)
                                                                                                                                           aaaCGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCGATGCGGCCCACCTCATGCTCACCGTTGTgcgtggcgtg    >  5:9498/1‑149 (MQ=255)
                                                                                                                                             aCGGCCATCAACGGGGGAAGTGTGCGCATCCTGATGAAAGAAGATCAAGACGCCGTGGTCAAGCGAAAGCAGTCGATGGGTATCAGCCCGAAGTTCGTGCATAGCTGAGATGCTGCACACCTCATGCTCACCGTTGTGCGTTGCGTGGa  <  8:364/149‑1 (MQ=255)
                                                                                                                                             aCGGCCATCAACGGGGCCAGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCGATGCGGCCCACCTCATGCTCAGCGTTGTgcgtggcgtg    >  2:189493/1‑147 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GTGAAAGCTGGTGAGGCCATGAAGTGGCTCCAGGACGCTGCAACCCTGGCCGCAAAAGAAGAGCTGCCGGTGTGCTGGACCACTGCCGTGGGCTTCCCGGTGATGCAGGCGTACAACTCCGTGGTAGCTCGCCGCGTCAAAACGGCCATCAACGGGGCCCGTGTGCGCCTCCTGATGAACGAAGAGCAAGACGCCCTGGACAAGCGCAAGCAGTCGCAGGGGATCAGCCCCAACTTCGTCCATAGCTGCGATGCGGCCCACCTCATGCTCACCGTTGTGCGTGGCGTGGA  >  NC_002947/2588655‑2588944

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: