Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,006,136 G→A A91A (GCG→GCA recQ → ATP‑dependent DNA helicase RecQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,006,1360GA98.0% 169.9 / NA 50A91A (GCG→GCArecQATP‑dependent DNA helicase RecQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (29/20);  minor base C (1/0);  total (30/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ACAGGTTTTACAAGAAACCTTTGGCTACCAACAGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCGGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCC  >  NC_000913/4005905‑4006381
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
aCAGGTTTTACAAGAAACCTTTGGCTACCAACAGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCg                                                                                                                                                                                                                                       >  1:139329/1‑250 (MQ=255)
  aGGTTTTGCAAGAAACCTTTGGATACCAACAGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCGCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTAGCAAATCCCTGCCTTATTGCTAAACGGCCTGACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTTGAGCTACTGCAAGCCAACGGCTTGGCAGCGGCGTGCCTTAACTCGACg                                                                                                                                                                                                                                    <  1:180596/251‑1 (MQ=255)
   ggTTTTACAAGAAACCTTTGGCTACCAACAGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGc                                                                                                                                                                                                                                   >  1:239495/1‑251 (MQ=255)
   ggTTTTACAAGAAAACTTTGGCTACCAACAGTTTCGCCCCGGCCAGGAAGAAATTATAGACACTGTGCTTTCCGGCCGCGATTGCCTCGTAGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGCTTTCGTTGATGCAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTTGACGc                                                                                                                                                                                                                                   >  2:304866/1‑251 (MQ=255)
         aCAAGAAACCTTTGGCTACCAACAGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAAcc                                                                                                                                                                                                                              <  1:191586/250‑1 (MQ=255)
                aCCTTTGGCTACCAACAGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCCAACCCGCGAAc                                                                                                                                                                                                                       >  1:292159/1‑250 (MQ=255)
                  cTTTGGCTACCAACAGTTTCGCCCCGGCCAGGCAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCACACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTACACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGCAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGc                                                                                                                                                                                                                    >  1:344961/1‑251 (MQ=255)
                     tGGCTACCAACAGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAAc                                                                                                                                                                                                                >  1:162999/1‑251 (MQ=255)
                     tGGCTACCAACAGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAAc                                                                                                                                                                                                                >  2:326017/1‑251 (MQ=255)
                       gCTACCAACAGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACtt                                                                                                                                                                                                              >  2:317914/1‑251 (MQ=255)
                        cTACCAACAGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTg                                                                                                                                                                                                             >  1:261519/1‑251 (MQ=255)
                                aGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGCTCAGGTGGATCAACTGCAAGCCACCGGCGTGGCAGCGGCGGGCCTTAACTCGCCGCAAACCCGCGAACAGCAACTTGAAGtgatg                                                                                                                                                                                                     >  2:180386/1‑251 (MQ=255)
                                 gTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGtgatga                                                                                                                                                                                                    <  1:86832/251‑1 (MQ=255)
                                  tttCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGAc                                                                                                                                                                                                   >  2:127195/1‑251 (MQ=255)
                                    tCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACa                                                                                                                                                                                                  <  1:100744/250‑1 (MQ=255)
                                      gCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCAGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAACTCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCAAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAAAAGCAACTTGAAGGGATGACAgg                                                                                                                                                                                                >  1:370678/1‑250 (MQ=255)
                                                ggAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCCACTTGAAGTGATGACAGGCTGCCGCACCg                                                                                                                                                                                     >  2:284852/1‑251 (MQ=255)
                                                 gaagaaATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCgg                                                                                                                                                                                    <  2:115016/251‑1 (MQ=255)
                                                       aTTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATACCCACTGTTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCTTTATTGCTAAACGGCCTTACCGTGGTTGGTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCCACTGCAAGCCAACGGCGTGGCAGCGACGTGCCTTCACTCGACACAAACCCGCGAACCGCAACTTGAAGTGATGACAGGCTGCCGCAcctgacaact                                                                                                                                                                              >  1:274206/1‑243 (MQ=255)
                                                                    tgCTTTCCGGCCGCGATTGCCTCGTCGTCCTGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGCTCAGGCGGCTCAACTGCAAGCCAACGGCGTGGCAGCGGCGCGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTtat                                                                                                                                                                 >  1:340423/1‑251 (MQ=255)
                                                                    ggCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCCCCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTtat                                                                                                                                                                 <  2:106113/250‑1 (MQ=255)
                                                                           cGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATCTCGccc                                                                                                                                                          <  1:51518/251‑1 (MQ=255)
                                                                              ccGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCgg                                                                                                                                                       <  1:333767/251‑1 (MQ=255)
                                                                                   aTTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAAcac                                                                                                                                                  >  1:301986/1‑249 (MQ=255)
                                                                                   aTTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGc                                                                                                                                                  <  1:284852/251‑1 (MQ=255)
                                                                                           gtcgtcATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGACTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGACAACCCGCGCACAGCAACTTGAAGTGATGACAGGCCGCCGCACCGGGCAAATTCGTCTGCTTTCTATCGCCCCGGAACGCCTGATGCt                                                                                                                                          >  1:170160/1‑251 (MQ=255)
                                                                                                    cccACTGGTGGCGTAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACtt                                                                                                                                 <  2:285366/251‑1 (MQ=255)
                                                                                                                      tCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGgcgc                                                                                                                >  1:39159/1‑250 (MQ=255)
                                                                                                                       cccTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGACTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCAc                                                                                                              >  2:57531/1‑251 (MQ=255)
                                                                                                                           ttGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCATGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCTACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCATATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCTCACTgg                                                                                                           <  2:379203/250‑1 (MQ=255)
                                                                                                                           ttGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAAGGCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTg                                                                                                            <  1:57053/249‑1 (MQ=255)
                                                                                                                             gCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGCTCAGGTGGATCACCTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAAt                                                                                                        >  1:193114/1‑251 (MQ=255)
                                                                                                                                aTCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATcc                                                                                                      <  2:121364/250‑1 (MQ=255)
                                                                                                                                    aaTCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCCGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGAGTGCCTTAACTCGACGCAAACCCGCGAAAAGCCACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGCTGCTGGATAACTTTCTTGAGCCTCTGGCGCACTGGCATCCGgtgt                                                                                                 >  1:247877/1‑251 (MQ=255)
                                                                                                                                         cTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACGGGAATCCGGTGttatt                                                                                             >  1:295036/1‑250 (MQ=255)
                                                                                                                                            ccTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGcc                                                                                         >  1:60867/1‑251 (MQ=255)
                                                                                                                                              ttATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGGTGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGCCTCCGGTGTTATTAGCCGt                                                                                       >  2:40906/1‑251 (MQ=255)
                                                                                                                                                  tGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGCTGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCGAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTtgat                                                                                   >  1:43047/1‑251 (MQ=255)
                                                                                                                                                  tGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTtgat                                                                                   >  1:84127/1‑251 (MQ=255)
                                                                                                                                                                  cGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTACCTCGACGCAAACCCGCTAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGTAATCCGGTGTTATTATCCGTTTATGAAGCGCACTGTAtct                                                                   <  2:261519/251‑1 (MQ=255)
                                                                                                                                                                      gttgttTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCGTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGGTATTAGCCGTTGATGAAGCGCACTGTATCTCCCa                                                               <  2:48853/251‑1 (MQ=255)
                                                                                                                                                                                         cGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCg                                            <  1:257820/250‑1 (MQ=255)
                                                                                                                                                                                              atgaAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGCTGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGGGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGGTTTCCGCCCgg                                       >  2:292809/1‑250 (MQ=255)
                                                                                                                                                                                                 aaaGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAat                                    >  1:162548/1‑250 (MQ=255)
                                                                                                                                                                                                            ggATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTc                         <  1:326017/250‑1 (MQ=255)
                                                                                                                                                                                                              aTCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGt                      <  1:254195/251‑1 (MQ=255)
                                                                                                                                                                                                                     gCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTg                <  2:89524/250‑1 (MQ=255)
                                                                                                                                                                                                                            aaCGGCGTGGCCGCGGCGGGACTTAACTCTACACCAAACCACCAACAGCCAACTGAAGAGATGACAGGCTGCCGCACCGGGCAAATTCATCAGCTTTCTATCGCCCCCGAACGCCTGATGCTGGATAACTTTCTTGAGAACCTGGCGCACTGGAATCCGGGGTTATTTGCCCCTAATGGAGCGCACTGTATCTCCCAATGGGGCTCAGCTTTCCGCCAGGAATATGACGCCCGCGGGCAGTTTCGcccagg        >  2:70637/1‑247 (MQ=255)
                                                                                                                                                                                                                               ggCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGtt     >  1:326175/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTccc  <  1:57531/251‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
ACAGGTTTTACAAGAAACCTTTGGCTACCAACAGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCGGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGCAACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCC  >  NC_000913/4005905‑4006381

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: