Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,516,192 C→G 6.1% Y19* (TAC→TAG ate → arginyl‑tRNA‑‑protein transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,516,1920CG 6.1% 207.6 / 4.2 66Y19* (TAC→TAGatearginyl‑tRNA‑‑protein transferase
Reads supporting (aligned to +/- strand):  ref base C (29/32);  new base G (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

CTTGACCAGCCCAACAGCGCTTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGCCCACACTGCCAGAACTGCAACGCCTGCGTGC  >  NC_002947/4516044‑4516340
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                                           >  4:207565/1‑149 (MQ=255)
                           gTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTg                                                                                                                           <  2:383739/149‑1 (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                                                                 <  2:416211/109‑1 (MQ=255)
                                                                                                                             aCTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATgg                                                                 >  1:416211/1‑109 (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                                <  3:234085/126‑1 (MQ=255)
                                                                                                                                          ccggCAGGTAGCTGCAGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGcc                                >  4:234085/4‑129 (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)
                                                                                                                                                    ccTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGCCCACACTGCCAGAACTGCAACGCCTGCGTGc  <  5:27864/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTTGACCAGCCCAACAGCGCTTCATGGGTTGCCTAGGCGAGTTCCCGTAGCTGGCTTACACTTAAAAAAAAGTCTCACCGAAGGGGTTGATCATGACAGAGTTGGCGCGGTTGAAGTTCTATGCCACTCAGCCCCACTCCTGCAGCTACCTGCCGGAGGAGCAGGCGACCACGCTGTTCCTCGACCCCAGCCAGCCGATGGACGTGCACGTGTACGCCGATTTGTCGGAAATGGGCTTTCGCCGCAGTGGCGACCACCTGTACCGCCCACACTGCCAGAACTGCAACGCCTGCGTGC  >  NC_002947/4516044‑4516340

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: