Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 969,604 G→A 6.1% A173T (GCC→ACC)  queA → S‑adenosylmethionine:tRNA ribosyltransferase‑isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947969,6040GA 6.1% 196.6 / 2.2 65A173T (GCC→ACC) queAS‑adenosylmethionine:tRNA ribosyltransferase‑isomerase
Reads supporting (aligned to +/- strand):  ref base G (33/28);  new base A (2/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

TTCGCTTCACCGAAGAGGTGCTGCCGTTGCTCGACCGTGTCGGCCACATGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTCGACAAGATCGAAGACCACC  >  NC_002947/969467‑969745
                                                                                                                                         |                                                                                                                                             
ttCGCTTCACCGAAGAGGTGCTGCCGTTGCTCGACCGTGTCGGCCACATGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTAcc                                                                                                                                    >  8:124952/1‑149 (MQ=255)
ttCGCTTCACCGAAGAGGTGCTGCCGTTGCTCGACCGTGTCGGCCACATGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGCAACCGCGccg                                                                                                                                       >  5:364472/1‑145 (MQ=255)
ttCGCTTCACCGAAGAGGTGCTGCCGTTGCTCGACCGTGTCGGCCACATGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGCAACCGCGccg                                                                                                                                       <  6:364472/146‑2 (MQ=255)
    cTTCACCGAAGAGGTGCTGCCGTTGCTCGACCGTGTCGGCCACATGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGCCGAGCGTGAGCGCTACCAGACTGTGTCCGCCCCGCGCGCCGGCGCGGTTGCAGCGCGTCCCGCAg                                                                                                                                >  4:65103/1‑149 (MQ=255)
       cACCGAAGAGGTGCTGCCGTTGCTCGACCGTGTCGGCCACATGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAgg                                                                                                                               <  2:190730/147‑1 (MQ=255)
          cGAAGAGGTGCTGCCGTTGCTCGACCGTGTCGGCCACATGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGACGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGc                                                                                                                          >  3:302791/1‑149 (MQ=255)
                    ctgcCGTTGCTCGACCGTGTCGGCCACATGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTcgacg                                                                                                                 <  4:155390/148‑1 (MQ=255)
                       ccGTTGCTCGACCGTGTCGGCCACATGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAg                                                                                                               <  3:212733/147‑1 (MQ=255)
                         gTTGCTCGACCGTGTCGGCCACATGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGgcgc                                                                                                           >  6:68907/1‑149 (MQ=255)
                                gACCGTGTCGGCCACATGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCgctgc                                                                                                        <  4:238520/145‑1 (MQ=255)
                                gACCGTGTCGGCCACATGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCgctgc                                                                                                        >  3:238520/1‑145 (MQ=255)
                                  ccGTGTCGGCCACATGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGaaa                                                                                                  <  1:361715/149‑1 (MQ=255)
                                     tggcGGCCACATGCCGCTGCCCCCCTTTCTCGCCCGCCCGGCCGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGa                                                                                               <  4:192440/146‑1 (MQ=255)
                                       tCGGCCACATGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATc                                                                                             <  1:248157/149‑1 (MQ=255)
                                         ggCCACATGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGc                                                                                           <  8:188172/149‑1 (MQ=255)
                                         ggCCACATGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGc                                                                                           <  4:358477/149‑1 (MQ=255)
                                           ccACATGCCGCTGCCGCCCTATATCGAACGCCCGGACGAGGGCGCAGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCAGTCGCGGTTGACGCGCCTACAGAAGGACTGAACTTCGACGAGGAGCTGCTGGAAAatatagac                                                                                          >  8:229876/1‑144 (MQ=255)
                                             acaTGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGcc                                                                                                                                             <  4:164826/95‑1 (MQ=255)
                                             acaTGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGcc                                                                                                                                             >  3:164826/1‑95 (MQ=255)
                                                  ccgctgccgcCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAgg                                                                                   <  5:68907/148‑1 (MQ=255)
                                                    gctgccgcCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCg                                                                                >  7:111520/1‑149 (MQ=255)
                                                          gccCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAgc                                                                          <  6:107617/149‑1 (MQ=255)
                                                             cTATATCGACCGCCCGGACGAGGGCGCCGACCTTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATc                                                                                             <  5:271643/127‑1 (MQ=255)
                                                             cTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATc                                                                                             >  6:271643/1‑127 (MQ=255)
                                                             cTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAgcgcg                                                                       >  3:89111/1‑149 (MQ=255)
                                                              tataTCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTcg                                                                                                                    >  7:292621/1‑103 (MQ=255)
                                                              tataTCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTcg                                                                                                                    <  8:292621/103‑1 (MQ=255)
                                                               atatCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAgc                                                                          >  3:305486/1‑144 (MQ=255)
                                                               atatCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAgc                                                                          <  4:305486/144‑1 (MQ=255)
                                                               atatCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGcc                                                                     >  1:221320/1‑149 (MQ=255)
                                                                   cGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCg                                                                 >  7:3699/1‑149 (MQ=255)
                                                                   cGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCg                                                                 >  5:78347/1‑149 (MQ=255)
                                                                    gACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCg                                                                 <  6:78347/148‑1 (MQ=255)
                                                                     aCCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGt                                                                >  5:18363/1‑148 (MQ=255)
                                                                          ccGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCt                                                          >  3:204288/1‑149 (MQ=255)
                                                                                               gAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCtt                                     >  5:111774/1‑149 (MQ=255)
                                                                                                  cgcTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCa                                  <  4:302791/149‑1 (MQ=255)
                                                                                                  cgcTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCa                                  <  2:262637/149‑1 (MQ=255)
                                                                                                    cTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTc                                    <  4:113624/145‑1 (MQ=255)
                                                                                                    cTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTc                                    >  3:113624/1‑145 (MQ=255)
                                                                                                     tACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGcc                               >  2:129732/1‑149 (MQ=255)
                                                                                                       ccAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCgg                             <  8:111520/149‑1 (MQ=255)
                                                                                                          gACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGc                          <  6:111774/149‑1 (MQ=255)
                                                                                                               tgtACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGt                      <  6:18363/148‑1 (MQ=255)
                                                                                                               tgtACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCACAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTc                     >  5:363101/1‑149 (MQ=255)
                                                                                                                gtACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCg                              >  3:183927/1‑139 (MQ=255)
                                                                                                                gtACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCg                              <  4:183927/139‑1 (MQ=255)
                                                                                                                 tACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGc                                                        <  8:103123/112‑1 (MQ=255)
                                                                                                                 tACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGc                                                        >  7:103123/1‑112 (MQ=255)
                                                                                                                  aCGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGCCAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTCGa                   >  7:198734/1‑148 (MQ=255)
                                                                                                                  aCGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTCGa                   >  7:204028/1‑148 (MQ=255)
                                                                                                                   cGCCCAGCGCGCCGTCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTCGACa                 <  4:204288/149‑1 (MQ=255)
                                                                                                                    gCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTCGAc                  >  2:327233/1‑147 (MQ=255)
                                                                                                                      ccAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTCGAc                  >  1:199753/1‑145 (MQ=255)
                                                                                                                      ccAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTCGAc                  >  3:125297/1‑145 (MQ=255)
                                                                                                                           gcgcCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTCGACAAGATCGaa         <  7:124952/149‑1 (MQ=255)
                                                                                                                           gcgcCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTCGACAAGATCGaa         <  2:221320/149‑1 (MQ=255)
                                                                                                                            cgcCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTCGACAAGATCg           >  8:300656/1‑146 (MQ=255)
                                                                                                                            cgcCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTCGACAAGATCg           >  8:305683/1‑146 (MQ=255)
                                                                                                                             gcCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTCGACAAGATCGAAGac      >  1:71645/1‑149 (MQ=255)
                                                                                                                              taggCGCGGCAACCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTCGACAAGATCGAAg        >  7:38645/3‑147 (MQ=255)
                                                                                                                              taggCGCGGCAACCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTCGACAAGATCGAAg        <  8:38645/145‑1 (MQ=255)
                                                                                                                               cGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTCGACAAGATCGAAGacc     >  3:130058/1‑149 (MQ=255)
                                                                                                                                     ggTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTCGACAAGATCGAAGaccacc  <  2:196263/146‑1 (MQ=255)
                                                                                                                                     ggTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTCGACAAGATCGAAGaccacc  >  1:196263/1‑146 (MQ=255)
                                                                                                                                         |                                                                                                                                             
TTCGCTTCACCGAAGAGGTGCTGCCGTTGCTCGACCGTGTCGGCCACATGCCGCTGCCGCCCTATATCGACCGCCCGGACGAGGGCGCCGACCGTGAGCGCTACCAGACTGTGTACGCCCAGCGCGCCGGCGCGGTTGCCGCGCCTACCGCAGGCCTGCACTTCGACGAGGCGCTGCTGGAAAAGATCGCCGACAAGGGCGTCGAGCGCGCCTTCGTCACCCTGCATGTGGGCGCCGGCACCTTCCAGCCGGTGCGGGTCGACAAGATCGAAGACCACC  >  NC_002947/969467‑969745

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: