Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,579,050 G→A A199V (GCC→GTC)  yhhW ← quercetin 2,3‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,579,0500GA98.1% 146.7 / NA 53A199V (GCC→GTC) yhhWquercetin 2,3‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (34/17);  minor base C (0/1);  total (35/18)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TTTACAATTATTAGATATACACTAAATCAATTGCATCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGGCTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGTTTTCTTC  >  NC_000913/3578807‑3579282
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                        
tttACAATTATTAGATATACACTAAATCAATTGCATCGCATTGTGCTAACGAAGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCAc                                                                                                                                                                                                                                   >  1:74467/1‑251 (MQ=255)
     aaTTATTAGATATAGACTCAATTAATTGCGTCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGAGACCGGAGCGAGCGAGGAACCCTGCAAACAAGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGTACAGTAACACTTCGCTATCTCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATGGACAGCCCCTCCCTGGCCGCGGCTTTCACGCCGt                                                                                                                                                                                                                              >  2:80458/1‑251 (MQ=255)
                       aaaacAATTGCATCGCATTGTGCGAACGACGCGGCAAACCCCTGGAAGCGTACAGCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTGGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGcc                                                                                                                                                                                                              <  1:67333/245‑1 (MQ=255)
                       aaaTCAATTGCATCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGcc                                                                                                                                                                                                              <  2:76846/249‑1 (MQ=255)
                         ctCAATTGCATCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGAAGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCATTTc                                                                                                                                                                                                          >  1:177226/2‑251 (MQ=255)
                         aTCAATTGCATCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCTCAGTTGCTATCCGCATGAATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTc                                                                                                                                                                                                          >  1:177231/1‑251 (MQ=255)
                         aTCAATTGCATCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTc                                                                                                                                                                                                          >  1:35754/1‑251 (MQ=255)
                                 cATCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGCGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTCTCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCCCCTg                                                                                                                                                                                                  >  2:156722/1‑251 (MQ=255)
                                  aTCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTgg                                                                                                                                                                                                 >  1:247734/1‑251 (MQ=255)
                                      cATTGTGCTAACGACGCGGCAAACCCCCGTAAGCGTACATCAGCACGCGACCGGAGCGAGCGGGGAAGCGGGCAAAAAGGCGGGCTTCAAAGATGCCGTTGAGTTTTAAACGGGCGGCAGATGGAACGGTAACACTTCGCTGGCGCTTTCCGCATGGATGGAGATGGGCTGCTCGTCCCAGTTTGCGAGACCATCGCTGGTTGAGACTTCCACGCGGTTAATGGTGACATTGCCTTTCTCCACCTGGATcc                                                                                                                                                                                             <  2:395072/251‑1 (MQ=255)
                                       atgttgttAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTGCACGCCGTTGATGGTGACATTGCCTTTCACCACCTGGATCCa                                                                                                                                                                                            <  1:341812/244‑1 (MQ=255)
                                       aTTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATcc                                                                                                                                                                                             >  1:87659/1‑250 (MQ=255)
                                                 cgacgCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGt                                                                                                                                                                                  >  1:396727/1‑251 (MQ=255)
                                                         cAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGATTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGACGAGACTTTCAAGCCGTTAATGGGGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGc                                                                                                                                                                          >  1:370213/1‑251 (MQ=255)
                                                                        gTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGACGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCACTCTGATGAACCgc                                                                                                                                                           >  1:265/1‑250 (MQ=255)
                                                                           catcaacccaagACCGGGGCGAGCGAAGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTACACCGGCGGCAGATCGATCAGTAAAACTTCGCTCTCGCGATCCGCATGGGTGGAGATTGCCTGCTCATCCCAGATGTCCTGACCATCCCTGTTCGCGACTTTCCCGCCGTTAATGGTGTCATTGCCTTTCACCACCTTGAGCCTGACGCGACGTTCAGCGGCGATTTGATGCACCGACTg                                                                                                                                                        >  2:135460/12‑251 (MQ=255)
                                                                                      gACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTc                                                                                                                                              >  1:323884/1‑250 (MQ=255)
                                                                                                  cgagGAAGCCTGCAAAAAGGCAGGCTTCAAAGCTGCCGTTGAGTTTTAAACCGGCGGCAGATCGCACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTAATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGc                                                                                                                                 >  2:27910/1‑251 (MQ=255)
                                                                                                   ggggagcccttcAAAAAGGCCGGCTTCCAAGACGCCGTTGCGTTTTAAACCGGCGTCAGATCGAACAGTAACACTTCGCTGTCGCTAACCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGATCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCACACGCGTCGATCAGCGACACTCTGATGCACCGCCTGCTCACCTTTCAGCACCGCCAAGCg                                                                                                                                >  2:287071/12‑251 (MQ=255)
                                                                                                              cAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGAAGCTCGAACAGTAACACTTCGCTGCCACTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGCGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTcc                                                                                                                      >  1:183256/1‑250 (MQ=255)
                                                                                                               aaaaaGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCTAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCAc                                                                                                                                                                                                                                   <  1:194944/140‑1 (MQ=255)
                                                                                                               aaaaaGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCAGCATGGCTGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCAc                                                                                                                                                                                                                                   >  2:194944/1‑140 (MQ=255)
                                                                                                                 aaaGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGTTCTTCCCAGATTGCCGGCCCACCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGTCGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCatat                                                                                                                  <  1:238039/251‑1 (MQ=255)
                                                                                                                         gcTTCAAAGAGGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCATTCTCATCTCAGATTGCCAGACCATCGCTGGTCGAGACTTGCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGCTCCATATCCTGATGc                                                                                                          <  2:91142/251‑1 (MQ=255)
                                                                                                                                     ccGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCGGTCGCTATCCGCATGGATGGAGATTGCCTGCTAATCCCAGATTGCCAGACCATCGCTGGTCGAGACTCTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGaa                                                                                              <  1:288995/251‑1 (MQ=255)
                                                                                                                                       gTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAAcc                                                                                            >  2:67771/1‑251 (MQ=255)
                                                                                                                                            gTTTTAAACCTGCGGCAGATCGAACAGTAACACTTCGCGGTCGTTATCCGAATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGAGCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATAGCCTGATGCACTTTCAATGAACCAtctc                                                                                       <  2:35754/251‑1 (MQ=255)
                                                                                                                                             ttttAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTcg                                                                                      >  1:173188/1‑251 (MQ=255)
                                                                                                                                                    ccGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATcc                                                                               >  1:360392/1‑251 (MQ=255)
                                                                                                                                                     cggcggCAGATCGAACAGTAACGCTTCGCTGCCGCTATCCGCATGGTGGGAGATTGCCTGCTCATCCCAGACTGCCAGACCATCGCTGGTCGAGCCTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCCCCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCg                                                                              <  2:131826/251‑1 (MQ=255)
                                                                                                                                                     cagcggcAGATCGAACAGTAACAATTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTAACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTCCAGTTCCATATCCTGCGGCACTTTCACTGACCCATCTCGCGCATCCg                                                                              >  2:17700/3‑251 (MQ=255)
                                                                                                                                                       gcggcAGATCGAACATTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGc                                                                            <  1:356514/251‑1 (MQ=255)
                                                                                                                                                                  aaCAGTAACCCTTCGCGGTCGCTAGCGGCATGGATGGAGATTGCCTGCGCAGCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAg                                                                 <  2:266508/251‑1 (MQ=255)
                                                                                                                                                                      gTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGt                                                             >  1:54872/1‑251 (MQ=255)
                                                                                                                                                                                 cTGTCGCTATCCGCATGGATGGAGATAGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGTTGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTAc                                                  >  2:273679/1‑251 (MQ=255)
                                                                                                                                                                                                 ggatggaGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACg                                   <  1:27910/250‑1 (MQ=255)
                                                                                                                                                                                                        gaTTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCAt                           >  1:307887/1‑251 (MQ=255)
                                                                                                                                                                                                         aTTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCACTGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCAt                           >  2:172954/1‑250 (MQ=255)
                                                                                                                                                                                                            gCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAAcg                       >  1:169053/1‑251 (MQ=255)
                                                                                                                                                                                                            gCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAAcg                       >  1:145197/1‑251 (MQ=255)
                                                                                                                                                                                                                gCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATACGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACgcggc                   >  2:380898/1‑251 (MQ=255)
                                                                                                                                                                                                                gCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACCTTCATGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCCGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGAACTTTCCATGAACCCTCTCGAGCCTCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGAGACGCTGTTCATCACgcggc                   >  2:258819/1‑251 (MQ=255)
                                                                                                                                                                                                                gCTCATCCCAGATGGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAAGGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCAGCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCTCTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCAGCGAAGCGACGCTGTTCATAAcgctg                    <  2:262646/250‑3 (MQ=255)
                                                                                                                                                                                                                gCTCATCCCAGAGTGCCAGACCATCGCTGGTCGAGACTGTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCATACGCGGCGGTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCTACGCTGTTCATAAcgcg                     <  2:193270/249‑1 (MQ=255)
                                                                                                                                                                                                                gCTCACCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACgcggc                   <  1:156722/251‑1 (MQ=255)
                                                                                                                                                                                                                  tCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGt                 >  1:245687/1‑251 (MQ=255)
                                                                                                                                                                                                                   cATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGt                 >  2:146201/1‑250 (MQ=255)
                                                                                                                                                                                                                    aTCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTa                >  1:148557/1‑250 (MQ=255)
                                                                                                                                                                                                                      cccAGATTGCCAGACCATCGCTGGTCGAGACTTGCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTAATa             <  2:323884/251‑1 (MQ=255)
                                                                                                                                                                                                                             tGCCACACCATCTCTGGCCGGGACTTTCGCGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCATACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTGCAGCAACGCCCGGCGGTACAGTTCCATATCCTGATGCACTTTCAATGGACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCGTGTTCGGCAGCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGttt       <  1:381534/250‑1 (MQ=255)
                                                                                                                                                                                                                                cAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGTTttctt   >  2:80116/1‑251 (MQ=255)
                                                                                                                                                                                                                                cAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGTTttctt   >  1:264312/1‑251 (MQ=255)
                                                                                                                                                                                                                                   aCCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGTTttcttc  >  1:61747/1‑249 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                        
TTTACAATTATTAGATATACACTAAATCAATTGCATCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGGCTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGTTTTCTTC  >  NC_000913/3578807‑3579282

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