Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,364,860 G→A Y61Y (TAC→TAT dsbD ← thiol‑disulfide exchange protein DsbD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,364,8600GA98.5% 237.7 / NA 68Y61Y (TAC→TATdsbDthiol‑disulfide exchange protein DsbD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (28/39);  minor base C (1/0);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGG  >  NC_000913/4364615‑4365104
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
gCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATagag                                                                                                                                                                                                                                                 <  1:162043/251‑1 (MQ=255)
 cGACCACTTCGCTTACCGGAACGGTTTTGGTTTCTGGCGGATCACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGGGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAgg                                                                                                                                                                                                                                                >  2:382786/1‑251 (MQ=255)
     cACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGgtagt                                                                                                                                                                                                                                            <  1:123056/251‑1 (MQ=255)
      aCTTCGCGTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGTTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCATTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGgtagta                                                                                                                                                                                                                                           <  2:651337/251‑1 (MQ=255)
               aacggaacggTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCTCTCGCCTGGTTTATGGTGACGGGAAGCGTTAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCTTCttt                                                                                                                                                                                                                                  <  2:46244/251‑1 (MQ=255)
                     acggTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTg                                                                                                                                                                                                                            >  2:654210/1‑251 (MQ=255)
                       ggTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGcc                                                                                                                                                                                                                          <  2:9058/251‑1 (MQ=255)
                         ttttGGTTTCTGGCGGATAACAGAAACAGGCATCAGCACAGCCCTGGCAGGTGAAAGGTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGGAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGGAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTgac                                                                                                                                                                                                                          >  1:568845/1‑247 (MQ=255)
                            tattttCTGACGGCTACCAGAACCCGGCATCAGCACAGCCATGGTCGGGGACAGTTAACGTCGCTCCCGCACTCGCATGGTTGATGGTCACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTCAAACTCATCTTCCTACCAGACTCCTTGCGGCAGCTGCAAGGCGGCAATCTTCGCGTGTTCCGGCGAAATGAGGATCTGTTTACGCTAGAGGTAGTAACAGTCTTTGATCTGCCAggg                                                                                                                                                                                                                      >  2:633087/4‑249 (MQ=255)
                                    ggCGGATAACAGAAACCGGTATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTaa                                                                                                                                                                                                              <  1:47476/250‑1 (MQ=255)
                                    ggCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTa                                                                                                                                                                                                               <  2:97768/249‑1 (MQ=255)
                                    ggCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGAGCTGCCAGGTCAGATTa                                                                                                                                                                                                               <  2:612771/249‑1 (MQ=255)
                                        gATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCa                                                                                                                                                                                                         <  1:579007/251‑1 (MQ=255)
                                           aaCAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGTACTCGGCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTATTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGt                                                                                                                                                                                                      <  2:462607/251‑1 (MQ=255)
                                                aaaCCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGtt                                                                                                                                                                                                 >  2:185752/1‑251 (MQ=255)
                                                  aCCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGtttt                                                                                                                                                                                               <  2:556054/251‑1 (MQ=255)
                                                   ccGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTg                                                                                                                                                                                              <  2:371559/251‑1 (MQ=255)
                                                   ccGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGACCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTg                                                                                                                                                                                              <  2:108408/251‑1 (MQ=255)
                                                         tCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCATGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCACTTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCATGTCCGATTAAGGTCATGTTGGTTTTGCTGaa                                                                                                                                                                                         >  2:122561/1‑250 (MQ=255)
                                                                      ggtaggtaACAGTTCACGTCGCTCCCGCCCGCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAACTCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTACCCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGaaaa                                                                                                                                                                            <  2:375321/250‑1 (MQ=255)
                                                                         aGGNGACAGTTAACGTCGCTCCCGCACTCGACTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCt                                                                                                                                                                        <  2:590348/251‑1 (MQ=255)
                                                                          ggTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTg                                                                                                                                                                                           >  2:392961/1‑231 (MQ=255)
                                                                          ggTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTg                                                                                                                                                                                           <  1:392961/231‑1 (MQ=255)
                                                                            tGACAGTTAACGTCGCTCCCACACTCGCATGGTTGATGGTGACGGGAAGAGTCAGCCGAACGCGGTAAATCTCGCATTTGCCGGAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCACTTTTCGCGTGTTCCGGCGTAATGCGGATCTGGTTACGATAGAGGGAGTAACCGTCTTTGATCTGCCAGGTCAGATTACGGTCATGTTGGGTTTGCTGAAATTCAAAAGCCAAGGcatg                                                                                                                                                                      >  2:647307/1‑247 (MQ=255)
                                                                             gACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGa                                                                                                                                                                     <  1:607425/250‑1 (MQ=255)
                                                                                   tAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGTTAAATCTCTCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCcgcg                                                                                                                                                               <  2:488401/250‑1 (MQ=255)
                                                                                   tAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCcgcg                                                                                                                                                               <  2:514547/250‑1 (MQ=255)
                                                                                        tCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCCTGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGAc                                                                                                                                                          >  2:189635/1‑250 (MQ=255)
                                                                                                     cGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACg                                                                                                                                              <  2:200639/249‑1 (MQ=255)
                                                                                                      gCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTcc                                                                                                                                           >  2:44143/1‑251 (MQ=255)
                                                                                                           gTTGATGGTGACGGGAAGCGTCAGCCGAGAGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCAAGACGCCTTGCGGCAGATGCACGTCGGCAATTTTGGCGTGTTCCGGCGTAATGCGGATCTGCTTACGATAGAGGTAGTAACCGTCTTTGATCTGACAGGTAAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTACCGgcg                                                                                                                                      >  2:223975/1‑251 (MQ=255)
                                                                                                                gggTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCGTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCg                                                                                                                                 <  1:189635/250‑1 (MQ=255)
                                                                                                                    gACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGaata                                                                                                                             <  1:415775/251‑1 (MQ=255)
                                                                                                                       gggAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCTCGTCGAATAATc                                                                                                                          <  2:302324/251‑1 (MQ=255)
                                                                                                                       gggAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGCCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATc                                                                                                                          >  2:512516/1‑251 (MQ=255)
                                                                                                                          aaGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCTTCGAATAATCCgg                                                                                                                       <  1:122561/251‑1 (MQ=255)
                                                                                                                                      aTCGCGGTCAATCTCGCTTTCGCCGTAAAACTCATCTTCATGACAGCCGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAATCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGTGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAACCGGAAg                                                                                                           <  2:300058/251‑1 (MQ=255)
                                                                                                                                      aTCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAg                                                                                                           <  1:602044/251‑1 (MQ=255)
                                                                                                                                         gcgGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGGCGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGTACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGtgc                                                                                                        <  1:399210/251‑1 (MQ=255)
                                                                                                                                          cgGTAAATCTCGCTTTCGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTTGACGTCGGCAATTTTCTCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCGTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAACCGGAAGtgct                                                                                                       <  1:505315/251‑1 (MQ=255)
                                                                                                                                           ggTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGCCAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGtgctg                                                                                                      >  2:326038/1‑251 (MQ=255)
                                                                                                                                                          tGCCGTAAAACTCATCTTCATGCCAGACGCCTTCCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGCTTAAGGTCATGTTGGTTTTGCTGAAAATCAACAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAATGTCCCGGCGCGTCGAATAATCAGGCAAAAATGGAAGTGCTGCAAAGTAGCAGGATc                                                                                       >  1:451941/1‑251 (MQ=255)
                                                                                                                                                             cGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGc                                                                                    <  1:504321/251‑1 (MQ=255)
                                                                                                                                                             cGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGc                                                                                    >  2:252604/1‑251 (MQ=255)
                                                                                                                                                                aaaaCTCATCTTCATGCCAGACGCCTGGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGt                                                                                  <  1:60761/250‑1 (MQ=255)
                                                                                                                                                                 aaaCTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTa                                                                                 >  2:29135/1‑250 (MQ=255)
                                                                                                                                                                   aCTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTaa                                                                                >  2:327593/1‑249 (MQ=255)
                                                                                                                                                                      cATCTTCATACCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGCATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAANATg                                                                           >  2:609067/1‑251 (MQ=255)
                                                                                                                                                                                  aGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATg                                                               >  2:9777/1‑251 (MQ=255)
                                                                                                                                                                                  aGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATg                                                               >  1:104068/1‑251 (MQ=255)
                                                                                                                                                                                      gCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGTAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGagag                                                           <  1:259754/251‑1 (MQ=255)
                                                                                                                                                                                                     cACGTCGGCAATTTGCACGTGTTCCGTCATAATTCGGATCTGTGTACGATAGAGGTATTAACCGTCTGTGATCGGCCAGGTCGGATTAAGGGCATGTTGGTTTTGAGGAAAATCAAAAGCAAAGGCTTGATCCGCGGAGACAAATTGTGAACGTCCCGGCTGGTCGAATAATGCGACAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTAAGCCATGAGAAGTAATCTGTGTCTcc                                            <  2:343716/251‑1 (MQ=255)
                                                                                                                                                                                                          cGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCgtgt                                        <  2:579478/250‑1 (MQ=255)
                                                                                                                                                                                                               aTTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACATATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCGAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGCTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAAc                                   >  1:457445/1‑250 (MQ=255)
                                                                                                                                                                                                                     gcgTGTTCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTaa                             >  1:269684/1‑250 (MQ=255)
                                                                                                                                                                                                                      cgTGTTCCGGCGTAATGCGGATATGCTTACGATCGAGGTAGTAACCGTCTTTGAGCTGCCAGGTCATGTTAAGGTCATGTTGGTTTTGCTGAAAATCACAAGCAAAGGCTTGATCCGCGGGGACCAATTGTGAACGTCCCGGAGCGCCGAATAATCCGACAAAAACGGAAGTGCTGAAAAGTAGCAGGATCAGCGTAAAGCTGCGTTTAGCCATGAGAGGTAATCTGTGTCTACGTGTGTAACAGGGaaa                            >  2:287459/1‑250 (MQ=255)
                                                                                                                                                                                                                          ttCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAAcc                         >  2:138581/1‑249 (MQ=255)
                                                                                                                                                                                                                           tCCGGCGTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAg                       <  1:543407/250‑1 (MQ=255)
                                                                                                                                                                                                                             cGGCATAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTTATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCCAATAATCCGGCAAAAGCGGAAGTGCTGCAAAGTAGCAGTATCTGCGTAAAGATGCGTTTAGCCATGAGAGGTAATCTGTGTCCCCGTGTGGAACAGGTAAAACCAGAAg                    <  2:101672/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gTAATGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTcc                >  2:185237/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aTGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATgaga                                                            >  1:627069/1‑204 (MQ=255)
                                                                                                                                                                                                                                    aTGCGGATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATgaga                                                            <  2:627069/204‑1 (MQ=255)
                                                                                                                                                                                                                                        ggATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAGCGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCgtgt                                        <  2:104429/220‑1 (MQ=255)
                                                                                                                                                                                                                                        ggATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAGCGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCgtgt                                        >  1:104429/1‑220 (MQ=255)
                                                                                                                                                                                                                                        ggATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCATAAGTTCCGGGGttt         <  1:309759/251‑1 (MQ=255)
                                                                                                                                                                                                                                        ggATCTGTTTACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGttt         >  2:518803/1‑251 (MQ=255)
                                                                                                                                                                                                                                                ttACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATgg  <  1:327593/250‑1 (MQ=255)
                                                                                                                                                                                                                                                ttACGATAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATgg  >  1:194194/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
GCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGG  >  NC_000913/4364615‑4365104

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