Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,761,620 T→G 9.2% K641Q (AAG→CAG)  yfjK ← radiation resistance protein; DEAD/H helicase‑like protein; CP4‑57 putative defective prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,761,6200TG 9.2% 1208.8 / 101.0 370K641Q (AAG→CAG) yfjKradiation resistance protein; DEAD/H helicase‑like protein; CP4‑57 putative defective prophage
Reads supporting (aligned to +/- strand):  ref base T (174/162);  new base G (14/20);  total (188/182)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-02

Only 100 of 381 total aligned reads displayed.

AAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTGCGACTGAACGTATGGGGAACCAGTC  >  NC_000913/2761348‑2761892
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                
aaaTCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAACACGCTTTa                                                                                                                                                                                                                                                                                >  2:293149/1‑275 (MQ=255)
            tcTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGCGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGCAAATTCCAGCAAGTCTATCAGNCGGTGTAACGTCTGAATGGGAAGGCCAATCTCTTCTAAGGAATTGACACATGGTGGTAGCTGGAAACTTTCGAACGCCAANTTGACAGGTGTGTAATGAACGTTGGNGGGGATAAACCGGTTTACAGCGTGGGGtt                                                                                                                                                                                                                                                                    >  3:163746/1‑275 (MQ=255)
             cTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtgg                                                                                                                                                                                                                                                                        >  4:417925/1‑270 (MQ=255)
                      atcgggcgcCTGATAAAAAGCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCGGAGTATCCCGAGGATACTGGGTCAAATCATCCACATCAGCCTCGTCGGAAACTTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGTGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtc                                                                                                                                                                                                                                                          <  4:310865/270‑1 (MQ=255)
                      atcgcgcgcCTGATAAACATCTGATCAACGTACCCACTCGATCTACTCCAAATGTCCTGAGTCTCCCGAAGATACTGGCTCAATTCCTCCACATCAGCCTCGTCAGAAACTTCCAGCCAGTCTCCCAGTCGGTGTAACGTCTGCATGGGAATGCCAATATCTTCTAAGGCCTTGACACCTTGTGGTAGATGGAAACTTTCGACCGCCAATTTGACATGTCTGTAATCCACTTTGCTGCTTCTACCACGCTTTCCCGCATGGTTTTTAACTACAtc                                                                                                                                                                                                                                                          >  4:84688/3‑275 (MQ=255)
                      aGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAgtg                                                                                                                                                                                                                                                                           >  1:211879/1‑256 (MQ=255)
                            ccccTAATAAACAACTGATCAACCGAACCAATCGACCTCCTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCANCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAg                                                                                                                                                                                                                                                <  3:209916/277‑1 (MQ=255)
                                              tCAACGTACCCAAACGGTTGACTCCAAAGGTCCTGAGGATCCCTAAGATACTGGCTCAACGCATCAACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAAGCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtagg                                                                                                                                                                                                                                  <  4:186433/275‑1 (MQ=255)
                                                                           gTCCTCAGTATCCCGAAGATACTGGCCCAACTCATCCACATCAGCCTCGTCGGCACATTCCAGCCAGTCTACCAGTCGGTGTAACGTCTGACTGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGACCGCCCATTTGACATGTGTGTAATCTACTTTGCATCGTATACCACGCTTTACAGCGTGGTGTTTACCTACAACCTCCATCAAACTGAGCCCCTTAGGTAGGTTTTAGCCGAATGAATTGGttttt                                                                                                                                                                                                      >  4:228614/1‑271 (MQ=255)
                                                                                          aaGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGt                                                                                                                                                                                                                            <  4:256433/237‑1 (MQ=255)
                                                                                                                      gCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCAcc                                                                                                                                                                                                                                       >  2:222032/1‑198 (MQ=255)
                                                                                                                          cgtcgGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCg                                                                                                                                                                                                                                                                           >  1:148323/1‑158 (MQ=255)
                                                                                                                                     tCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtg                                                                                                                                                                                                                                                                         <  7:997495/149‑1 (MQ=255)
                                                                                                                                     tCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtg                                                                                                                                                                                                                                                                         <  8:329477/149‑1 (MQ=255)
                                                                                                                                     tCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtg                                                                                                                                                                                                                                                                         <  6:513820/149‑1 (MQ=255)
                                                                                                                                       cAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtggt                                                                                                                                                                                                                                                                       <  5:908648/149‑1 (MQ=255)
                                                                                                                                        aGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtggtg                                                                                                                                                                                                                                                                      >  7:653202/1‑149 (MQ=255)
                                                                                                                                         gCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGgtgt                                                                                                                                                                                                                                                                     <  5:747731/149‑1 (MQ=255)
                                                                                                                                          cAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtggtata                                                                                                                                                                                                                                                                    <  7:540371/149‑4 (MQ=255)
                                                                                                                                           aaGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtggtatac                                                                                                                                                                                                                                                                   <  7:363356/149‑5 (MQ=255)
                                                                                                                                              tCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCtt                                                                                                                                                                                                 >  2:340111/1‑212 (MQ=255)
                                                                                                                                                tACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTa                                                                                                                                                                                                                                                              >  7:1143280/1‑149 (MQ=255)
                                                                                                                                                     gTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtcc                                                                                                                                                                                                                                                         >  6:454453/1‑149 (MQ=255)
                                                                                                                                                      tCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGc                                                                                                                                                                                                                                                                            <  6:226401/129‑1 (MQ=255)
                                                                                                                                                          tgtAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCTAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCa                                                                                                                                           <  4:369081/254‑1 (MQ=255)
                                                                                                                                                            tAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTg                                                                                                                                   <  2:98681/260‑1 (MQ=255)
                                                                                                                                                               cGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGc                                                                                                                                                                                                                                               <  5:1651742/149‑1 (MQ=255)
                                                                                                                                                                gTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTAc                                                                                                                                                                                                                                                                               >  2:608480/1‑116 (MQ=255)
                                                                                                                                                                gTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGa                                                                                                                                                                                   >  4:115097/1‑208 (MQ=255)
                                                                                                                                                                gTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtggta                                                                                                                                                                                                                                                                      >  5:18826/1‑124 (MQ=255)
                                                                                                                                                                 tCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGttt                                                                                                                                                                                                                                                                   >  6:1272548/1‑127 (MQ=255)
                                                                                                                                                                  cTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTaa                                                                                                                                                                                                                                                                 <  7:702377/128‑1 (MQ=255)
                                                                                                                                                                  cTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGa                                                                                                                                                                                                                                            <  6:163069/149‑1 (MQ=255)
                                                                                                                                                                    gAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGaa                                                                                                                                                                        >  2:552351/1‑215 (MQ=255)
                                                                                                                                                                     aaTGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtg                                                                                                                                                                                                                                                                         <  8:499135/117‑1 (MQ=255)
                                                                                                                                                                       tGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtg                                                                                                                                                                                                                                                                         >  7:867271/1‑115 (MQ=255)
                                                                                                                                                                         ggAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtggt                                                                                                                                                                                                                                                                       <  7:922249/115‑1 (MQ=255)
                                                                                                                                                                           aaTGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggta                                                                                                                                                                                                                                <  1:63542/152‑1 (MQ=255)
                                                                                                                                                                           aaTGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtag                                                                                                                                                                                                                                   <  8:1639554/149‑1 (MQ=255)
                                                                                                                                                                              gCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggta                                                                                                                                                                                                                                >  7:1585712/1‑149 (MQ=255)
                                                                                                                                                                                cAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                                                                                                                                                              <  6:758160/149‑1 (MQ=255)
                                                                                                                                                                                 aaTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAggg                                                                                                                                                                                                                             >  6:1408204/1‑149 (MQ=255)
                                                                                                                                                                                 aaTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAggg                                                                                                                                                                                                                             >  6:1124604/1‑149 (MQ=255)
                                                                                                                                                                                  aTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                                                                                                                                                              <  8:279578/147‑1 (MQ=255)
                                                                                                                                                                                      cttctAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtggtg                                                                                                                                                                                                                                                                      <  7:438497/103‑1 (MQ=255)
                                                                                                                                                                                       ttctAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGc                                                                                                                                                                                                                       <  5:418305/149‑1 (MQ=255)
                                                                                                                                                                                          tAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCa                                                                                                                                                                                                                                         <  5:240153/128‑1 (MQ=255)
                                                                                                                                                                                          tAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGa                                                                                                                                                                                                                    >  5:1641394/1‑149 (MQ=255)
                                                                                                                                                                                                aTTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTAt                                                                                                                                                                                                              >  5:471040/1‑149 (MQ=255)
                                                                                                                                                                                                aTTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTAt                                                                                                                                                                                                              <  7:272085/149‑1 (MQ=255)
                                                                                                                                                                                                 ttGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTTCATAATCGATCTGATTTTTCAGATCTTCCGAATAACTTTCAGCACTTAAAt                                                                                >  2:409490/1‑274 (MQ=255)
                                                                                                                                                                                                   gACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGaa                                                                                                                                                                                                                   <  8:6664/141‑1 (MQ=255)
                                                                                                                                                                                                     cacCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCa                                                                                                                                                                                               <  2:745739/159‑1 (MQ=255)
                                                                                                                                                                                                        cTGGTGGTAGATGGAAACTTTCGAACGCCAATATGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCa                                                                        >  4:55803/1‑275 (MQ=255)
                                                                                                                                                                                                         tggtggTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGa                                                                                                                                                                                                     <  8:1343307/149‑1 (MQ=255)
                                                                                                                                                                                                          ggtggtAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGAGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCATCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTCCAAATCATTGTTGGGGGCAATGGAAAGCGTACGTTTTTCCCCt                                                                                                                                                      >  4:688295/1‑195 (MQ=255)
                                                                                                                                                                                                            tggtAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCAcc                                                                                                                                                                                                                                       <  2:463737/112‑1 (MQ=255)
                                                                                                                                                                                                                   tGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCaaa                                                                                                                                                                                             >  8:370283/1‑147 (MQ=255)
                                                                                                                                                                                                                    ggAATCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATATCCGAAGGTATTGGTTATGAGCTTCAAATCa                                                                                                                                                                                          >  8:910941/1‑149 (MQ=255)
                                                                                                                                                                                                                    ggAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTa                                                                                                                                                                                                                                                                                >  7:1454346/1‑63 (MQ=255)
                                                                                                                                                                                                                    ggAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTaa                                                                                                                                                                                                                                                                 <  7:330884/78‑1 (MQ=255)
                                                                                                                                                                                                                    ggAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCa                                                                                                                                                                                               <  8:1410788/144‑1 (MQ=255)
                                                                                                                                                                                                                    ggAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCa                                                                                                                                                                                          >  5:709349/1‑149 (MQ=255)
                                                                                                                                                                                                                     gAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCaa                                                                                                                                                                                              >  4:604752/1‑144 (MQ=255)
                                                                                                                                                                                                                      aaaCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCGCCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAtt                                                                                                                                                                                        >  7:388263/1‑149 (MQ=255)
                                                                                                                                                                                                                      aaaCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAtt                                                                                                                                                                                        >  5:1097358/1‑149 (MQ=255)
                                                                                                                                                                                                                       aaCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAttg                                                                                                                                                                                       >  5:1587023/1‑149 (MQ=255)
                                                                                                                                                                                                                          tttCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAttgttg                                                                                                                                                                                    >  5:834260/1‑149 (MQ=255)
                                                                                                                                                                                                                          tttCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAttgttg                                                                                                                                                                                    >  6:1585111/1‑149 (MQ=255)
                                                                                                                                                                                                                           ttCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGa                                                                                                                                                                                   <  6:620474/149‑1 (MQ=255)
                                                                                                                                                                                                                           ttCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGa                                                                                                                                                                                   <  8:1079938/149‑1 (MQ=255)
                                                                                                                                                                                                                             cGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAgg                                                                                                                                                                                 <  5:1374879/149‑1 (MQ=255)
                                                                                                                                                                                                                             cGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAgg                                                                                                                                                                                 >  8:1550639/1‑149 (MQ=255)
                                                                                                                                                                                                                              gAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAggg                                                                                                                                                                                <  5:1230090/149‑1 (MQ=255)
                                                                                                                                                                                                                                   ccAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATg                                                                                                                                                                           <  8:185411/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  tgtAATCGACTTTGCTGCGTATACCACNCTTTACAGNNTGGTGTTTAACTACATCCTCCATNAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCa                                                                                                                                           <  3:63387/166‑1 (MQ=255)
                                                                                                                                                                                                                                                      aTCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGaaa                                                                                                                                                                       >  8:776109/1‑134 (MQ=255)
                                                                                                                                                                                                                                                      aTCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTccc                                                                                                                                                        >  7:1377178/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      aTCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTccc                                                                                                                                                        <  8:505018/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       tCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAg                                                                                                                                                                                                                  <  5:1159584/90‑1 (MQ=255)
                                                                                                                                                                                                                                                       tCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAg                                                                                                                                                                                                                  >  6:1159584/1‑90 (MQ=255)
                                                                                                                                                                                                                                                          aCTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTc                                                                                                                                                    <  6:1315491/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           cTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCa                                                                                                                                                   <  7:1689995/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           cTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCa                                                                                                                                                   <  8:12898/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               gctgcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGaaa                                                                                                                                               <  8:8857/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tgcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATc                                                                                                                                             >  7:1282824/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   cGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCa                                                                                                                                           <  7:502091/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tataCCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCCCTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTGCGACTGAACGTATGGGGAAc       >  1:107723/1‑279 (MQ=255)
                                                                                                                                                                                                                                                                      acaCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAgg                                                                                                                                                                                 >  8:1193048/3‑108 (MQ=255)
                                                                                                                                                                                                                                                                       caCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAtt                                                                                                                                                                                        <  6:992161/99‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCATGCTTAGCACCTTAGGTAGGGTAAAGCCGAAGGTATTGGGTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGGACGTTTCTCCCCTTCAGAAATCCACTGTCGt                                                                                                                                    >  7:1042621/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGt                                                                                                                                    >  8:438131/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTg                           >  2:664367/1‑254 (MQ=255)
                                                                                                                                                                                                                                                                           aCGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                                                                                                                                                              >  5:340498/1‑58 (MQ=255)
                                                                                                                                                                                                                                                                           aCGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                                                                                                                                                              <  6:340498/58‑1 (MQ=255)
                                                                                                                                                                                                                                                                           aCGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTc                                                                                                                                                                                                <  6:1655702/88‑1 (MQ=255)
                                                                                                                                                                                                                                                                           aCGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTa                                                                                                                                                                  >  6:865786/1‑118 (MQ=255)
                                                                                                                                                                                                                                                                             gCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGaaa                                                                                                                                               >  2:18358/1‑135 (MQ=255)
                                                                                                                                                                                                                                                                               tttACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAtt                                                                                                                                                                                        >  5:1337518/1‑92 (MQ=255)
                                                                                                                                                                                                                                                                               tttACAGCGGGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTGCGACTGAACGTATGGGGAACCAGTc  <  2:110904/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                
AAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTGCGACTGAACGTATGGGGAACCAGTC  >  NC_000913/2761348‑2761892

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