Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,820,019 T→C 20.9% E664G (GAG→GGG)  alaS ← alanyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,820,0190TC20.9% 122.4 / 2.6 65E664G (GAG→GGG) alaSalanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (33/15);  new base C (11/4);  total (46/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTGAGCGCACTTTATCAGCCAGATTATTGCTATCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACC  >  NC_000913/2819805‑2820246
                                                                                                                                                                                                                       |                                                                                                                                                                                                                                   
tACTGAGCGCACTTTATCAGCCAGATTATTGCTATCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAAACCGGAAGAGaaaaag                                                                                                                                                                                                                       >  2:44882/1‑226 (MQ=255)
     aGCGCACTTTATCAGCCAGATTATTGCTATCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGCACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACCCGGATGAGAAATCCCCCa                                                                                                                                                                                                                   >  1:40318/1‑228 (MQ=255)
     aGCGCACTTTATCAGCCAGATTATTGCTATCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGCTTCAGAGCTGCTGCGAAACAGCCCACTATCACCCGTCCGGCTGGCGTGAGTCACCCCCCCCACCCCGGTAGAGAAATCccccca                                                                                                                                                                                                                  >  2:280889/1‑227 (MQ=255)
     aGCGCACTTTATCAGCCAGATTATTGCTATCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGCTGATGCGGAACAGACCAAACTCACCAGTGCGGCTGGCGTGAGTACAGCCACCAAACACGGAAGAGAAATAGCCCaa                                                                                                                                                                                                                  >  2:506798/1‑228 (MQ=255)
     aGCGCACTTTATCAGCCAGATTATTGCTATCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCGACACAACACGGTAGAGAAATAGCCCAt                                                                                                                                                                                                                  >  1:660177/1‑229 (MQ=255)
     aGCGCACTTTATCAGCCAGATTATTGCTATCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACCCAACCCGGTAGAGAAATCGCCCa                                                                                                                                                                                                                   >  1:156950/1‑228 (MQ=255)
     aGCGCACTTTATCAGCCAGATTATTGCTATCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACACACATCACCAGCACGCCTGGCCGGAGAACCGACACACAACCCCGCAGAAAAAGCAcccca                                                                                                                                                                                                                  >  2:90774/1‑227 (MQ=255)
          aCTTTATCAGCCAGATTATTGCTATCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGATGTTGCGGAACAGCCCAATATCACCCGTGCGGCCGGCGGGAGGACCGCAACCCAACCCGGTAGAGaaacacccccggccca                                                                                                                                                                                                             >  2:451413/1‑216 (MQ=255)
          aCTTTATCAGCCAGATTATTGCTATCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGCCCAATATCACCAGTGCGGCTGGCGTGCGTACCGCCACACCACTCGGTAGAGAAATCGCCCATGCTc                                                                                                                                                                                                              >  1:18612/1‑228 (MQ=255)
          aCTTTATCAGCCAGATTATTGCTATCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCa                                                                                                                                                                                                             >  1:119192/1‑229 (MQ=255)
               aTCAGCCAGATTATTGCTATCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGGGAAATCGCCCATGCTCAGCAc                                                                                                                                                                                                         >  1:536495/1‑228 (MQ=255)
                        attattGCTATCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGCACGCCTGCAGCAGTACCCGCTTCAGAGCTGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGCGTACCCCCACACATCTCGGTAGAGAAATCGCCCATGCTCCGCACCCGAACgcgc                                                                                                                                                                                               >  1:717175/1‑229 (MQ=255)
                        attattGCTATCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGCTTACGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGGGAGGCACTCCACACAACCCGGTAAAAACACTCCCCATTCTCAACAcgccgaacc                                                                                                                                                                                                 >  2:106846/1‑221 (MQ=255)
                        attattGCTATCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTAACCGATTCCGAGATGATGCGGACCAGACCAATATCACCAGTGCGGCTGGCGGGGGGACCCCCACACAACCCGTTAGAAAAATCGCCCCTGCCCacacccccacaccc                                                                                                                                                                                                >  2:170255/1‑215 (MQ=255)
                                 aTCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCTGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGGGAGTACCGCCACACACCTCGGGAAGGAAATCGCCCATGCTCAGCACGCGTaaacaccccaaaaaa                                                                                                                                                                                      >  2:122653/1‑215 (MQ=255)
                                  tCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACCCGGTCGAGAAATCGCCCATGCCCAGCACGCCGACGCGCCCATCATAc                                                                                                                                                                                      >  2:519338/1‑228 (MQ=255)
                                  tCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGGACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACCCGGGAGAGAACACGCCCATGCTCAGAACGCGTACGCGCTACCCATACt                                                                                                                                                                                     >  1:217721/1‑229 (MQ=255)
                                  cccccTTTCAGGAGGTGGGCGCCTTCCCTAACGCGAGCGGTGTTGGCATGAACGGGGGCGATAGCCCCTTCCCCGGTTACCGTCTCGATACGCCGAACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGCCCAATACCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACt                                                                                                                                                                                     <  2:210726/226‑1 (MQ=255)
                                       tttCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACCCAACCCGGTAGAGAAATCGCCCATCCTCACACCGCGTCCGCGCTCACCCAACTtccc                                                                                                                                                                                 >  2:308504/1‑226 (MQ=255)
                                         tCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGCTGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGGGAAATCGCCCACGCCCAGCACGCGAACGCCCccatcacaatttcccc                                                                                                                                                                               >  2:216439/1‑218 (MQ=255)
                                                aTGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCCCAGCACCCGTACGCGCTCCCCAAACTTCTCGCCGAACAg                                                                                                                                                                        >  2:737766/1‑228 (MQ=255)
                                                           cGCTTAAGCGAGCGCTGTCTGCATGAACGGGGGCGATAGCCCCTTCCCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCa                                                                                                                                                              <  1:519338/227‑1 (MQ=255)
                                                                            tCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTATAGAAATCGCCCATGCTCAGCACCCGTACGCGCTCACCATACTTCTCCCCGAACAGCGCCATCCCACCCTTCCCTTTCCCCGCtt                                                                                                                                           >  2:426374/1‑229 (MQ=255)
                                                                            tCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATGCCACCTTTGCCTTTAGCCGCtt                                                                                                                                           >  1:180712/1‑229 (MQ=255)
                                                                                                     cTTCCCCGGTTACCGCTCCGATACGACGAACGCCTGCAGCAGTACCGGATTCAGAGATGATGCGGAACAGACCAATATCCCCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGa                                                                                                                   <  1:170255/228‑1 (MQ=255)
                                                                                                             gTTCCCGTTTCGATACAACGAACGCCGGCAGCAGTACCCGATTCAGGTATGAGGCGGACTAGACCAATATCCCCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAAt                                                                                                          <  1:737767/229‑1 (MQ=255)
                                                                                                                  cGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAAAAGCGCCATCGCCCCTTTCCCTTTCGCCGCTTCGAGATCCATGATGTTGGGTTCGATCGGCAAAtttctc                                                                                                     >  1:150929/1‑225 (MQ=255)
                                                                                                                         aaacgccgAACGCCGGCGGCGGGCCCCTGGTCAGGGATGCGGCGGAACAGACCAATTCCCCCAGGGCGGCTGGCGTGAGTGCCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCt                                                                                              <  2:119192/227‑1 (MQ=255)
                                                                                                                         aTACGACGAACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACAAATATCCCCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCt                                                                                              <  1:451413/229‑1 (MQ=255)
                                                                                                                               cgaACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCCCTTCGAGATCCATGATGTTGGTTTCGATCGGCCAATCGCGCCGTATCtgtgc                                                                                          >  1:732131/1‑226 (MQ=255)
                                                                                                                                  aCGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCCTCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCCTGATGTTGGTTTCGATCGGCCAATTGCGCCGAATCTGTGTGtttc                                                                                     >  2:221338/1‑227 (MQ=255)
                                                                                                                                    tcgtcgagcagTCCCGGATGCAGAGATGAGGCGGAACACACACGATCAGCACGCGCGCGTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTGCCCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCAcc                                                                                   <  2:18612/223‑1 (MQ=255)
                                                                                                                                    gcctggaagctgtaccGATTCAGGGAGGATGGGGAACAGACCAATATCCCCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCAcc                                                                                   <  1:506798/215‑1 (MQ=255)
                                                                                                                                    ccctgcacggttccccGATCCAGGGATGAGCCGGAACAGACCAATATCGCCAGTGCGCTGGGCGTGAGTTCCGCCACCCAACCCGGTAGGGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATCCTGCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCAcc                                                                                   <  2:115717/216‑1 (MQ=255)
                                                                                                                                    cccggcagcaGGCCCCGATTCAGAGAGGAGGCGGAACAGACCAATACCACCAGTGCGGCGGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCAcc                                                                                   <  1:280889/225‑1 (MQ=255)
                                                                                                                                         cagcagTACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCGCACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGCTGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCAcacggtt                                                                              >  2:318928/1‑223 (MQ=255)
                                                                                                                                                    gATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCCCCAGGCCTTCGCCCGCAc                                                                   >  1:672696/1‑229 (MQ=255)
                                                                                                                                                    gATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCACTCATCATACTTCTCGCCGAACAGCGCCATCACACCTTTCCCCTTCGCCGCCTCGGACTCCATGATTTTGGTTTCCACCGGCCAATTGCGCCGAACCTGTGTGTTCAACAGGACTTCGAACGcca                                                                   >  2:229833/1‑227 (MQ=255)
                                                                                                                                                         agagTTGGTGCGGAACAGACCAATACCGCCAGTGCGGCTGGGGTGAGTGCCGCCCCACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCACGAAt                                                               <  2:156950/228‑1 (MQ=255)
                                                                                                                                                         agagATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCgag                                                                                                                                       >  1:43552/1‑156 (MQ=255)
                                                                                                                                                           agaTGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGGGTTCACCAGGTCTTCGACCGCACGCAtctc                                                            >  1:702357/1‑229 (MQ=255)
                                                                                                                                                               gatgCGAAACAGCCCAAGACCCCCAGGGGGGCTGGCGTGAGTACCGCCACCCAACTCGGTAGAGAAATCGCCCATGCCCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCACGAATCtcttct                                                        <  1:712583/229‑1 (MQ=255)
                                                                                                                                                                atgCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCCCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGGTGTTGGTTTCGATCGGCAAATTGCGCCGAACCAGTGGTTTAACCAGGTTTTGTACCGCACGCATCtcttc                                                         >  2:304240/1‑227 (MQ=255)
                                                                                                                                                                     gAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTc                                                  >  1:367211/1‑229 (MQ=255)
                                                                                                                                                                           aCCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCCCGAATCTCTTCTGGTTTCATCGCt                                            >  1:504700/1‑229 (MQ=255)
                                                                                                                                                                           aCCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTGGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCCCGAACTCCTTTGGGTTTAATCGCt                                            >  2:431424/1‑229 (MQ=255)
                                                                                                                                                                                tatCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCTTGATGTTGGTTTCCATCGGCAAATTGTGACGAATCTATGTGTTCAcc                                                                                   >  2:494005/1‑185 (MQ=255)
                                                                                                                                                                                tatCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGCTCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCCCTTCGtt                                       >  1:473097/1‑229 (MQ=255)
                                                                                                                                                                                tatCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAACTTGCGACGCATCTGTGTGTTCACCAGGTCTTCGACCGCACGCATCTCTACTGGTTTCTTCGCTTCCtt                                       >  1:700094/1‑229 (MQ=255)
                                                                                                                                                                                tatCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCCCGCATCTCTTCTGGTTTCATCGCTTCGt                                        >  1:725871/1‑228 (MQ=255)
                                                                                                                                                                                tatCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGt                                        >  1:446889/1‑228 (MQ=255)
                                                                                                                                                                                     ccAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGCGATCCATGATGTTGGTTTCGATCGGCCAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTtttgtg                                  >  1:173866/1‑229 (MQ=255)
                                                                                                                                                                                     ccAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGCCGAATCTTGGTGTTCACAAGGCTTTTGACCGCACGAATTCCTTCTGGTTTCATCGTTTTtgtgtgtg                                  >  2:250651/1‑229 (MQ=255)
                                                                                                                                                                                         gTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTAGAGATCCATGATGTTGGTTTCGATCGGCAAATTCCGACGAATCTTGGTGTTCACCAGGGCTTCGACCGCACGAATCCCTTCTGGTTTCAACGCttgcgtgggggcaa                               >  2:392951/1‑217 (MQ=255)
                                                                                                                                                                                            cGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAg                             >  1:399129/1‑228 (MQ=255)
                                                                                                                                                                                            cGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGt                            >  1:31509/1‑229 (MQ=255)
                                                                                                                                                                                              gacggtctgggtcccGCCACACAAATCGGTAGAAAACACGCCCATGCCTAGCACGCGTTCGCGATCATCCGTCTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCg                          <  2:269512/216‑1 (MQ=255)
                                                                                                                                                                                              gCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCAAGCGTACGCGCTCATCATACTTATCGCCGAAAAGCGCCATCGCAACTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTATTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTACGTTGTGTgagactga                           >  1:269512/1‑224 (MQ=255)
                                                                                                                                                                                                      gAGTCCCCCCACACAACTCGGGAGGGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCgca                   <  1:440749/228‑1 (MQ=255)
                                                                                                                                                                                                           cccccACACAACCCGGTAGAGAAATCCCCCATGCTCAGCACGCGTACGCGCTCATCATTCTTCTTGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCAcc              <  1:44882/225‑1 (MQ=255)
                                                                                                                                                                                                           ccGCCACACAACCCGGTAGAGAAACCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCAcc              <  1:216439/228‑1 (MQ=255)
                                                                                                                                                                                                              ccACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTg           <  1:250651/228‑1 (MQ=255)
                                                                                                                                                                                                                  ccaACCCGGTAGATAAACGCCCCATGCCCACGACGCTTACCCGCTCATCATTCGTCTCGCCAAACAGCGCCATCGCTCCTTTCGCTTTGGCCGCTTCGAGTCCCGTGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGtt      <  2:180712/228‑1 (MQ=255)
                                                                                                                                                                                                                    aaCTCGGTAGACAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTa     <  1:221338/228‑1 (MQ=255)
                                                                                                                                                                                                                      cTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCCCGAATCTCTTTTGGTTTCATCGCTTCGTTGTGTGTGAAGTCGAAGCGCAGAACCTTGTCGTTAAcc  >  1:655919/1‑229 (MQ=255)
                                                                                                                                                                                                                       |                                                                                                                                                                                                                                   
TACTGAGCGCACTTTATCAGCCAGATTATTGCTATCGCCTTTCAGCAGATGCGCGACTTCGCTTAAGCGATCGCTGTCTGCATGAACGGTGGCGATAGCACCTTCTCCGGTTACCGCTTCGATACGACGAACGCCTGCAGCAGTACCCGATTCAGAGATGATGCGGAACAGACCAATATCACCAGTGCGGCTGGCGTGAGTACCGCCACACAACTCGGTAGAGAAATCGCCCATGCTCAGCACGCGTACGCGCTCATCATACTTCTCGCCGAACAGCGCCATCGCACCTTTCGCTTTCGCCGCTTCGAGATCCATGATGTTGGTTTCGATCGGCAAATTGCGACGAATCTGTGTGTTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACC  >  NC_000913/2819805‑2820246

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: