Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,129,874 C→T L14L (CTG→TTG)  metL → fused aspartate kinase/homoserine dehydrogenase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,129,8740CT97.1% 211.6 / ‑1.7 69L14L (CTG→TTG) metLfused aspartate kinase/homoserine dehydrogenase 2
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (29/38);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

TTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCC  >  NC_000913/4129633‑4130119
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
ttACGATGGCGGAATCATTAGGGGGAGTGGAAAGTTTAATCTCCCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAAt                                                                                                                                                                                                                                              >  1:212428/1‑251 (MQ=255)
      tGGCGGAATCATTAGGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCAGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGCTTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTtggtg                                                                                                                                                                                                                                        >  2:60250/1‑251 (MQ=255)
                    gggggAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGAATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCt                                                                                                                                                                                                                          <  1:605177/251‑1 (MQ=255)
                         ggTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTTGTGGCATTAGTCTGGCTGAtgt                                                                                                                                                                                                                     <  2:643928/250‑1 (MQ=255)
                             gAAAGTTTAATCTCTCACGCCGCAACCATGACACAGGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAg                                                                                                                                                                                                                 <  1:280600/251‑1 (MQ=255)
                               aaGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACATGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAgtg                                                                                                                                                                                                               <  1:326116/251‑1 (MQ=255)
                                   ttAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCTAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTCGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTa                                                                                                                                                                                                            <  1:612503/250‑1 (MQ=255)
                                    tAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGCATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTAtt                                                                                                                                                                                                          >  1:68951/1‑251 (MQ=255)
                                      acctcgcACGCCGCAACGATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCGGGGATCTCCGAGACGCTGCTGCGTAGCTACACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTGAAAATGAGTGTGAGTGCGCAGGCAGGGGCGAAAGGTAGTCAGTTGCATAAATGTGGTGGCAGTAGTCTGGCGGATGGGAATTGTTATTTg                                                                                                                                                                                                        <  2:621896/245‑1 (MQ=255)
                                        ctctcACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCg                                                                                                                                                                                                      >  1:20352/1‑251 (MQ=255)
                                          ctcACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTg                                                                                                                                                                                                        <  2:338896/247‑1 (MQ=255)
                                               gccgcAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTcgcg                                                                                                                                                                                               <  1:243141/251‑1 (MQ=255)
                                                  gcAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGc                                                                                                                                                                                            >  2:245605/1‑251 (MQ=255)
                                                     aCCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCAtt                                                                                                                                                                                         <  2:546341/251‑1 (MQ=255)
                                                               aTGCAGGCATGGCACCAGCAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGCTGCATAAATTtggtgg                                                                                                                                                                                                                                       >  1:184966/1‑195 (MQ=255)
                                                               aTGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGCTGCATAAATTtggtgg                                                                                                                                                                                                                                       <  2:184966/195‑1 (MQ=255)
                                                                 gCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTa                                                                                                                                                                              <  1:55365/250‑1 (MQ=255)
                                                                   aGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGCAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTActc                                                                                                                                                                           >  1:36632/1‑251 (MQ=255)
                                                                            aCCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTg                                                                                                                                                                   >  1:567236/1‑250 (MQ=255)
                                                                             ccAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGTAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGAc                                                                                                                                                                 <  1:249512/251‑1 (MQ=255)
                                                                              cAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCTGGCATTATGGCGGAGTACTCTCAGCCTggcg                                                                                                                                                                >  2:357019/1‑248 (MQ=255)
                                                                                     gcgTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATatgat                                                                                                                                                         >  2:136626/1‑251 (MQ=255)
                                                                                      cgTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTa                                                                                                                                                                              <  2:220962/229‑1 (MQ=255)
                                                                                      cgTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTa                                                                                                                                                                              >  1:220962/1‑229 (MQ=255)
                                                                                        tgctgcCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGAtggg                                                                                                                                                      >  1:591137/1‑250 (MQ=255)
                                                                                         gctgcCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGAtggtg                                                                                                                                                     <  2:543461/251‑1 (MQ=255)
                                                                                            gccgccGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGGCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGtt                                                                                                                                                  <  2:577116/251‑1 (MQ=255)
                                                                                              cgccgGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCGGGAAAATGGCTTCCGGGCTGCAAAAAAGGGGTAAAACTGAGTGTGATTGCGCAGGCAGGGGCGAAAGTTCGTCAGTTGCCTAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCCGACGGCATGATGGTGGTTTc                                                                                                                                                >  2:642469/1‑251 (MQ=255)
                                                                                               gccgGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTcc                                                                                                                                               <  1:60250/251‑1 (MQ=255)
                                                                                                     atcgccGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTccgtcgcc                                                                                                                                         <  2:12104/247‑1 (MQ=255)
                                                                                                         ccGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTccgccgccg                                                                                                                                        <  2:313647/248‑1 (MQ=255)
                                                                                                         ccGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTccgccgcc                                                                                                                                         <  1:136020/247‑1 (MQ=255)
                                                                                                          cGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGt                                                                                                                                      >  2:462734/1‑249 (MQ=255)
                                                                                                           gagaCGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGc                                                                                                                                   >  1:587615/1‑251 (MQ=255)
                                                                                                           gagaCGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGc                                                                                                                                   <  1:363962/251‑1 (MQ=255)
                                                                                                             gaCGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGcac                                                                                                                                 <  1:202252/251‑1 (MQ=255)
                                                                                                               cGCTACTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGcacc                                                                                                                                >  1:640680/1‑250 (MQ=255)
                                                                                                                  tgctgcGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGCAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGCCGTCAGTTGCATAAATTTGTCGGCAGTAGTCTGGCTGATGTGACGTGTTATTTGCGTGTCGCGGGCCTTATGGCGGCGTACTTTCAGCCTGACGACATGCTGGCGGTTTCCGCCGCCGGGAGCACCACt                                                                                                                             >  2:307379/1‑250 (MQ=255)
                                                                                                                   gctgcGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGAGATGATGGTGGTTTCCGCCGCCGGTAGCACCACTaa                                                                                                                           <  2:311550/251‑1 (MQ=255)
                                                                                                                            gtcCACCGGGATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCATTTGCATAAATTGGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTCCTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGAt                                                                                                                  <  2:19689/250‑1 (MQ=255)
                                                                                                                                aCCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGATGCAAACAAGGGGTAAAAATGAGAGTGATTGCGCAGGCAGGGGCGAAAGGTCGCCAGTTGCATAAATTTGGCGGCAGTAGTCTGGATGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACAACTAACCAGTTGAttcac                                                                                                              >  2:311025/1‑248 (MQ=255)
                                                                                                                                  gggTATTGAAGATGGGGAAGATTTAATTGCCGACCGGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGTAGTACTCTCAGCCTGACGATATGATGGGTGTTTCCGCCGCCGGTAGCACCTCTAACCAGTTGATTAACTg                                                                                                            <  2:84165/250‑1 (MQ=255)
                                                                                                                                           aGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACt                                                                                                   <  1:195284/251‑1 (MQ=255)
                                                                                                                                             aTGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTTCTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTaa                                                                                                 <  1:312439/251‑1 (MQ=255)
                                                                                                                                                  gAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGAAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTGATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAg                                                                                            <  2:178003/251‑1 (MQ=255)
                                                                                                                                                    aGATTTAATTTCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTTCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCTTGTCGCGGGCATTATGGCGGAGTCCTCTCAGCCTGACGATATCATGGTGGTTTCCTTCGCCGGTAGCACCACTAACCa                                                                                                                        <  2:75872/221‑1 (MQ=255)
                                                                                                                                                      aTTTAAGTGCCGACCTGGAAAATGCCTTCCGGGCGTCAAACAAGGGGAAAAAATGAGTGTGAATGAGCAGGCAGGGGCGAAAGGTCGTCAGTAGGATAAATTTGGTGGCAGTGGTCTGGCTGATGTGAATTGTTATTTGCGTTTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGACGGTGGTTTCCGCCGCAGGTAGCACCCCTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGaacg                                                                                        <  2:44771/251‑4 (MQ=255)
                                                                                                                                                         tAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGAt                                                                                      <  1:136626/250‑1 (MQ=255)
                                                                                                                                                               ccGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGACGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATAGtctcc                                                                               >  2:112467/1‑250 (MQ=255)
                                                                                                                                                                       gAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGAGTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCAt                                                                       <  2:340508/251‑1 (MQ=255)
                                                                                                                                                                         aaaTGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCa                                                                     <  1:165533/251‑1 (MQ=255)
                                                                                                                                                                                 tCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTca                                                               <  1:97909/249‑1 (MQ=255)
                                                                                                                                                                                  ccGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTcaac                                                             >  2:55027/1‑250 (MQ=255)
                                                                                                                                                                                   cGGGCTGCCAACAAGGGGTAACAATGAGAGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTACCAGCCTGACGATCTGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAACCTAAGCCAGACCGAACGTATCTCTGCGCATCCGGTTcaaca                                                            >  2:115138/1‑250 (MQ=255)
                                                                                                                                                                                      gCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGTGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTATCACCACTAACCAGTTGATTAACTGGTTGAAACTCATCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACg                                                        >  1:616242/1‑251 (MQ=255)
                                                                                                                                                                                         gCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTg                                                     >  1:402880/1‑251 (MQ=255)
                                                                                                                                                                                         gCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCGGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCt                                                      <  2:383271/250‑1 (MQ=255)
                                                                                                                                                                                           aaaNAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGGGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGAGCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCg                                                   <  2:591137/251‑1 (MQ=255)
                                                                                                                                                                                              cAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCg                                                >  2:463854/1‑251 (MQ=255)
                                                                                                                                                                                                    gTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATc                                           >  2:513985/1‑250 (MQ=255)
                                                                                                                                                                                                       aaaaTGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCGATAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTACAACAAACGCTGCGTCGCTATCAGTg                                       <  2:487715/251‑1 (MQ=255)
                                                                                                                                                                                                        aaaTGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGc                                      >  1:246076/1‑251 (MQ=255)
                                                                                                                                                                                                            gAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATc                                  >  1:510132/1‑251 (MQ=255)
                                                                                                                                                                                                             ggtgtGATTGCGCAGGCAGGGGCGAAAGGCCGTCTGTTGCATAAATTTGGTAGCAATAGTCTGGTGGAGGTGAAGTGTTATTTGCGTGTCGCGGGCTTTATGACGGAGTACGCTCAGCCTGACGATAGGATGGGGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCAGCAGGTTCAACAACCGCTGCGTCCCTAACAGTGCGATCt                                 <  2:245340/250‑1 (MQ=255)
                                                                                                                                                                                                                              aggGGCGAAAGGTCGTCAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGATAc                >  2:383266/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gtcgtcAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTgaaga      >  2:78012/1‑250 (MQ=255)
                                                                                                                                                                                                                                           cgtcAGTTGTATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGAGCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGc   <  2:80123/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gtcAGTTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGcc  <  1:78012/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gtcAGTTGAATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTACCTGGTTGAAACTAAGCCAGACCGACCGTCTCTCTGCGCATCAGGTTAAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTTGCGGTCTGCAACCCGCTGCAGAAGcc  >  2:335171/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
TTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGTATCTCCACCGGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCC  >  NC_000913/4129633‑4130119

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