Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,082,030 C→A 6.9% E746* (GAA→TAA) ‡ alaS ← alanine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,082,0300CA 6.9% 185.4 / 3.1 58E746* (GAA→TAA) ‡alaSalanine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (24/30);  new base A (2/2);  total (26/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.50e-01

TTGTTCTTCAGCTGATCGACCAGGGCCAGCAGGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTGCTCTTCGGCAGCGTTCAGGTAAGCCAGCGCAGCGGCACCGGTTACG  >  NC_002947/5081882‑5082178
                                                                                                                                                       |                                                                                                                                                    
ttGTTCTTCAGCTGATCGACCAGGGCCAGCAGGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTc                                                                                                                                                      <  4:80811/149‑1 (MQ=255)
  gTTCTTCAGCTGATCGACCAGGGCCAGCAGGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTc                                                                                                                                                      <  5:118108/147‑1 (MQ=255)
    tcttcAGCTGATCGACCAGGGCCAGCAGGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAgc                                                                                                                                                  >  8:82643/1‑149 (MQ=255)
      ttcAGCTGATCGACCAGGGCCAGCAGGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAgctg                                                                                                                                                >  2:274803/1‑149 (MQ=255)
        cAGCTGATCGACCAGGGCCAGCAGGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAgctgc                                                                                                                                               <  1:277542/148‑1 (MQ=255)
             gATCGACCAGGGCTAGCAGGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTtctc                                                                                                                                          >  7:286271/1‑148 (MQ=255)
                 gACCAGGGCCAGCAGGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAgctgc                                                                                                                                               >  3:222056/1‑139 (MQ=255)
                 gACCAGGGCCAGCAGGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAgctgc                                                                                                                                               <  4:222056/139‑1 (MQ=255)
                  aCCAGGGCCAGCAGGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGGAACAgctgcag                                                                                                                                             >  1:248848/1‑138 (MQ=255)
                  aCCAGGGCCAGCAGGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGGAACAgctgcag                                                                                                                                             <  2:248848/140‑3 (MQ=255)
                   ccAGGGCCAGCAGGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTgg                                                                                                                                   >  5:99739/1‑149 (MQ=255)
                      gggCCAGCAGGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCgg                                                                                                                                >  8:180249/1‑149 (MQ=255)
                       ggCCAGCAGGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGt                                                                                                                               >  5:80497/1‑149 (MQ=255)
                       ggCAAGCAGGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGt                                                                                                                               >  7:304791/1‑149 (MQ=255)
                         ccAGCAGGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCAGCAGCTGGAACAgctgcagg                                                                                                                                            <  6:40932/134‑4 (MQ=255)
                         ccAGCAGGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCAGCAGCTGGAACAgctgcagg                                                                                                                                            >  5:40932/1‑131 (MQ=255)
                            gcagGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCGTCGCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTgc                                                                                                                            <  1:274803/149‑1 (MQ=255)
                             cagGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCt                                                                                                                         <  6:80497/149‑1 (MQ=255)
                                ggCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCT‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCaa                                                                                                                      <  7:180249/149‑1 (MQ=255)
                                ggCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCaa                                                                                                                      <  8:170433/149‑1 (MQ=255)
                                  ccTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGc                                                                                                                    <  8:286271/149‑1 (MQ=255)
                                   cTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGca                                                                                                                   >  1:149405/1‑149 (MQ=255)
                                        cATCCTGTCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGcc                                                                                                              >  2:145809/1‑149 (MQ=255)
                                        cATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGcc                                                                                                              >  2:187465/1‑149 (MQ=255)
                                               cccATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGc                                                                                                        <  4:184205/148‑1 (MQ=255)
                                               cccATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCt                                                                                                       <  6:33516/149‑1 (MQ=255)
                                                ccATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGTCGGCTGTGCTTGAGCGCAACCGCCAGCtg                                                                                                      <  2:142496/149‑2 (MQ=38)
                                                   tCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTc                                                                                                   <  8:220773/149‑1 (MQ=255)
                                                     cAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGa                                                                                                 <  8:304791/149‑1 (MQ=255)
                                                      aGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGa                                                                                                 >  6:72027/1‑148 (MQ=255)
                                                        gcgggcggCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCGTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCa                                                                                              <  7:82643/149‑1 (MQ=255)
                                                                    gACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGt                                                                                  <  2:149405/149‑1 (MQ=255)
                                                                       cTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTAc                                                                               >  3:72105/1‑149 (MQ=255)
                                                                          ggCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTAccc                                                                             <  5:72027/148‑1 (MQ=255)
                                                                               ccTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGa                                                                         <  1:42463/147‑1 (MQ=255)
                                                                               ccTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGa                                                                         >  2:42463/1‑147 (MQ=255)
                                                                                 ttGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAg                                                                     >  4:88106/1‑149 (MQ=255)
                                                                                      cTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGtt                                                                                                                              <  8:193382/87‑1 (MQ=255)
                                                                                      cTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGtt                                                                                                                              >  7:193382/1‑87 (MQ=255)
                                                                                      cTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTtgtg                                                                <  5:350203/149‑1 (MQ=255)
                                                                                       tCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGc                                                               >  7:93019/1‑149 (MQ=255)
                                                                                       tCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGc                                                               >  3:61215/1‑149 (MQ=255)
                                                                                             gCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTc                                                         <  1:288315/149‑1 (MQ=255)
                                                                                                      gagagaTCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCt                                                 >  1:218717/1‑149 (MQ=255)
                                                                                                      gagagaTCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTg                                                >  2:172694/1‑149 (MQ=255)
                                                                                                      gagagaTCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTg                                                >  8:300162/1‑149 (MQ=255)
                                                                                                              gTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTGCTCTTCgg                                        <  6:99739/149‑1 (MQ=255)
                                                                                                                 cccG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTGCTCTTCGGCa                                      <  7:300162/148‑1 (MQ=255)
                                                                                                                       gcggcgCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTGCTCTTCGGCAGCGtt                                 >  7:182384/1‑149 (MQ=255)
                                                                                                                          gcgcTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTGCTCTTCGGCAGCGTTCAg                              <  1:145809/149‑1 (MQ=255)
                                                                                                                          gcgcTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTGCTCTTCGGCAGCGTTCAg                              <  4:72105/149‑1 (MQ=255)
                                                                                                                          gcgcTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTGCTCTTCGGCAGCGTTCAg                              <  4:61215/149‑1 (MQ=255)
                                                                                                                           cgcTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTGCTCTTCGGCAGCGTTCAgg                             <  1:172694/149‑1 (MQ=255)
                                                                                                                            gcTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTGCTCTTCGGCAGCGTTCAgg                             >  6:235862/1‑148 (MQ=255)
                                                                                                                                  ggcCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTGCTCTTCGGCAGCGTTCAGGTAAGCca                      <  1:135933/149‑1 (MQ=255)
                                                                                                                                     cTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTGCTCTTCGGCAGCGTTCAg                              <  2:290230/138‑1 (MQ=255)
                                                                                                                                     cTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTGCTCTTCGGCAGCGTTCAg                              >  1:290230/1‑138 (MQ=255)
                                                                                                                                               cAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTGCTCTTCGGCAGCGTTCAGGTAAGCCAGCGCAGCGGCAc          <  5:235862/149‑1 (MQ=255)
                                                                                                                                                 gCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTGCTCTTCGGCAGCGTTCAGGTAAGCCAGCGCAGCGGCACCg        <  8:182384/149‑1 (MQ=255)
                                                                                                                                                 gCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTGCTCTTCGGCAGCGTTCAGGTAAGCCAGCGCAGCGGCACCg        <  2:218717/149‑1 (MQ=255)
                                                                                                                                                       ccAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTGCTCTTCGGCAGCGTTCAGGTAAGCCAGCGCAGCGGCACCGGTTAc   <  3:88106/148‑1 (MQ=255)
                                                                                                                                                       ccAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTGCTCTTCGGCAGCGTTCAGGTAAGCCAGCGCAGCGGCACCGGTTACg  >  5:231625/1‑149 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
TTGTTCTTCAGCTGATCGACCAGGGCCAGCAGGGCCTTGCCATCCTGCCCATCCAGGCGGGCGGCGAGGACCTTGGCACCCTTGACCTCAACGGCCGCGTTGGAGAGATCGTCCCCG‑‑GCGGCGCTGGCGGCCTTGGCCTGCAGCTGTTCCAGCTGCTTCTCCAGCTGGCGGTTGCGCTCAAGCACAGCCGACAGCTTGTCGATCAGGTTGTCGCGGTTACCCTTGACCAGTTGTGCGGCTTCCTTGACCTGCTCTTCGGCAGCGTTCAGGTAAGCCAGCGCAGCGGCACCGGTTACG  >  NC_002947/5081882‑5082178

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: