Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,516,189 C→G 5.8% S18R (AGC→AGG)  ate → arginyl‑tRNA‑‑protein transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,516,1890CG 5.8% 226.1 / 4.1 69S18R (AGC→AGG) atearginyl‑tRNA‑‑protein transferase
Reads supporting (aligned to +/- strand):  ref base C (31/34);  new base G (2/2);  total (33/36)
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.38e-01

CACCTTGACCAGCCCAACAGCGCTTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGCCCACACTGCCAGAACTGCAACGCCTGC  >  NC_002947/4516041‑4516336
                                                                                                                                                    |                                                                                                                                                   
cACCTTGACCAGCCCAACAGCGCTTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGc                                                                                                                                                     <  2:240454/149‑1 (MQ=255)
  ccTTGACCAGCCCAACAGCGCTTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTa                                                                                                                                                   <  3:21634/149‑1 (MQ=255)
   cTTGACCAGCCCAACAGCGCTTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTAc                                                                                                                                                  >  7:51083/1‑149 (MQ=255)
   cTTGACCAGCCCAACAGCGCTTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTAc                                                                                                                                                  >  4:308111/1‑149 (MQ=255)
      gACCAGCCCAACAGCGCTTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTg                                                                                                                                               >  4:94093/1‑149 (MQ=255)
        ccAGCCCAACAGCGCTTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGcc                                                                                                                                             <  2:404332/149‑1 (MQ=255)
         cAGCCCAACAGCGCTTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCg                                                                                                                                            >  6:256632/1‑149 (MQ=255)
            cccAACCGCGCGTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTAAAAAAAAAGTCTCACCGAAGGTGTTGATCATGACAGAGTTGGCGCGGTTTAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCggag                                                                                                                                         <  8:343343/149‑1 (MQ=255)
            cccAACAGCGCTTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCggag                                                                                                                                         <  2:418370/149‑1 (MQ=255)
            cccAACAGCGCATCATTGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGTGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCggag                                                                                                                                         <  5:221146/149‑1 (MQ=255)
              cAACAGCGCTTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCggagga                                                                                                                                       <  2:101041/149‑1 (MQ=255)
               aaCAGCGCTTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAATAAAAAGTCTCACCGAAGTGGTTGATAATGACAGAGTTGGCGCGGTTGAAGTTCTATGGCACTCAGCCCCACTCCTGCAGCTACCTGCCGgaggag                                                                                                                                      >  2:41414/1‑149 (MQ=255)
                aCAGCGCTTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCGGCAGGTAGCTGCAggagtggg                                                                                                                                     <  3:325814/149‑5 (MQ=255)
                  aGCGCTTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAg                                                                                                                                   >  8:232541/1‑149 (MQ=255)
                    cgcTTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGc                                                                                                                                 <  8:51083/149‑1 (MQ=255)
                    cgcTTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGc                                                                                                                                 <  6:275215/149‑1 (MQ=255)
                         cATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGAcc                                                                                                                             <  7:183992/148‑1 (MQ=255)
                              gTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTg                                                                                                                       <  2:383739/149‑1 (MQ=255)
                              gTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTg                                                                                                                       >  4:207565/1‑149 (MQ=255)
                                  ccTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTcc                                                                                                                   <  7:232541/149‑1 (MQ=255)
                                  ccTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTc                                                                                                                    <  3:306951/148‑1 (MQ=255)
                                    tAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTc                                                                                                                 >  6:170407/1‑149 (MQ=255)
                                          aGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGAccc                                                                                                            >  4:216024/1‑148 (MQ=255)
                                                     cTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGAt                                                                                                <  5:359972/149‑1 (MQ=255)
                                                        gCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGa                                                                                             <  5:256632/149‑1 (MQ=255)
                                                          ttACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACg                                                                                           >  7:119048/1‑149 (MQ=255)
                                                              acTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCATGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCa                                                                                       >  1:194487/1‑149 (MQ=255)
                                                               cTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCAc                                                                                      >  7:237889/1‑149 (MQ=255)
                                                                ttAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCCGCCAGCCGATGGACGTGCACg                                                                                     >  4:248557/1‑149 (MQ=255)
                                                                ttAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACg                                                                                     >  5:411369/1‑149 (MQ=255)
                                                                  aaaaaaaaGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACgt                                                                                    <  6:379643/148‑1 (MQ=255)
                                                                   aaaaaaaGACTCCCCGAAGGTGTTGATCAGGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACgt                                                                                    <  4:359953/147‑1 (MQ=255)
                                                                         aGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCg                                                                            <  3:216024/149‑1 (MQ=255)
                                                                             tcACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGGCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCATCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGAttt                                                                        >  4:247757/1‑149 (MQ=255)
                                                                               aCCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGt                                                                      >  1:131654/1‑149 (MQ=255)
                                                                                ccGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTc                                                                     <  7:80314/149‑1 (MQ=255)
                                                                                ccGAAGGGGATGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTc                                                                     <  7:380210/149‑1 (MQ=255)
                                                                                 cGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCAGATTTGTCg                                                                    >  4:124453/1‑149 (MQ=255)
                                                                                 cGAAGCGGTTGATCATGACAGAGTTGGCGCGGGTGAAGTGCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCg                                                                    <  4:100392/149‑1 (MQ=255)
                                                                                         ttGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATggg                                                            <  8:237889/149‑1 (MQ=255)
                                                                                           taTCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGc                                                           <  5:355358/147‑1 (MQ=255)
                                                                                            aTCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCt                                                          <  8:119048/148‑1 (MQ=255)
                                                                                            aTCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCt                                                          <  6:411369/148‑1 (MQ=255)
                                                                                               aTGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTc                                                       <  3:308111/148‑1 (MQ=255)
                                                                                                             gcgGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGac                                         >  3:245165/1‑148 (MQ=255)
                                                                                                             gcgGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTCCCGGAGGAGCACGCGACCACGCTGATACTCGACCCCAGCCAGCCGATCGACGTGCACGTCTACGCCCATATGTCGTAAATGGGCTTTCGACGCAgatgcgcc                                         >  5:414532/1‑141 (MQ=255)
                                                                                                                 ttGAAGTTCTATGCCACTCAGCCCCACACCTGCAGCTACCTGCCTGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCACCCAGCCGATGGACGTGCACGTGAACGCCGATTTGTTCGAAATGGGCTTTCGCCGCAGTGGCGaccac                                      <  1:41414/147‑1 (MQ=255)
                                                                                                                     aGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTAc                                <  5:170407/149‑1 (MQ=255)
                                                                                                                       ttCTATGCCACTCAGTCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGATCACCTGTACCg                              >  4:130816/1‑149 (MQ=255)
                                                                                                                         cTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGcc                            >  4:209286/1‑149 (MQ=255)
                                                                                                                         cTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGc                             <  3:248557/148‑1 (MQ=255)
                                                                                                                           aTGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGAGGCGACCACCTGTACCGCCca                          >  4:226531/1‑149 (MQ=255)
                                                                                                                              ccACTCACCCCCACTCATGCAGCTAACTTCCGTCGGAGCAGGCTACCACGCTGTTCATCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCTCAGTGGCGACCACCTGTACCGCCCacac                       <  3:51749/149‑1 (MQ=255)
                                                                                                                                aCTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGCCCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGCCCACACTg                     >  5:37440/1‑149 (MQ=255)
                                                                                                                                aCTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATgg                                                             >  1:416211/1‑109 (MQ=255)
                                                                                                                                aCTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATgg                                                             <  2:416211/109‑1 (MQ=255)
                                                                                                                                    aGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGCCCACACTGcc                   <  3:226531/147‑1 (MQ=255)
                                                                                                                                     ctgctccTCCGGCAGGTAGCTGCAGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGCCCACACTGCCAg                 >  7:187311/7‑148 (MQ=255)
                                                                                                                                      ccccACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGCCCACACTGCCAGaa               >  1:410255/1‑149 (MQ=255)
                                                                                                                                      ccccACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGCCCACACTGCCAGaa               >  7:180616/1‑149 (MQ=255)
                                                                                                                                       cccACTCCTGCAGCTACCTGCGGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGCCCACACTGCCAGAAc              >  5:271470/1‑149 (MQ=255)
                                                                                                                                          aCTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGCCCACACTGCCAGAACTg            >  6:45602/1‑148 (MQ=255)
                                                                                                                                          aCTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGCCCACACTGCCAGAACTg            <  5:45602/148‑1 (MQ=255)
                                                                                                                                             ccggCAGGTAGCTGCAGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGcc                            >  4:234085/4‑129 (MQ=255)
                                                                                                                                             ccggCAGGTAGCTGCAGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGcc                            <  3:234085/126‑1 (MQ=255)
                                                                                                                                              cTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGCCCACACTGCCAGAACTGCAAc        <  2:131654/148‑1 (MQ=255)
                                                                                                                                              cTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGCCCACACTGCCAGAACTGCAACg       >  6:27864/1‑149 (MQ=255)
                                                                                                                                              cTGCAGCTACCTGCCGGAGGAGAAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTCTACCCCGAGTTGGCGGAAATGGGCGATCGCCGCAGTGGCGACCACCACTACCGCCCACACTGCCAGAACTGCAACg       >  7:66501/1‑149 (MQ=255)
                                                                                                                                               tGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGCCCACACTGCCAGAACTGCAACGc      >  4:177706/1‑149 (MQ=255)
                                                                                                                                                   gCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGCCCACACTGCCAGAACTGCAACGCCTGc  <  3:94093/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CACCTTGACCAGCCCAACAGCGCTTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGCCCACACTGCCAGAACTGCAACGCCTGC  >  NC_002947/4516041‑4516336

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: ‑