Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,152,494 G→T A87S (GCG→TCG)  metC → cystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,152,4940GT94.4% 105.2 / NA 38A87S (GCG→TCG) metCcystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (24/10);  total (24/14)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGG‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGAC  >  NC_000913/3152247‑3152730
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              
aaaGCATGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGgt‑cgg                                                                                                                                                                                                                                            >  1:335053/1‑251 (MQ=255)
    cTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGc                                                                                                                                                                                                                                        >  1:297943/1‑251 (MQ=255)
        aTACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGtt                                                                                                                                                                                                                                    <  2:293128/251‑1 (MQ=255)
                                gcagcaAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGAGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATc                                                                                                                                                                                                           >  2:194313/1‑251 (MQ=255)
                                       aaaaTACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAgg                                                                                                                                                                                                    >  2:272274/1‑251 (MQ=255)
                                         aaTACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAAGCGCGCCAATGGAGAGTTGTTCGATGGACGGCGCGGAACGTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAAGTCC‑ATTCTTGCTTTTATCGAACAGGGc                                                                                                                                                                                                  <  2:245980/251‑1 (MQ=255)
                                             cacTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTA‑ACCCATTTCTCCGTACAACAAGCGATGTGTGAA‑CTGTAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATc                                                                                                                                                                                              <  1:241880/251‑1 (MQ=255)
                                               ctctCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGACGCCAAAAAACACGCGACACGTAATCGCGCCAATGGCGAGTTGTTCTATGGACGGCGCGGAACGTTA‑ACCCATTTCTCCTTACAACAAGCGCTGTGTGAA‑CTGGCAGGAGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCAt                                                                                                                                                                                            >  2:209793/1‑251 (MQ=255)
                                                               tAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACCGGGCGATCCTGTGTTGATGACCAac                                                                                                                                                                             >  2:27980/1‑250 (MQ=255)
                                                                   gTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGc                                                                                                                                                                        >  1:69325/1‑251 (MQ=255)
                                                                                      cGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCCATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGCGTCAGGAt                                                                                                                                                     >  2:353233/1‑251 (MQ=255)
                                                                                         tGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTc                                                                                                                                                  >  1:204637/1‑251 (MQ=255)
                                                                                                 gACAGTGTATCAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTAGTTCTATGGACGACGCGGAACGTTA‑ACCCATTTCTCCTTACAAGAAGCGATGTGTGCA‑CTGGAGGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAGTTCC‑ATTCTTGCTTTTATCGAAAGGGGCGCTCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGCAGTTCTGTAGCa                                                                                                                                           >  2:290478/1‑250 (MQ=255)
                                                                                                      tgtAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATcc                                                                                                                                     >  2:99199/1‑251 (MQ=255)
                                                                                                       gtAGAAGCAAAAAAACCCGGGACACGTAATCGCGCCAATGGAGATTTTTTCTAGGGACGGCGGGGAAC‑TTATCCCCATTTCTCCTGAAATCAAGCTATGTGTGAA‑CTGGAGGGGGGCGCGGTTTGCGGGTTTTTTCCGTGCTGGT‑CGGCGGTGGTTGCTAATTCC‑ATTCTGCCTTTTATGGAACAGGGCGATAAGTTATAGTGGACCAACACCGCCTATTAACCGAGTCAGGATTTCTGTAGCAAAATCCt                                                                                                                                    <  1:210633/251‑1 (MQ=255)
                                                                                                                    aaCACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGc                                                                                                                       >  1:28409/1‑251 (MQ=255)
                                                                                                                     acacGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCg                                                                                                                      >  2:378643/1‑251 (MQ=255)
                                                                                                                     acacGCGACACGTAATCGCGCCAATGGAGAGTTGTACTATGGCAGGCGCGCAACGCTA‑ACCCATTTATCCATACAACAAGCGCTGTGTGAC‑CTGGAAGGTGGCGCAGGCTGAGTGCTATTTCCATGCGGGT‑CGGCAGCGGTTGCTAATTAC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTAAGCAAACTGGGCg                                                                                                                      >  1:392830/1‑251 (MQ=255)
                                                                                                                         gcgACACGTATTCGCGCCAATGGTGAGTTTTTCTATGGACGGCGCGACACGTTA‑CCCCATTTCTCCTTACAACAAGCGATGT‑TGAAACTGGAAGGTGGCGCAGTCTGCGTGTTATTCCCCGGCGGGT‑CGGCAGCG‑TTGGTTATTCCCATTCTTGTTTTTATCGAACAGTGCGATAATGTTTTGATGACCAACCCCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAAc                                                                                                                  <  1:327098/251‑1 (MQ=255)
                                                                                                                                  aaTCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATg                                                                                                         >  1:141413/1‑251 (MQ=255)
                                                                                                                                         ccAATGGAGAGTTTTTCTATGGACGGCGCGGACCGTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGCAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGACCAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTTTGTAGCAAAATCCTCACCAAACAGGGCGTAACGACATCATGGTTTGAt                                                                                                  >  2:220255/1‑251 (MQ=255)
                                                                                                                                                 gaggtgtactctGGACGGCGCGGAACGTTA‑AACCATTTCTCCTTACAACAAGCGAAGTGTGCA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGATAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGTTGATGTGTTGATGACCAAAACCGCCTATGAACCGAGTCAGGAGTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTACCGACATCATGGTTTCATCCGCTGAt                                                                                          >  2:337436/12‑251 (MQ=255)
                                                                                                                                                   gttgtgcaagggggggCGCGGAAAGTTA‑ACCCATTCCTCCTAAAAAAAAGCGATGGGTGAA‑CGGGAAGGGGGCGCAGGGTGCGTGC‑ATTTTCCGGGGGGGTCGGCAGCGGTGGCAAATTCA‑ATTCGTGCTTGTAACGAACGGGGGGAGAATGTGGGGATGACAAACACCGCCGAGGAACCGAGTCAGGATTTCTGGAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGGTTGATCCGCTGATTg                                                                                        <  2:143667/237‑1 (MQ=255)
                                                                                                                                                     tgttCTATGGACGGGGCGGAACGTTA‑ACCCATTTCTCCTTACAACA‑GCGATGTTTGAA‑CGGGAGGGGGGCGCAGTCTGCGTGCTATTCCCCTGCGGGT‑CGGCAGCGGTTGCTATTTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCAGGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGt                                                                                      <  1:44062/250‑1 (MQ=255)
                                                                                                                                                       ttCTATGGACGGCGCGGAACGTTA‑ACCCATTCCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTTTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTg                                                                                     <  1:307661/250‑1 (MQ=255)
                                                                                                                                                             ggACGGCGGGGACCTTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CGGGAGGGGGGCGCAGGCTGCGTTCTTTTTCCCTGGGGGT‑CGGCAGGGGTGGCTAATTCC‑ATTCTTGCTTTTAGCGAACGGGGCGTTCATGTTTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGatat                                                                              <  1:164306/251‑1 (MQ=255)
                                                                                                                                                                 ggCGCGGAACGTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGtt                                                                          >  2:189356/1‑251 (MQ=255)
                                                                                                                                                                         ccGTTA‑CCCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCGGCGGGC‑CGGCAGCGGTTGCTAATCCC‑ATTCTTGCTTTTACCGAACGGGGCGATCATGTGTTGATGCCCACCACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCt                                                                  <  1:272274/250‑1 (MQ=255)
                                                                                                                                                                           gTTA‑CCCCTTTTCTCCTTACAACAAGCGAGGTGTGAA‑CTGAAGGGTGGCGCGGGCGGCTTGTTATTTCCCTGCGGGC‑GGGCAGGGGTTGCAAATCCC‑ATTTTTGTTTTTACCGAACGGGGGGATCATGTGTTGATGCCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTg                                                                 <  2:141413/250‑1 (MQ=255)
                                                                                                                                                                           gTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGc                                                                >  1:300490/1‑251 (MQ=255)
                                                                                                                                                                            gtaaccccttttCTCCTTACAACAAGCGATGGGTGAA‑CTGGAGGGTGGCGCAGGCTGCGTGTTATTTCCCTGCGGGG‑CGGCGGCGGTGTCTAATCCC‑ATTCTTGCTTTTATCGAACGGGGCGATAATGTGTTGAGGACCAACACCGCCTATGAACTGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTg                                                                 <  2:286082/241‑1 (MQ=255)
                                                                                                                                                                             tA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTCAGCATCTGc                                                                >  2:248260/1‑249 (MQ=255)
                                                                                                                                                                                             tACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATc                                               >  1:288276/1‑251 (MQ=255)
                                                                                                                                                                                                                      aagtggccgcAGGCTGCGTGCTATTTCCCGGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTGGCTTTTACCGAACGGGGCGATCATGTGTGGAGGACCACCACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTcc                       <  1:146443/244‑1 (MQ=255)
                                                                                                                                                                                                                       aGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGATTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGAGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCACACACTAAAATCGTGTTTCTGGAATCGCCAGGCTcc                       >  2:310973/1‑250 (MQ=255)
                                                                                                                                                                                                                              gcAGGCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGACCCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGCTATCGTTAAGCATCTGCAGCCACACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCAt               >  2:280532/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gCTGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGa            >  1:345291/1‑250 (MQ=255)
                                                                                                                                                                                                                                    tGCGTGCTATTTCCCTGCGGGT‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACAGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGt         >  2:356052/1‑251 (MQ=255)
                                                                                                                                                                                                                                            attccccTGGGGGT‑CGGCAGCGGTTGCTAATCCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCacgac  <  1:344558/246‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              
AAAGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTA‑ACCCATTTCTCCTTACAACAAGCGATGTGTGAA‑CTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGG‑CGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGAC  >  NC_000913/3152247‑3152730

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