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,4940GT93.9% 147.4 / ‑5.4 51A87S (GCG→TCG) metCcystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  major base T (29/17);  minor base C (1/2);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTCCGGCGATTG  >  NC_000913/3152247‑3152743
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         
aaaGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGCACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGgtcgg                                                                                                                                                                                                                                                        >  2:50344/1‑251 (MQ=255)
        aTACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGtt                                                                                                                                                                                                                                                <  2:241303/251‑1 (MQ=255)
                 tGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTcc                                                                                                                                                                                                                                       >  2:354250/1‑251 (MQ=255)
                          cAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCt                                                                                                                                                                                                                              <  1:99206/251‑1 (MQ=255)
                                   gcaAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGa                                                                                                                                                                                                                      >  1:235482/1‑250 (MQ=255)
                                    caAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAAc                                                                                                                                                                                                                    <  2:184225/251‑1 (MQ=255)
                                    caAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTATTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAAc                                                                                                                                                                                                                    <  2:225073/251‑1 (MQ=255)
                                            acacTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGAt                                                                                                                                                                                                            >  1:177101/1‑251 (MQ=255)
                                               ctctCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCAt                                                                                                                                                                                                         >  2:13597/1‑251 (MQ=255)
                                                   cGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATgtgt                                                                                                                                                                                                     <  1:298599/251‑1 (MQ=255)
                                                    ggCGCGGTAAATCGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACATAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAGCTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCAAACAGGGCGATCATGTGtt                                                                                                                                                                                                    >  2:221829/1‑251 (MQ=255)
                                                         ggTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTtgatg                                                                                                                                                                                                >  1:369157/1‑250 (MQ=255)
                                                          gTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGAc                                                                                                                                                                                              <  1:224369/251‑1 (MQ=255)
                                                                aGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAaca                                                                                                                                                                                         <  1:206113/250‑1 (MQ=255)
                                                                aGCGTGATTCAGCGCGCTTCTTAGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAacac                                                                                                                                                                                        >  1:108479/1‑251 (MQ=255)
                                                                               gcctcctcggtcgTCTTTGCCAGTGAAGAAGCCAAAAAACACCCGACACGTAATCGCGCAAATGGAGAGATGTTCTATGGACGCCACGGAACGATAACCCATTTCTCCTTACAACAAGCGATGTGTGCACTGGTAGGCGGCGCAGCCTGCGTGCTATTTCCCTCCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTCTCGAACAGTGCGATCATGTGTTGATGACCAACCACGCCTCTGACCCGa                                                                                                                                                                          >  2:302608/13‑250 (MQ=255)
                                                                                    ttcGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTATTGATGACCAACACCGCCTATGAACCGAGTCAgg                                                                                                                                                                    >  2:97707/1‑251 (MQ=255)
                                                                                      cgcgggTCTTTGACAGTGTAGAAGCCAAAAAACACGGGACACGTAAGCGCGCCAATGGAGAGTGGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAAAAAGCGATGTGTGAACTGGAAGGGGGCGCAGGCTGCGTGATATTCCCCGGCGGGTCGGCAGCGTTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTGGATGACCAACACCGCCTATGAACCGAGTCAGGAt                                                                                                                                                                  <  1:387853/247‑1 (MQ=255)
                                                                                         tGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTc                                                                                                                                                               >  2:52605/1‑251 (MQ=255)
                                                                                          ggTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGCGTTGTTCTATGGACGGCGCGGACCGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCt                                                                                                                                                              >  2:34086/1‑251 (MQ=255)
                                                                                             ttttGACAGGGTAAAAGAAAAAAAACACGCGACACGTAATCGCGCCAAGGGAGAGTTGTTCTATGGACGGCGCGGAACGTTACCCCATTTCTCCTTACAACAAGGGATGTGGGAACTGGAAGGTGGCGCAGGCTGCGTGCTTTTCCCCGGCGGGTGGGCAGCGTTTGCTAATTCCATTCTTTCTTTTATCGAACAGGGCGATCTTGTGTTGAGTACCAACACCGCCTATGAACCGAGTCAGGATTTCTGt                                                                                                                                                            <  2:29544/249‑1 (MQ=255)
                                                                                                               cAAAAAACACGCGACACGTAACCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTCCAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAAc                                                                                                                                         >  2:377144/1‑251 (MQ=255)
                                                                                                                     cccggcgACCCGTAACCGCCCCAAGGGAGAGTTTTTCAATGGCGGGCGGGGAAGGTAACCCCATTCCTCCTTCCAACAGGCGATGTGGGACCGGGAGGGGGGCGCAGGCTGGGTGTTTTTCCCCTGCGGGTCGGCGGCGTTTGCTAATTCCATTCTTGCTTTTATCGAACGGGGCGATCTTGTGTTGATGACAAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCg                                                                                                                                   <  2:9009/247‑1 (MQ=255)
                                                                                                                     acacGCAAAACGTAATCGAGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGCTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGGGGCGCAGGCTGCGTGCTATTTCCCTACGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGCGCGATCATGCGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTggccg                                                                                                                                   >  2:450529/1‑248 (MQ=255)
                                                                                                                       ccgcgACACGTAATCGCGCCAAGGGAGAGTTGTTCTAGGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACGGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCGGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTATTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTa                                                                                                                                 <  1:97707/250‑1 (MQ=255)
                                                                                                                            acacGTAATCGCGCCAATGGAGGGTTGTTCTATGGACGGCGCGGGACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGCTCATGTGTTGATGACCAACACCGGCTATGAACCGAGTCAGGCTTTCTGTTGCAAAATCCTCAGCAAACTGGGCGTACCGAc                                                                                                                            >  2:254657/1‑251 (MQ=255)
                                                                                                                               gttAATCGCGACATTGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAAAAAGCGATGTGTGAACGGGAGGGTGGCGCGGGCTGCGTGTTATTTCCCGGCGGGGCGGCAGCGGTTTCTAATTCCATTCTTCCTTTAATCGAACAGGGCGATCTTGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcatc                                                                                                                         <  2:384212/249‑1 (MQ=255)
                                                                                                                                  aaTCGCGCCAATGGCGAGTAGTTCTATGGTCGGCGCGGACCGTTAACCCATTTCCCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCACGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGCGTTGATGAACACCACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAACTCCTCAGCAAACTGGGCGTAACGACATCATg                                                                                                                      >  2:425426/1‑251 (MQ=255)
                                                                                                                                    ccgcgCCAATGGAGAGTTTTTCTATGGCCGGCGCGAAACTTTACCCCATTTCTCCTTACAACAAGCGATGTGTGAACGGGAAGGTGGCGCAGGCTGCGTGCTATTCCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTGGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGt                                                                                                                    <  2:24319/250‑1 (MQ=255)
                                                                                                                                     cgcgCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGtt                                                                                                                   >  1:151850/1‑251 (MQ=255)
                                                                                                                                       cgcCAATGGCGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTg                                                                                                                 >  2:210894/1‑251 (MQ=255)
                                                                                                                                              ggAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACCGGAAGCTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGCCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCt                                                                                                          >  1:457705/1‑251 (MQ=255)
                                                                                                                                              ggAGAGTTGTTCTAGGGACGGCGCGGAATGTTAACCCATTTCTCCTTACAACAAGCGAGGTGTGAACTGGAGGGGGGCGAAGGCTGCGTGCTTTTCCCCTGCGGGTGGGCAGCGTTTGTGATTTCCATTCTTTCTTTTAGCGAACGGGGCTATCAGTTTTTGAGGACCAACACCGCCTATGAACCGAGTCAGGATGTCTGTAGCAAAATCCGCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCt                                                                                                          <  1:29593/251‑1 (MQ=255)
                                                                                                                                                  agTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGAtt                                                                                                      >  2:259514/1‑251 (MQ=255)
                                                                                                                                                   gttgttCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTg                                                                                                     >  1:197377/1‑251 (MQ=255)
                                                                                                                                                      gttCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTg                                                                                                  >  2:237932/1‑251 (MQ=255)
                                                                                                                                                             ggACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGatat                                                                                           >  1:160287/1‑251 (MQ=255)
                                                                                                                                                              gcgggagcggaccgtgaccccATTTCACCTTAAAACAAGCGATGTGGGAACTGGAGGGTGGCGCAGTCGGCTGGTTATTCCCCGGCGGGCGGGCAGCGGTTGCTAATTCCATTCTTGCTTTTACCGAACAGGGCGATCATGTGTGGATGACAAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGatata                                                                                          <  2:81798/233‑2 (MQ=255)
                                                                                                                                                                      ggAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTACTTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGc                                                                                   >  2:306254/1‑250 (MQ=255)
                                                                                                                                                                             tAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGc                                                                             >  1:319583/1‑249 (MQ=255)
                                                                                                                                                                             aaccccATTTCTCCTTACAACAAGGGAGGGGGGAACGGGAAGGTGGCGCAGGCTGCGGGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATCCCATTCTTGCTTTTATCGAACAGGGCGATCAGGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAg                                                                           <  1:139006/248‑1 (MQ=255)
                                                                                                                                                                              aaCCCATTTCTCCTTACAACAAGCGATGTGTGAGCTGGAAGGTGGCGCAGGCTGCTTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGAGGCCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAg                                                                           <  1:247573/250‑1 (MQ=255)
                                                                                                                                                                                       ctcCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAACCACt                                                                  >  2:380794/1‑250 (MQ=255)
                                                                                                                                                                                         ccTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTaa                                                                >  2:377014/1‑250 (MQ=255)
                                                                                                                                                                                                     cGATGTGTGAACTGGAAGGTGGCGCGGGCTGCGTGCTATTTCCCTGCGGGTCGGCGGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTg                                                   <  1:366401/251‑1 (MQ=255)
                                                                                                                                                                                                              aacgagAAGGTGGCGCAGGCTGCTTCCTATTTCCCGTCGGGCAGGCTGCGGTGGTTAATTCCATTCTTGCTTTCACCGAACAGGGCGATAATTTTTTAATGCCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGTGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGACTCGCCa                                          <  2:93946/246‑1 (MQ=255)
                                                                                                                                                                                                                 tGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGc                                       >  2:393658/1‑251 (MQ=255)
                                                                                                                                                                                                                   gAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTc                                     >  2:137123/1‑251 (MQ=255)
                                                                                                                                                                                                                                gCTGCGTGCTTTTCCCCGGCGGGTGGGCAGGGTTTTCTAATTCCATTCTTTCTTTTATCAAAAAGGGCAATCATGTGTTGAGGCCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGaa                        <  2:249460/251‑1 (MQ=255)
                                                                                                                                                                                                                                         tATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCacgac               >  1:32845/1‑251 (MQ=255)
                                                                                                                                                                                                                                           tttCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGt             >  1:177036/1‑251 (MQ=255)
                                                                                                                                                                                                                                             tCCCTGCGGGGCGGCAGCGGTCGCTAATTCCCTTCTTCCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACGCCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTc           >  2:225862/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      ggcggcggcGGTTGTTAATTCCATTCGTGTTTTTAGCGAACGGGGCGAGCATGTGTTGATGAGAAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGACAGGCTCCATCACCATGGAAGTCCACGACGGTCCGGCGAttg  <  2:398260/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         
AAAGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTCCGGCGATTG  >  NC_000913/3152247‑3152743

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