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→GCA) recQATP‑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

ACAGGTTTTACAAGAAACCTTTGGCTACCAACAGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCGGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGC‑CACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCC  >  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)
                     tGGCTACCAACAGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AAc                                                                                                                                                                                                                >  1:162999/1‑251 (MQ=255)
                     tGGCTACCAACAGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AAc                                                                                                                                                                                                                >  2:326017/1‑251 (MQ=255)
                       gCTACCAACAGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACtt                                                                                                                                                                                                              >  2:317914/1‑251 (MQ=255)
                        cTACCAACAGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTg                                                                                                                                                                                                             >  1:261519/1‑251 (MQ=255)
                                aGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGCTCAGGTGGATCAACTGCAAGCCACCGGCGTGGCAGCGGCGGGCCTTAACTCGCCGCAAACCCGCGAACAGC‑AACTTGAAGtgatg                                                                                                                                                                                                     >  2:180386/1‑251 (MQ=255)
                                 gTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGtgatga                                                                                                                                                                                                    <  1:86832/251‑1 (MQ=255)
                                  tttCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGAc                                                                                                                                                                                                   >  2:127195/1‑251 (MQ=255)
                                    tCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACa                                                                                                                                                                                                  <  1:100744/250‑1 (MQ=255)
                                      gCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCAGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAACTCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCAAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAAAAGC‑AACTTGAAGGGATGACAgg                                                                                                                                                                                                >  1:370678/1‑250 (MQ=255)
                                                ggAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑CACTTGAAGTGATGACAGGCTGCCGCACCg                                                                                                                                                                                     >  2:284852/1‑251 (MQ=255)
                                                 gaagaaATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCgg                                                                                                                                                                                    <  2:115016/251‑1 (MQ=255)
                                                       aTTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATACCCACTGTTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCTTTATTGCTAAACGGCCTTACCGTGGTTGGTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCCACTGCAAGCCAACGGCGTGGCAGCGACGTGCCTTCACTCGACACAAACCCGCGAACCGC‑AACTTGAAGTGATGACAGGCTGCCGCAcctgacaact                                                                                                                                                                              >  1:274206/1‑243 (MQ=255)
                                                                    tgCTTTCCGGCCGCGATTGCCTCGTCGTCCTGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGCTCAGGCGGCTCAACTGCAAGCCAACGGCGTGGCAGCGGCGCGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTtat                                                                                                                                                                 >  1:340423/1‑251 (MQ=255)
                                                                    ggCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCCCCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTtat                                                                                                                                                                 <  2:106113/250‑1 (MQ=255)
                                                                           cGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATCTCGccc                                                                                                                                                          <  1:51518/251‑1 (MQ=255)
                                                                              ccGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCgg                                                                                                                                                       <  1:333767/251‑1 (MQ=255)
                                                                                   aTTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAAcac                                                                                                                                                  >  1:301986/1‑249 (MQ=255)
                                                                                   aTTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGc                                                                                                                                                  <  1:284852/251‑1 (MQ=255)
                                                                                           gtcgtcATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGACTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGACAACCCGCGCACAGC‑AACTTGAAGTGATGACAGGCCGCCGCACCGGGCAAATTCGTCTGCTTTCTATCGCCCCGGAACGCCTGATGCt                                                                                                                                          >  1:170160/1‑251 (MQ=255)
                                                                                                    cccACTGGTGGCGTAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACtt                                                                                                                                 <  2:285366/251‑1 (MQ=255)
                                                                                                                      tCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGgcgc                                                                                                                >  1:39159/1‑250 (MQ=255)
                                                                                                                       cccTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGACTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCAc                                                                                                              >  2:57531/1‑251 (MQ=255)
                                                                                                                           ttGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCATGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCTACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCATATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCTCACTgg                                                                                                           <  2:379203/250‑1 (MQ=255)
                                                                                                                           ttGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAAGGCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTg                                                                                                            <  1:57053/249‑1 (MQ=255)
                                                                                                                             gCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGCTCAGGTGGATCACCTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAAt                                                                                                        >  1:193114/1‑251 (MQ=255)
                                                                                                                                aTCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATcc                                                                                                      <  2:121364/250‑1 (MQ=255)
                                                                                                                                    aaTCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCCGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGAGTGCCTTAACTCGACGCAAACCCGCGAAAAGC‑CACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGCTGCTGGATAACTTTCTTGAGCCTCTGGCGCACTGGCATCCGgtgt                                                                                                 >  1:247877/1‑251 (MQ=255)
                                                                                                                                         cTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACGGGAATCCGGTGttatt                                                                                             >  1:295036/1‑250 (MQ=255)
                                                                                                                                            ccTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGcc                                                                                         >  1:60867/1‑251 (MQ=255)
                                                                                                                                              ttATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGGTGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGCCTCCGGTGTTATTAGCCGt                                                                                       >  2:40906/1‑251 (MQ=255)
                                                                                                                                                  tGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGCTGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCGAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTtgat                                                                                   >  1:43047/1‑251 (MQ=255)
                                                                                                                                                  tGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTtgat                                                                                   >  1:84127/1‑251 (MQ=255)
                                                                                                                                                                  cGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTACCTCGACGCAAACCCGCTAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGTAATCCGGTGTTATTATCCGTTTATGAAGCGCACTGTAtct                                                                   <  2:261519/251‑1 (MQ=255)
                                                                                                                                                                      gttgttTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCGTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGGTATTAGCCGTTGATGAAGCGCACTGTATCTCCCa                                                               <  2:48853/251‑1 (MQ=255)
                                                                                                                                                                                         cGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGC‑CACGATTTCCg                                            <  1:257820/250‑1 (MQ=255)
                                                                                                                                                                                              atgaAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGCTGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGGGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGC‑CACGGTTTCCGCCCgg                                       >  2:292809/1‑250 (MQ=255)
                                                                                                                                                                                                 aaaGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGC‑CACGATTTCCGCCCGGAat                                    >  1:162548/1‑250 (MQ=255)
                                                                                                                                                                                                            ggATCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGC‑CACGATTTCCGCCCGGAATATGCCGCGCTc                         <  1:326017/250‑1 (MQ=255)
                                                                                                                                                                                                              aTCAACTGCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGC‑CACGATTTCCGCCCGGAATATGCCGCGCTCGGt                      <  1:254195/251‑1 (MQ=255)
                                                                                                                                                                                                                     gCAAGCCAACGGCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGC‑CACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTg                <  2:89524/250‑1 (MQ=255)
                                                                                                                                                                                                                            aaCGGCGTGGCCGCGGCGGGACTTAACTCTACACCAAACCACCAACAGCCAACT‑GAAGAGATGACAGGCTGCCGCACCGGGCAAATTCATCAGCTTTCTATCGCCCCCGAACGCCTGATGCTGGATAACTTTCTTGAGAACCTGGCGCACTGGAATCCGGGGTTATTTGCCCCTAATGGAGCGCACTGTATCTCCCAATGGGGCTCA‑GCTTTCCGCCAGGAATATGACGCCCGCGGGCAGTTTCGcccagg        >  2:70637/1‑247 (MQ=255)
                                                                                                                                                                                                                               ggCGTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGC‑CACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGtt     >  1:326175/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gTGGCAGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGC‑CACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTccc  <  1:57531/251‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
ACAGGTTTTACAAGAAACCTTTGGCTACCAACAGTTTCGCCCCGGCCAGGAAGAAATTATCGACACTGTGCTTTCCGGCCGCGATTGCCTCGTCGTCATGCCCACTGGTGGCGGAAAATCCCTTTGCTATCAAATCCCTGCCTTATTGCTAAACGGCCTTACCGTGGTTGTTTCACCGCTGATTTCGTTGATGAAAGATCAGGTGGATCAACTGCAAGCCAACGGCGTGGCGGCGGCGTGCCTTAACTCGACGCAAACCCGCGAACAGC‑AACTTGAAGTGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGC‑CACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCC  >  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: ‑