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.87.2% 155.3 / 17.8 47coding (140/1161 nt)metBO‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base . (22/19);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01
*NC_0009134,128,8120G.87.2% 156.7 / 15.6 47coding (141/1161 nt)metBO‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase
Reads supporting (aligned to +/- strand):  ref base G (2/4);  new base . (22/19);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAAC  >  NC_000913/4128577‑4129053
                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                                 
agaagTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTACCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAAc                                                                                                                                                                                                                                  <  1:338601/251‑1 (MQ=255)
          aTGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGAtgtg                                                                                                                                                                                                                        <  1:369477/251‑1 (MQ=255)
           tGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCCACGCGCGATGTgg                                                                                                                                                                                                                       >  1:162317/1‑251 (MQ=255)
                  aTGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCg                                                                                                                                                                                                                >  1:118975/1‑251 (MQ=255)
                    gTATTGAAGACCATAAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCCTCACATGACGCGTAACCAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGCTTTAATGAACCGCGCGCGCATGTTTACT‑‑CGTCGCGGCACCCCAACGCGCGATGAGGTTCAGCggg                                                                                                                                                                                                              >  2:467131/1‑249 (MQ=255)
                           cGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGc                                                                                                                                                                                                       >  1:488265/1‑251 (MQ=255)
                                      ccaATGTTCTCTCTGGTGCCTGACATTTCACCGACAAAGCCCATTGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCTTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTTCCACCGATCCATCATTCCAGCACCTACAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCTTGCGCTGGCAGAACTGGAAg                                                                                                                                                                                            <  2:7664/250‑1 (MQ=255)
                                                   tGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctggggc                                                                                                                                                                               >  1:321270/1‑248 (MQ=255)
                                                      tGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGt                                                                                                                                                                            >  1:404737/1‑251 (MQ=255)
                                                                  cACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACATATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCACCGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATAc                                                                                                                                                                >  1:175954/1‑251 (MQ=255)
                                                                     cGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCgg                                                                                                                                                             >  2:189585/1‑251 (MQ=255)
                                                                        cAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCAt                                                                                                                                                          <  2:120300/251‑1 (MQ=255)
                                                                          aaGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGCTGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTTGTGCTGGTGCAGTACTTACTAATACCGGCATGt                                                                                                                                                        >  1:187632/1‑251 (MQ=255)
                                                                               cAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTAGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCcgcg                                                                                                                                                   >  2:62498/1‑251 (MQ=255)
                                                                               cAGGGAACTCCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACGTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCcgcg                                                                                                                                                   <  1:465916/251‑1 (MQ=255)
                                                                                   gAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTTCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGACGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTc                                                                                                                                               >  1:322395/1‑251 (MQ=255)
                                                                                   gAACTTCATCACATGACGCGTAAACAGGCCACCAGCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTGTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACGAATACCGGCATGTCCGCGATTc                                                                                                                                               <  2:533410/251‑1 (MQ=255)
                                                                                          atcaCATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGt                                                                                                                                        >  2:127593/1‑251 (MQ=255)
                                                                                                         aaCAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGGTCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCttttt                                                                                                                         >  1:507797/1‑251 (MQ=255)
                                                                                                                         aGTGCGTAGCGGGTTAACTGACGACGAACAGTATGGGTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCAGGGAACGACAGTCTTTTTTAAACCTGGCGATctg                                                                                                         >  2:299796/1‑251 (MQ=255)
                                                                                                                                aGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTg                                                                                                  <  2:537045/251‑1 (MQ=255)
                                                                                                                                 gCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTgc                                                                                                 >  1:100620/1‑251 (MQ=255)
                                                                                                                                  cGGGTTAAATGACGCCGAACAGTATGGGTGCGTGGTCCCACCGATCCAGCTGCCCAGCACCTATAACTTTACCGGATTTAATGACCCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCACCGCGCGATGTGGTTCAGCGTGCGCTGTCAGAACTGGAAGGTGGTGCTGGTGCAGTACTAACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTGTGAATCCTGGCTATCTGCTGGTTgcg                                                                                                <  2:88659/251‑1 (MQ=255)
                                                                                                                                            gacgacGAACAGTATGGTTGCGTTGTCCAACCGATCCATCTTTCCAGCACCTATAACTTTCCCGGATTTAATGAACCGCGCGCGCATGATTCCTCGCGTCGCGGCACCCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCGGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGTAACCTGGCGATCTTCTGGTTGCGCCGCACgc                                                                                        >  2:42371/1‑250 (MQ=255)
                                                                                                                                               gacgaACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTg                                                                                     >  1:52731/1‑251 (MQ=255)
                                                                                                                                                  gaACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTAcg                                                                                <  1:129176/251‑1 (MQ=255)
                                                                                                                                                  aaaCATTATGGTTTCGTTGTCCCCCCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGAGTACT‑‑CGTCGCGGCATCCCAACGCGCGATGTGGTTTAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTAct                                                                                <  2:289696/250‑2 (MQ=255)
                                                                                                                                                   aaCAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCTCCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTAcgg                                                                               <  2:163894/251‑1 (MQ=255)
                                                                                                                                                               gggTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACAGCGCGCGCATGATTACT‑‑CGTAGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGACACCTGGCGATCTGCTGGTTGCGCCGCACGACTGATCcggcgggtgtcatcc                                                                   >  2:195716/3‑242 (MQ=255)
                                                                                                                                                                ggtggacccAACGATCCATCTTTCCAGCACCTATAACTGTACCGGATGTAATGAACCGCGCGCGCATGATTACTCG‑G‑CGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGTAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGc                                                                  <  2:238888/245‑1 (MQ=255)
                                                                                                                                                                cGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGATACCTGTCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGGAGCTATccc                                                                  >  2:258171/1‑249 (MQ=255)
                                                                                                                                                                     tCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGAGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGCTTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTAcggcgggtgactttccctgct                                                             >  2:21564/1‑236 (MQ=255)
                                                                                                                                                                        cACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCCATCGCCCGTTCgg                                                          >  1:271798/1‑250 (MQ=255)
                                                                                                                                                                             aTCCATCTTTCCAGCACCTATAACGTTACCGGATTTAATGAACCGCGCGCGCAGGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTGCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCGGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTc                                                     <  2:100620/251‑1 (MQ=255)
                                                                                                                                                                                      cccAGCACCTATACATGTACCGGATATAATGAACCGCGCGCGCATGATTACTCG‑G‑CGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACGGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGGCTTTTTGAATCCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGGCTGGCGAAAc                                            <  2:162317/250‑1 (MQ=255)
                                                                                                                                                                                         aGCACCTATAACTTTACCGGATTTAATGAACCGCGCGTGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCGGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGCCTTTTTGAAACCTGGCGATCTGCTGGTGGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAAcg                                           <  1:42371/251‑1 (MQ=255)
                                                                                                                                                                                           cACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGt                                       <  1:387655/251‑1 (MQ=255)
                                                                                                                                                                                             ccTATAACCTTACCGGAATTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGAGGCAACCCAACGCGTGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAAGACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGTAACGCGGTTg                                     <  2:1206/251‑1 (MQ=255)
                                                                                                                                                                                                 taACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAACTATCACCTGTTCGCCAGTATGGCGAAACGCGGTTGCTAt                                 >  1:37441/1‑251 (MQ=255)
                                                                                                                                                                                                         ccGGATTTAATGAACCGCGCGCGCATGATTACT‑‑TGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATAGCCTGTTCGACAGTCTGGGGAAACGCGGGTGCTATcgccggtt                         >  1:359167/1‑246 (MQ=255)
                                                                                                                                                                                                                       ccGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGtt                                                             <  1:508083/201‑1 (MQ=255)
                                                                                                                                                                                                                         gcgcgcgcATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGc         <  1:127593/251‑1 (MQ=255)
                                                                                                                                                                                                                             gcgcAGGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGGCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATg     <  2:321270/251‑1 (MQ=255)
                                                                                                                                                                                                                              cgcATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGGACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGGCTGCTACGGCGGTAGCTATCGGCTGTTCGACAGGCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTGGATCAAGGCGag      <  2:52731/251‑2 (MQ=255)
                                                                                                                                                                                                                               gcATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTGAGCGTGCGCTTGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGaa   <  1:120377/251‑1 (MQ=255)
                                                                                                                                                                                                                                cATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAAc  <  1:68854/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aTGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATTGTCTGTTCGACCGACTTGCGCAAAGCGGGTGCTTTCTCGGGTTGTTTGTTTGTCATGGCGATGaca  >  1:285710/1‑248 (MQ=255)
                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                                 
AGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAAC  >  NC_000913/4128577‑4129053

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: ‑