Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,261,924 G→A G86S (GGC→AGC)  dusA → tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,261,9240GA98.0% 162.6 / NA 49G86S (GGC→AGC) dusAtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/24);  minor base C (1/0);  total (25/24)
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

CCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGGGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGT  >  NC_000913/4261630‑4262215
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                    
ccccATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGt                                                                                                                                                                                                                                                                                                 >  1:60403/1‑300 (MQ=255)
                      aatCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGgcgc                                                                                                                                                                                                                                                                           <  2:357945/300‑1 (MQ=255)
                                tttGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCaaa                                                                                                                                                                                                                                                                 <  2:348696/300‑1 (MQ=255)
                                 ttGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGTAGTGTTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACATGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCaaa                                                                                                                                                                                                                                                                 <  2:1973/299‑1 (MQ=255)
                                  tGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCTCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGTCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGc                                                                                                                                                                                                                                                               <  2:292780/300‑1 (MQ=255)
                                  tGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGc                                                                                                                                                                                                                                                               <  2:327753/300‑1 (MQ=255)
                                  tGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGc                                                                                                                                                                                                                                                               <  1:377434/300‑1 (MQ=255)
                                       aTGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTNGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCa                                                                                                                                                                                                                                                          <  1:129843/300‑1 (MQ=255)
                                              gTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAgcgc                                                                                                                                                                                                                                                   >  1:548853/1‑300 (MQ=255)
                                                                          aaccaGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGt                                                                                                                                                                                                                       >  1:255062/1‑300 (MQ=255)
                                                                          aaccaGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCAGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAAACTGaacgt                                                                                                                                                                                                                       >  2:365814/1‑297 (MQ=255)
                                                                                                        tCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAAc                                                                                                                                                                                           >  1:12260/1‑298 (MQ=255)
                                                                                                                 tGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTgg                                                                                                                                                                                <  2:562905/300‑1 (MQ=255)
                                                                                                                      ggTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGTTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCg                                                                                                                                                                            <  1:71577/299‑1 (MQ=255)
                                                                                                                       gTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTCTCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCAGGTGCAGAAAGGCATGTTTGGTgtgt                                                                                                                                                                           >  2:143005/1‑296 (MQ=255)
                                                                                                                       gTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGCGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGACTGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCgt                                                                                                                                                                           >  2:494998/1‑299 (MQ=255)
                                                                                                                       gTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCgt                                                                                                                                                                           >  2:548658/1‑299 (MQ=255)
                                                                                                                            gCGTTGCACCAATGCTCGACTGGACGGACAGACATTTCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGa                                                                                                                                                                     >  1:488129/1‑300 (MQ=255)
                                                                                                                                   aCCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTaa                                                                                                                                                              <  1:348358/300‑1 (MQ=255)
                                                                                                                                      aaTGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGttctgtctcttatacacatctcgcag                                                                                                                                                           <  1:384576/300‑25 (MQ=255)
                                                                                                                                        tGCTCGACGGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGTTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATgcgc                                                                                                                                                         <  2:119055/300‑1 (MQ=255)
                                                                                                                                                    cggacAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGcc                                                                                                                                              >  1:528980/1‑299 (MQ=255)
                                                                                                                                                                        ttCTTGCGTCTGCTTTCCCGCAATACTTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTNATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATgcgc                                                                                                                         <  1:114585/300‑1 (MQ=255)
                                                                                                                                                                           ttGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGAt                                                                                                                      <  2:328890/300‑1 (MQ=255)
                                                                                                                                                                                  tGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGgtgt                                                                                                               <  1:525956/300‑1 (MQ=255)
                                                                                                                                                                                         ccGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTcc                                                                                                       >  2:120125/1‑300 (MQ=255)
                                                                                                                                                                                                   ttGCTGTATACCGAAANGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAANGCGATGCGCGATGTGGTGNCGATTCcggtgacggt                                                                                              >  1:521445/1‑299 (MQ=255)
                                                                                                                                                                                                     gttGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGCAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAGGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGAATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGaa                                                                                           >  2:250852/3‑300 (MQ=255)
                                                                                                                                                                                                     gCTGTATACCGAAATGGTGACCACAGGGGCGATTCTTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAGGAACATCCGGTAGCGTTGAAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGTAGTGTGCAAAGTTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGATGCACGTATGAACGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGaa                                                                                           >  2:251566/1‑300 (MQ=255)
                                                                                                                                                                                                     gCTGTACACCGAAATGGTGAACACAGGAGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAATGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGAGCTGGCGCAGTGTGCAAAGCTGGAAGAAGCGTGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTAGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGATGTCGACTCaggtgacgtcgaa                                                                                           >  2:251747/1‑295 (MQ=255)
                                                                                                                                                                                                                  aTGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGAGTATTGGc                                                                              >  2:156182/1‑300 (MQ=255)
                                                                                                                                                                                                                    ggTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCAt                                                                            >  2:596383/1‑300 (MQ=255)
                                                                                                                                                                                                                       gACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATcgac                                                                        <  1:218314/300‑1 (MQ=255)
                                                                                                                                                                                                                        aCCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATcgac                                                                        <  1:107790/300‑1 (MQ=255)
                                                                                                                                                                                                                                 gctattattCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGAc                                                               <  2:90427/297‑1 (MQ=255)
                                                                                                                                                                                                                                 gCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGAc                                                               <  1:600357/300‑1 (MQ=255)
                                                                                                                                                                                                                                         tCACGGTCAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCNGATTCcggtgacggtg                                                                                             <  1:510326/263‑1 (MQ=255)
                                                                                                                                                                                                                                              gTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTtc                                                  >  1:367274/1‑300 (MQ=255)
                                                                                                                                                                                                                                                   ggTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGc                                             <  2:518566/300‑1 (MQ=255)
                                                                                                                                                                                                                                                       aTTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCNATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGAtt                                         >  1:132041/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                 gTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATAc                               <  1:412738/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gTGAAGAAGAACATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCATAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGttt                          >  1:42934/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                        gaagaagaaCATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGAACGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTcc                        >  1:565538/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                         acgaagaaCATCCGGTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACACGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTATGCGATTTCATCAATACCGTTTcc                        >  2:298052/3‑299 (MQ=255)
                                                                                                                                                                                                                                                                                      ggTAGCGTTGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGANGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgtg          <  1:442971/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             ttGCAGCTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCACGTCTNACCGGGTGCAGAACGGCATGTTTGGTGCGTGTNTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCAGGCAAAGGCCAGTGTGGGAt    >  2:458153/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                             ttGCAACTGAGTGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGCCTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGAACAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGCGAt    >  2:133265/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                             ttGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTACAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAGTTTCTATGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga       >  2:132816/1‑296 (MQ=255)
                                                                                                                                                                                                                                                                                              tGCAACTGAGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTTAGATGt  <  1:532551/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                    
CCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGTGACCACAGGGGCGATTATTCACGGTAAAGGTGATTATCTGGCGTACAGTGAAGAAGAACATCCGGTAGCGTTGCAACTGGGCGGTAGCGATCCGGCGGCGCTGGCGCAGTGTGCAAAGCTGGCAGAAGCGCGCGGATATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGT  >  NC_000913/4261630‑4262215

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