Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,000,803 C→G 6.2% V65L (GTG→CTG)  PP_2624 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,000,8030CG 6.2% 181.3 / 2.4 63V65L (GTG→CTG) PP_2624hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/29);  new base G (2/2);  total (32/31)
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.62e-01

AGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGA  >  NC_002947/3000658‑3000949
                                                                                                                                                 |                                                                                                                                                  
aGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCt                                                                                                                                                 <  2:149320/149‑1 (MQ=255)
     gggCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTg                                                                                                                                             >  5:348652/1‑148 (MQ=255)
       gCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAat                                                                                                                                          <  6:262358/149‑1 (MQ=255)
         aCCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGAACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACAGCGTGTTGaacga                                                                                                                                        >  7:24841/1‑146 (MQ=255)
          ccTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGAACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACAGCGTGTTGaacgac                                                                                                                                       <  8:24841/149‑5 (MQ=255)
                 cgcgATCCGGTCCGGCCTGAATTCGTTGATGGTGCTGAACTGCAGGTGGACGGCCTGGTCTGTATCAGGGGCGGTGCCGAGCCTGGACGGTACCGTCATGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTAt                                                                                                                                <  1:88591/149‑1 (MQ=255)
                  gcgATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGc                                                                                                                                                  >  1:276260/1‑130 (MQ=255)
                  gcgATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGc                                                                                                                                                  >  5:30759/1‑130 (MQ=255)
                  gcgATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGc                                                                                                                                                  <  2:276260/130‑1 (MQ=255)
                  gcgATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGc                                                                                                                                                  <  6:30759/130‑1 (MQ=255)
                   cgATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATc                                                                                                                               >  1:305387/1‑148 (MQ=255)
                     aTCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGaa                                                                                                                                           >  8:97616/1‑134 (MQ=255)
                     aTCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGaa                                                                                                                                           <  7:97616/134‑1 (MQ=255)
                     aTCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGAt                                                                                                                            >  8:85904/1‑149 (MQ=255)
                               gCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGAcc                                                                                                                  >  7:9983/1‑149 (MQ=255)
                                  cGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTGATAGATGCGcctgacctga                                                                                                               >  4:57158/1‑149 (MQ=255)
                                  cGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCCTCACCTCCACGCACACATTGGCGTCGTACAACACGCTGTTGAAAGTATGTTTATCGATCCGCATTAcatta                                                                                                               >  4:61875/1‑145 (MQ=255)
                                   gAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACAGCGTGTTGaac                                                                                                                                          >  8:174585/1‑120 (MQ=255)
                                   gAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACAGCGTGTTGaac                                                                                                                                          <  7:174585/121‑2 (MQ=255)
                                   gAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGAc                                                                                                              <  7:85904/149‑1 (MQ=255)
                                       tCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTc                                                                                                          <  1:149281/149‑1 (MQ=255)
                                         gTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGAc                                                                                                                   >  5:117357/1‑138 (MQ=255)
                                         gTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGAc                                                                                                                   <  6:117357/138‑1 (MQ=255)
                                          ttGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGtt                                                                                                                                              >  3:276487/1‑110 (MQ=255)
                                          ttGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGtt                                                                                                                                              <  4:276487/110‑1 (MQ=255)
                                              tGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAgc                                                                                                   <  1:341184/149‑1 (MQ=255)
                                                      aaCTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCAc                                                                                               >  6:257127/1‑145 (MQ=255)
                                                      aaCTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCAc                                                                                               <  5:257127/145‑1 (MQ=255)
                                                      aaCTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACgg                                                                                             <  2:305387/147‑1 (MQ=255)
                                                      aaCTGCAGGTGGACGACCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatgta                                                                                                                                       >  7:323753/1‑105 (MQ=255)
                                                      aaCTGCAGGTGGACGACCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatgta                                                                                                                                       <  8:323753/105‑1 (MQ=255)
                                                         tGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGt                                                                                         >  2:206745/1‑148 (MQ=255)
                                                           cAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGGGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCa                                                                                       >  6:47147/1‑148 (MQ=255)
                                                            aGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATc                                                                                     >  1:150255/1‑149 (MQ=255)
                                                              gTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGaa                                                                                                                                           <  1:263952/93‑1 (MQ=255)
                                                              gTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGaa                                                                                                                                           >  2:263952/1‑93 (MQ=255)
                                                                ggACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCt                                                                                           >  5:200186/1‑139 (MQ=255)
                                                                ggACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCt                                                                                           <  6:200186/139‑1 (MQ=255)
                                                                        aggtCGGTATGCGGGGCGGTGCCGAGCCTGGACGGTACCTTCTCTTGCAGGCTCACATTGGCGTCGTTCAACTCGCTGTTGTATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTc                                                                         <  3:61875/146‑1 (MQ=25)
                                                                        aGCTCGGTATCCGGGGCGGTGCCGCGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGGTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTc                                                                         <  3:57158/149‑1 (MQ=255)
                                                                                    gtgtcGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGccc                                                                     <  5:59249/137‑1 (MQ=255)
                                                                                    ggggCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCccagcaccag                                                             >  2:307712/1‑149 (MQ=255)
                                                                                    ggggCGGTGCCGAGCCAGGACGGTAACGTCACGTCCAGGCTCACATTGGCGTCGTTAAACACGCTGTTGAATGTATGTTTATAGATTCGCATTACCATACCGTCCACAACCGCACGGCTGTCATACAGGCCGAACTCACccaataacg                                                              >  6:59249/1‑142 (MQ=255)
                                                                                       gCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCAt                                                          <  8:9983/149‑1 (MQ=255)
                                                                                             ccGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTtgtg                                                    >  8:254847/1‑149 (MQ=255)
                                                                                                  ccTGGACGTGACGGTACCGTCCAGGCTCGCATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCgg                                                 >  1:318929/1‑147 (MQ=18)
                                                                                                   cTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAg                                              <  2:150255/149‑1 (MQ=255)
                                                                                                      gACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGc                                             <  5:47147/147‑1 (MQ=255)
                                                                                                         ggTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTAt                                                                                                                                >  3:255471/1‑61 (MQ=255)
                                                                                                         ggTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTAt                                                                                                                                <  4:255471/61‑1 (MQ=255)
                                                                                                          gTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCAcc                                       <  7:254847/149‑1 (MQ=255)
                                                                                                                  aCGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTg                               <  6:348652/149‑1 (MQ=255)
                                                                                                                     tCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCAGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGca                             >  3:19686/1‑147 (MQ=255)
                                                                                                                            tCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCt                      >  3:265856/1‑148 (MQ=255)
                                                                                                                            tCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTg                     >  3:188870/1‑149 (MQ=255)
                                                                                                                            tCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTg                     >  7:176719/1‑149 (MQ=255)
                                                                                                                              acaTTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCt                      <  4:265856/146‑1 (MQ=255)
                                                                                                                               caTTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTgg                    >  4:141395/1‑147 (MQ=255)
                                                                                                                                      gtcgtTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTGAAGCGCACTGCTGCCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTTGTAGCgggta           <  7:232595/149‑3 (MQ=255)
                                                                                                                                      gtcgaTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAATCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGAATCTTGTGCGGCAGCTCCACCTCAGACTGCTCATCGCTGGTAGCgggg            >  8:232595/1‑148 (MQ=255)
                                                                                                                                         gtTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACa         >  7:124014/1‑148 (MQ=255)
                                                                                                                                         gtTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACAt        <  1:206745/149‑1 (MQ=255)
                                                                                                                                             aaCACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCGAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACtt    <  8:124014/149‑1 (MQ=255)
                                                                                                                                               cacGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGa  <  8:176719/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGA  >  NC_002947/3000658‑3000949

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: