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,0500GA95.8% 216.5 / NA 72A199V (GCC→GTC) yhhWquercetin 2,3‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (41/28);  minor base C (1/2);  total (42/30)
Fisher's exact test for biased strand distribution p-value = 5.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATTAGATATACACTAAATCAATTGCATCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGGCTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGTTTTCTTCAGGCATGATCCAGA  >  NC_000913/3578814‑3579296
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
ttattaGATATACACTAAATCAATTGCATCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGtt                                                                                                                                                                                                                                           >  1:398522/1‑250 (MQ=255)
 tagtaGATATACACTAAATCAATTGCATCGCATTGTGCTAACGACGCGGCAAACCCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCACCCCAGATGGCCAGCCCATCGCTGGTCGAGACTTTCACGCCGTTa                                                                                                                                                                                                                                          <  1:462470/247‑1 (MQ=255)
  attaGATATACACTAAATCACTTGCATCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCACAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGttcat                                                                                                                                                                                                                                        >  2:122956/1‑248 (MQ=255)
         ataCACTAAATCAATTGCATCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACa                                                                                                                                                                                                                                 <  1:100082/251‑1 (MQ=255)
            cacTAAATCAATTGCATCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTg                                                                                                                                                                                                                              >  1:357496/1‑251 (MQ=255)
                 aatcaaccgcatccccttctgCTAACGACGCGGCAAAACACCGGAACCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAACAAGGCCGGCTTCACACCTGCCGTTGCGTTTTCCACCGACGAAAGATCGACCAGTAACACTTCGCTCTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGCGCCTTTCACGCCGTTAATGGTGAcctttccttt                                                                                                                                                                                                                         >  2:196937/20‑251 (MQ=255)
                  aTCAATTGCATCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTc                                                                                                                                                                                                                        >  1:237428/1‑251 (MQ=255)
                  aTCAATGGCATCGAATTGTTCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTTTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTc                                                                                                                                                                                                                        <  2:68/251‑1 (MQ=255)
                   tCAATTGCATCGCATTGTGCTAACGACGCGGAAAAACCACGGTAGCGTACATCCGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTAGCTGTAGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGAATTTCACGCCGTTAATGGTGACATTGCCTTTc                                                                                                                                                                                                                        >  2:64124/1‑250 (MQ=255)
                    cAATGGCATCGCATTGTGCTAACGACGCGCCAAAACCCCGGAAGCGTCCATCAGCGCGCGACCGAAGCGAGCGAGGAAGCTTGCAAAAGGGCAGGCTTCAAATATGCCGTTGAGTTTTAAACCGGCGCCAGAGCGAAAATTAACCCCTCGCTGTCGCTATCCGCATGGATGGAGATTGCCGGCGCATCCCAGATTGCCAGTCCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTca                                                                                                                                                                                                                       <  1:60224/250‑1 (MQ=255)
                        tGCATCGCATTGTGCTAACGACGCGGCAAACCCCCGGAAGCGTACACCAGCACGGGCCCGGAGCGCGCGAGGAAGCTGGCAAAAAGGCAGGCTGAAAAGAGGCCGTTGAGTTTTAAGCCGGCGGCAGAGCGAAGAGTAACTCTTCGCTGTCGCTATCCGCTGGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTcaccac                                                                                                                                                                                                                   <  1:206718/250‑1 (MQ=255)
                                    tgCTAACAACGCGCCAAAACCCCGGAAGCGGACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGTTTGCCTGCCCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGAc                                                                                                                                                                                                       <  1:330860/250‑1 (MQ=255)
                                        accgacgCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAATCCTGAAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGAcgcg                                                                                                                                                                                                    <  1:260656/247‑1 (MQ=255)
                                        aaCGACGCGGCAAAACCCCGGAAGCGTACATCAGCCCGCGACCGGAGCGGGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGACCAGTAACACTTCGCTGTCGTTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGTTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACgcgg                                                                                                                                                                                                   <  2:164348/250‑1 (MQ=255)
                                        aaCGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGCGCGAGGAAGCCTGCACAAAGGCAGGCTTCAAACATACCGTTGAGTTTTACACCGGCGGCTGATCGAACCGTACCACTTCTCTGTCGCTATCCGCATGGATGGAGCTTCCCTGCTCATCCCAGATTGACAGACCATCGCTGGTCGCGACTTTCACGCCGTTACTGGAGACATTGCCTTTCAACACCTGGATCCATAcgcgac                                                                                                                                                                                                  >  1:424929/1‑249 (MQ=255)
                                        aaCGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACgcggc                                                                                                                                                                                                  <  1:387236/251‑1 (MQ=255)
                                        aaCGAAGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGCGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGCGCTTTAAACCGGCGGCAGATCGACCAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCCTCCCAGATTGCCAGACCCTCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACgcggc                                                                                                                                                                                                  >  2:237384/1‑251 (MQ=255)
                                            acgCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCTGCGTTc                                                                                                                                                                                              >  2:106565/1‑251 (MQ=255)
                                                                   aCATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGAc                                                                                                                                                                        <  1:260089/250‑1 (MQ=255)
                                                                         gCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATTCACCGACTGCTCAt                                                                                                                                                                 <  1:317098/251‑1 (MQ=255)
                                                                                  cGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTg                                                                                                                                                                                                                                    >  1:391115/1‑175 (MQ=255)
                                                                                  cGGAGCGAGCGAGGAAGCCTGAAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTg                                                                                                                                                                                                                                    <  2:391115/175‑1 (MQ=255)
                                                                                   ggAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCa                                                                                                                                                        >  1:282218/1‑250 (MQ=255)
                                                                                      gcgagcgaGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGc                                                                                                                                                    <  1:167615/251‑1 (MQ=255)
                                                                                       ggggcgagGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCAGGGATGGAGATTGCCTGCTCATCCCAGATGGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGcc                                                                                                                                                   <  2:98872/248‑1 (MQ=255)
                                                                                        gagcgagGAAGCCTGCAAAACGGCAGGCTTCAAAGACGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGCTGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGccc                                                                                                                                                  >  2:333250/1‑251 (MQ=255)
                                                                                                aaGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTAc                                                                                                                                          >  2:142166/1‑251 (MQ=255)
                                                                                                aaGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTAGAGACTTTCACGCCGTTAATGGGGACATTGCCTTTCACAACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTAc                                                                                                                                          >  1:338484/1‑251 (MQ=255)
                                                                                                   ccTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTAAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCAGCCCAGATTGCCAGACCAGCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGt                                                                                                                                       <  2:402174/251‑1 (MQ=255)
                                                                                                      gCAAAAAGGCAGCCTTCAAAGAGTCCGTTGAGTTTCAACCCGGGGGCAGAGCGAACAGTAACCCTTCGCTGTCGCTATCCGCCTGGATGGGGATTGCCTGCTCATCCCAGATTGCCAGTCCATCGCTGGTCGCGCCTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCACCGCCCAGCGGTACAGTTcc                                                                                                                                    <  2:1000/251‑1 (MQ=255)
                                                                                                       cAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTcc                                                                                                                                    >  1:417770/1‑250 (MQ=255)
                                                                                                            aggcaggcTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATcc                                                                                                                              >  1:58777/1‑251 (MQ=255)
                                                                                                               caggcTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTg                                                                                                                            >  2:14758/1‑250 (MQ=255)
                                                                                                                aggcTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGAt                                                                                                                          >  1:263212/1‑251 (MQ=255)
                                                                                                                  gcTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACg                                                                                                                                                     >  2:233330/1‑222 (MQ=255)
                                                                                                                  gcTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCGGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACg                                                                                                                                                     <  1:233330/222‑1 (MQ=255)
                                                                                                                     tCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACt                                                                                                                     >  1:125591/1‑251 (MQ=255)
                                                                                                                      cAAAGATGCCGTTGAGTTTTAACCCGGCGGCAGATCGAACAGTAACCCTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACt                                                                                                                     <  1:388717/250‑1 (MQ=255)
                                                                                                                      cAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACtt                                                                                                                    >  2:252230/1‑251 (MQ=255)
                                                                                                                                gTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCg                                                                                                                                                                          >  2:113646/1‑187 (MQ=255)
                                                                                                                                gTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCg                                                                                                                                                                          <  1:113646/187‑1 (MQ=255)
                                                                                                                                   gAGTTTTAAACCGGCGGCAGATCGCACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCAtc                                                                                                       >  2:95508/1‑251 (MQ=255)
                                                                                                                                         tAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGCCTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACAATCTCgcgc                                                                                                 >  2:123301/1‑251 (MQ=255)
                                                                                                                                            aCCGGAGGCAGATCGAACAGTAACAATTCGATGTCGCTATCCGCATGGATGGCGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAACAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATc                                                                                              >  2:220585/1‑251 (MQ=255)
                                                                                                                                             ccGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATcc                                                                                             >  1:15380/1‑251 (MQ=255)
                                                                                                                                             ccGGCGGCAGATCGAACAGTAAAACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCCCCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGAACAGTTCCATATCATGATGCACTTTCAATGAACCATCTCGCGCATcc                                                                                             >  2:87095/1‑251 (MQ=255)
                                                                                                                                                     aGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCgag                                                                                       >  2:458457/1‑249 (MQ=255)
                                                                                                                                                             cAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCACGCGCATCCGGCGAGAGCACCAGCt                                                                             >  2:303387/1‑251 (MQ=255)
                                                                                                                                                                tAACACTTCGCTGTAGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGCTCCAGACGCGGCGTTCAGCGGAAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGtt                                                                          >  1:148311/1‑251 (MQ=255)
                                                                                                                                                                    acTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGcc                                                                      >  1:386596/1‑251 (MQ=255)
                                                                                                                                                                       tCGCTGTCGCTATCCGAATGGATGGAGATTGCCTGCTCATCCCAGCTTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTCCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCAACGACTGCTCATCTTTCAGCAACGCCCCGCGGTACAGTTACATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACAAGCTGTTTGCCCTg                                                                   >  2:266679/1‑251 (MQ=255)
                                                                                                                                                                           tGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACg                                                               >  1:431144/1‑251 (MQ=255)
                                                                                                                                                                            gTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACgg                                                              >  1:29405/1‑251 (MQ=255)
                                                                                                                                                                            gTCGCTATCCGCAGGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTAc                                                                <  1:94271/249‑1 (MQ=255)
                                                                                                                                                                              cGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCa                                                            <  1:248096/251‑1 (MQ=255)
                                                                                                                                                                                 tATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCg                                                         >  1:3285/1‑251 (MQ=255)
                                                                                                                                                                                    ccGCATGGATGGAGATTGCCTGCTCATCCCAAAGTGCTATACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCATCTGTAGCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGACCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCCCCAGCTGTTTGCCCTGTACGGCATCGAAg                                                      <  1:90951/251‑1 (MQ=255)
                                                                                                                                                                                          ggatggaGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATAGCCTGATGCTCTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACg                                                 <  1:300188/250‑1 (MQ=255)
                                                                                                                                                                                                     gCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAAc                                      <  1:252230/250‑1 (MQ=255)
                                                                                                                                                                                                      ccTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAAcgc                                    >  1:39752/1‑251 (MQ=255)
                                                                                                                                                                                                               cccAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTact                            >  1:342176/1‑248 (MQ=255)
                                                                                                                                                                                                                      ggCCAGACCATCGCGGGCCGAGACTGTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGACTCAGACGCGGCGTGCAGCGGCAATCTGATGCACCGCCTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTACCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGttt                     <  2:128531/249‑1 (MQ=255)
                                                                                                                                                                                                                      ggCCAGACCATCGCCGGTCGAGCCTTTCACGCCGTTAATGGTGACATTGCCTTTCACCCCCTGGATCCTGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTGCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCGGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTAATTCCGttt                     <  2:128528/249‑1 (MQ=255)
                                                                                                                                                                                                                        ccAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGTTttct                  >  2:451398/1‑251 (MQ=255)
                                                                                                                                                                                                                          ggACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGTTttcttc                <  2:349970/250‑1 (MQ=255)
                                                                                                                                                                                                                            aCCATCGCTGGCCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGTTTTCTTCAg              <  1:106565/251‑1 (MQ=255)
                                                                                                                                                                                                                               aTCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGTTTTCTTCAGGc            >  1:407293/1‑250 (MQ=255)
                                                                                                                                                                                                                                     tgTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGCGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTCATCCCGTTTTCTTCAGGCATGATcc     >  2:36160/2‑251 (MQ=255)
                                                                                                                                                                                                                                     ggTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGTTTTCTTCAGGCATGATcc     <  1:122956/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCAACTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGTTTTCTTCAGGCATGATcc     >  1:386915/1‑250 (MQ=255)
                                                                                                                                                                                                                                       tCGAGACTTTCACGCCGCTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGTTTTCTTCAGGCATGCTcc     >  1:151149/1‑249 (MQ=255)
                                                                                                                                                                                                                                        cGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGGTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGTTTTCTTCAGGCATGATCCAGa  <  2:386596/251‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
TTATTAGATATACACTAAATCAATTGCATCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGGCTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGTTTTCTTCAGGCATGATCCAGA  >  NC_000913/3578814‑3579296

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