Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,128,813 Δ2 bp coding (142‑143/1161 nt) metB → O‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,128,8110C.98.4% 289.5 / ‑1.7 64coding (140/1161 nt)metBO‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (36/27);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01
*NC_0009134,128,8120G.98.4% 289.7 / ‑3.6 64coding (141/1161 nt)metBO‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (36/27);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACA  >  NC_000913/4128561‑4129054
                                                                                                                                                                                                                                                           ||                                                                                                                                                                                                                                                  
ctccaTCTATACGCAAAGCAGTTTCGATGCCCAGCTGGATTGACGTCCATTAACACAATGCTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACCCGTAAACCGGCCACCCTCGCAGTGCGTAGCGGGTTAACTGACGACGCACAGTATGGTTGCGTTGTACCACCGATCCATCTTTGCACCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACtc                                                                                                                                                                                                                                                     >  1:50282/5‑251 (MQ=255)
  cAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACtcgt                                                                                                                                                                                                                                                   >  2:147227/1‑250 (MQ=255)
       tataCGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACtcgtcgcgg                                                                                                                                                                                                                                              >  2:336229/1‑245 (MQ=255)
        ataCGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACtcgtcg                                                                                                                                                                                                                                                 <  1:185671/247‑4 (MQ=255)
                agaagTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAAc                                                                                                                                                                                                                                   <  1:170416/251‑1 (MQ=255)
                 gaagTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAAcg                                                                                                                                                                                                                                  >  1:44277/1‑251 (MQ=255)
                  aagTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAACCCCAGGAAACTTCATCACACGACGAGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAAAGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATATTTCCAGCACC‑TATAACTTTACAGGATTTAATGAACCGCGCGCGCCTGATTACT‑‑CGTCGCGGCAACCCAAcgc                                                                                                                                                                                                                                 >  1:312314/1‑251 (MQ=255)
                        aGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGAtg                                                                                                                                                                                                                           >  1:391626/1‑251 (MQ=255)
                        aGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGAtg                                                                                                                                                                                                                           >  2:414336/1‑251 (MQ=255)
                                aGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTCTGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTc                                                                                                                                                                                                                     >  2:337391/1‑249 (MQ=255)
                                        tGACGTACATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTAGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCt                                                                                                                                                                                                           >  1:440122/1‑251 (MQ=255)
                                            gTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGc                                                                                                                                                                                                        >  1:22630/1‑250 (MQ=255)
                                                 ttAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAACCAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCCGCGTGCGCTGGCAGAACt                                                                                                                                                                                                  >  1:439109/1‑251 (MQ=255)
                                                  tAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTg                                                                                                                                                                                                 >  2:235830/1‑251 (MQ=255)
                                                             ttACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGTGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGgtggtg                                                                                                                                                                                       <  1:185940/250‑1 (MQ=255)
                                                               aCTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctg                                                                                                                                                                                    <  2:137442/251‑1 (MQ=255)
                                                                ctctGGTGCCTGACCTTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctgg                                                                                                                                                                                   >  1:368891/1‑251 (MQ=255)
                                                                   tGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctggtgc                                                                                                                                                                                >  2:350822/1‑251 (MQ=255)
                                                                    ggTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAATGtggtgctggt                                                                                                                                                                                  <  1:357219/248‑1 (MQ=255)
                                                                    ggTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCa                                                                                                                                                                               <  1:187792/251‑1 (MQ=255)
                                                                     gTGCCTGACAGTTCACCGACAAAGCCCAGGGAACGTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAg                                                                                                                                                                              <  1:342316/251‑1 (MQ=255)
                                                                       gCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGGGCAGt                                                                                                                                                                             >  2:90008/1‑250 (MQ=255)
                                                                            aCATTTCCCAGACAAAGACCAGGGAACTTCATCACATGACGCGTAAACAGGCCCCCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAAGTGGTGCTGGTGCAGtacttac                                                                                                                                                                       >  2:153047/1‑251 (MQ=255)
                                                                            aCATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGtacttac                                                                                                                                                                       <  1:235830/251‑1 (MQ=255)
                                                                              aTTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGtacttact                                                                                                                                                                      >  1:403287/1‑250 (MQ=255)
                                                                                  cACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATa                                                                                                                                                                  <  2:120512/250‑1 (MQ=255)
                                                                                     cGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGACCAGTATGGATGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGACGCGGCAACCCAACGCGCGATGAGGTTCAGCGTGCTCTGGCAGAACTGGAAGGTGGTGCTGGAGCAGAACTTACTAATACCgg                                                                                                                                                              >  2:40567/1‑251 (MQ=255)
                                                                                     cGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGTCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGACCTGGAAGGTGGTGCTGGTGCAGTCCTTCCTAATACCgg                                                                                                                                                              >  2:405/1‑251 (MQ=255)
                                                                                      gACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGc                                                                                                                                                             >  1:6117/1‑251 (MQ=255)
                                                                                         aaaGTCCAGGGAACTTCATCCCATGACGCGTAAACAGGCCACCATCGAAGTGCGTAGCGGGTTAAATGACGACGAACAGGGTGGTTGCGCTGTCCCACCGCTCAATCTTTCCAG‑ACCGTATAACTTTAACGGACTTAATGAACCGCGAGCGCATGATTACT‑‑AGTCGCGGAAACCCAAAGCGCGCTGTGGTTCAGCGTACGCTGGCAGAAATGGAAGGTGGTGCTGGTGCAGAACATACTCATACAGGcctg                                                                                                                                                          >  1:393045/1‑248 (MQ=255)
                                                                                           atcccATGGAACTTCATCGCATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAAGGACGACGAACAGTATTGTCGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTGTACCGGATCTAATGAACCGCGCGCTCAGGATTACT‑‑CGTCGCGGCATCCCATCGCGCGAGGTGGTTCAGCGTGCGCTGGCAGAACCGGAAGGTGGTGCTGGTGCAGTGTTTCCTAATACCTGCATGTc                                                                                                                                                        <  2:50282/249‑1 (MQ=255)
                                                                                           aGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTc                                                                                                                                                        <  2:402936/251‑1 (MQ=255)
                                                                                           aGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTc                                                                                                                                                        >  1:367091/1‑251 (MQ=255)
                                                                                               cAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCATCCCAACGCGCGATGTGGTTCAGCGGGCGCTGGCATAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGTATGGCcgcg                                                                                                                                                    <  2:443050/251‑1 (MQ=255)
                                                                                                  ggAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGAtt                                                                                                                                                 <  2:53756/251‑1 (MQ=255)
                                                                                                   gAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAAGACCGGCATGTCCGCGATTc                                                                                                                                                <  2:44277/251‑1 (MQ=255)
                                                                                                       ttCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGAATGTCCGCGATTCACCt                                                                                                                                            >  2:366672/1‑251 (MQ=255)
                                                                                                        tcatcaCATGACGCGTAAACAAGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCCCCGATCCATCTTTCCAGCACC‑TATAACTGTACCGGATTTAATGAACCGCGCGCGCATGATGACT‑‑CGTCGCGGTAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCt                                                                                                                                            <  2:105021/250‑1 (MQ=255)
                                                                                                              caTGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGAGGTTCAACTGGTAACg                                                                                                                                     >  2:140566/1‑251 (MQ=255)
                                                                                                                 gACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGAc                                                                                                                                   <  1:414336/250‑1 (MQ=255)
                                                                                                                                                    ggTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGcgcc                                                                                               >  2:172693/1‑251 (MQ=255)
                                                                                                                                                     gTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGcgccg                                                                                              >  1:278736/1‑251 (MQ=255)
                                                                                                                                                      ttAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGcgccgc                                                                                             >  1:308796/1‑251 (MQ=255)
                                                                                                                                                               gacgaACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCt                                                                                    <  1:147227/251‑1 (MQ=255)
                                                                                                                                                                 cgaACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTAc                                                                                  <  1:177589/251‑1 (MQ=255)
                                                                                                                                                                      aGTATGGTTGCGTTGGCCCACCGATGCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGGATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTAcggcgg                                                                             <  2:391626/251‑1 (MQ=255)
                                                                                                                                                                        tATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGt                                                                            >  1:356409/1‑250 (MQ=255)
                                                                                                                                                                         aTGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTa                                                                           <  2:368891/250‑1 (MQ=255)
                                                                                                                                                                              tGCGTTGTCCCACCGATCCATCTTTTCAGCACC‑TATAACTTTACCGGGTTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGCAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACCGATAcggcgggaactca                                                                      >  2:431273/1‑243 (MQ=255)
                                                                                                                                                                              tGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGCTGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGGAGCTATc                                                                     >  1:196794/1‑251 (MQ=255)
                                                                                                                                                                                cGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGc                                                                   >  1:286459/1‑251 (MQ=255)
                                                                                                                                                                                   tGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGccattccccg                                                                >  2:232378/1‑242 (MQ=255)
                                                                                                                                                                                     tCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCAAGACTGCTAcggcggcagatctcgcccgtt                                                              >  1:459137/1‑236 (MQ=255)
                                                                                                                                                                                      ccaaCCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGtt                                                              <  2:403287/247‑1 (MQ=255)
                                                                                                                                                                                        cACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGC‑‑CGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCAGCGATTCACCTGGTAACGACCGTCTTTTTGAAAACTGGCGATCTGCTGGTTGCGCCGCACGACTGCTAcggcggcaccacaaccctatttac                                                           >  2:88740/1‑233 (MQ=255)
                                                                                                                                                                                        cACCGATCCATATTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCAAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTAACCTGGTAAAGCCCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTAcggcgggagacatcgaccgtttg                                                            >  2:419958/1‑233 (MQ=255)
                                                                                                                                                                                            gATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAAACTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGAAGCTATCGCCTGTTCGAccgg                                                       >  2:88480/1‑248 (MQ=255)
                                                                                                                                                                                            gATCCATCTTTCCAGCACC‑TATAACGTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGt                                                       <  2:293691/251‑1 (MQ=255)
                                                                                                                                                                                                       ccAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAAcg                                            <  1:350822/251‑1 (MQ=255)
                                                                                                                                                                                                            aCC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGtt                                       <  1:419267/251‑1 (MQ=255)
                                                                                                                                                                                                             cc‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGGTCGGCAAGCTGGCGGAACACGGTTg                                      >  1:122189/1‑251 (MQ=255)
                                                                                                                                                                                                                 ataACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCCATCGCCTGTTTGGCCGTCCGGCGAAACGCGGTGGCtt                                   >  1:443741/1‑250 (MQ=255)
                                                                                                                                                                                                                   aaCTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTTGACAGTCTGGGGAAACGCGGTTGCTATc                                 >  1:176709/1‑251 (MQ=255)
                                                                                                                                                                                                                         aCCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCgtgt                           <  2:356409/251‑1 (MQ=255)
                                                                                                                                                                                                                          ccGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGCTCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTTGCGAAACGCGGTTGCTATcgcggggt                          >  1:51843/1‑247 (MQ=255)
                                                                                                                                                                                                                                  aaTGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTtgtttgtt                   <  1:405/250‑1 (MQ=255)
                                                                                                                                                                                                                                          gcgcgcgcATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGc          <  1:230336/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  aTGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACGTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACa  <  2:352707/251‑1 (MQ=255)
                                                                                                                                                                                                                                                           ||                                                                                                                                                                                                                                                  
CTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACC‑TATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACA  >  NC_000913/4128561‑4129054

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