Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,390,331 C→T Q10Q (CAG→CAA psd ← phosphatidylserine decarboxylase proenzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,390,3310CT90.7% 110.5 / 2.0 43Q10Q (CAG→CAApsdphosphatidylserine decarboxylase proenzyme
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (20/19);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

CGGCAGGCATGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCTTCCTCAACCGCTTCCCGG  >  NC_000913/4390093‑4390574
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
cGGCAGGCAGGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTGCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAg                                                                                                                                                                                                                                          <  2:308532/250‑1 (MQ=255)
 ggcaggcaTGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGtt                                                                                                                                                                                                                                        >  2:81055/1‑251 (MQ=255)
    aggcaTGACCAGTACATTCGGAGCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGTTTTTGCGCCTCTTTCATGACGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGTCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTaa                                                                                                                                                                                                                                     <  1:342414/251‑1 (MQ=255)
          gACCAGTAAATTCGGATCGGTATCGATTGGGGGTTCTTCGTCACGCAGCGGACGGACAAAGAATTTGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGaa                                                                                                                                                                                                                                <  2:392058/250‑1 (MQ=255)
           aCCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAAtt                                                                                                                                                                                                                              >  1:25966/1‑251 (MQ=255)
           aCCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAAtt                                                                                                                                                                                                                              >  1:82447/1‑251 (MQ=255)
                               aTCGATTGGGCGTACTTCGTCACGCAGCGGCCGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGTCCTTGTAGTATTTAACGAACAGATCGTTCACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGc                                                                                                                                                                                                          >  1:459469/1‑251 (MQ=255)
                                 cGATTGGGGGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAGCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCa                                                                                                                                                                                                        <  2:277245/251‑1 (MQ=255)
                                      gggCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGtt                                                                                                                                                                                                   >  2:197187/1‑251 (MQ=255)
                                              ttCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGCAAGTTTAAATGAATTTAACAAGGTAGCCTCAAGGCCa                                                                                                                                                                                                        >  2:45941/1‑238 (MQ=255)
                                              ttCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCa                                                                                                                                                                                                        <  1:45941/238‑1 (MQ=255)
                                                    ccGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTTTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATACCCATTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGAt                                                                                                                                                                                     <  1:54602/250‑1 (MQ=255)
                                                                gacAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATg                                                                                                                                                                         <  1:238144/251‑1 (MQ=255)
                                                                 acAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGc                                                                                                                                                                        >  2:98419/1‑251 (MQ=255)
                                                                  cAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTATGCGCCTATTTCTTGTCGACATTGTAGTATTTAACGAACAGATCGATAACAAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACGCGCCAGGCGAGTAAGCCATCGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAAATTAACACGGTAGCCTCAGGGCCATTGTTTTGTAGTTACTGATCCGGcataactgac                                                                                                                                                                       >  2:247029/1‑242 (MQ=255)
                                                                   aaaGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCg                                                                                                                                                                      >  1:112308/1‑251 (MQ=255)
                                                                   aaaGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCg                                                                                                                                                                      >  1:447091/1‑251 (MQ=255)
                                                                                 aGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGaaaaaa                                                                                                                                                         >  1:368887/1‑250 (MQ=255)
                                                                                      cGGTAGCTGGTGGTGTCCTGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCATTTTTGTCAGCCATCCTGCCCTCTTGCTTGCTCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGa                                                                                                                                                   <  1:97100/251‑1 (MQ=255)
                                                                                           gCTGGCGGTGGCCGGTTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGAtt                                                                                                                                              <  2:439292/251‑1 (MQ=255)
                                                                                              ggCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTa                                                                                                                                           <  2:410010/251‑1 (MQ=255)
                                                                                                      ccGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACgg                                                                                                                                    <  1:197187/250‑1 (MQ=255)
                                                                                                                     ctTTCATGTCGTCCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTgtc                                                                                                                      <  1:484306/249‑1 (MQ=255)
                                                                                                                                cccTTGTAGTATTTAACGAACAGATCGATAACCCGTTTTGTCAGCCACCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATc                                                                                                         <  2:458053/250‑1 (MQ=255)
                                                                                                                                aCCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATc                                                                                                         >  2:360858/1‑251 (MQ=255)
                                                                                                                                     gtagtaTTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGCTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGaa                                                                                                     >  1:236008/1‑250 (MQ=255)
                                                                                                                                       agtaTTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGaaaa                                                                                                   <  1:98419/250‑1 (MQ=255)
                                                                                                                                          aTTTAACAGACAGTTCGGTAAGCAGTGTTGGCAGCCATCCGGACCGCTTGCTTGCTCCCAAACTCGCCAGACGAGTAAGACATAGTTTCGACCGACTATATTCAAGCGAAAGTTTAAATTACTTTAACAAGGGAGCCTCCAGGCCATGGTTTGGTCGTTCCTGTTCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGAAGCTTCTAACGACGGTTAGCTTACTTGTCAGTCAAACGTAGCAGaaaattgt                                                                                               >  2:353532/1‑247 (MQ=255)
                                                                                                                                               aCGAACATCTCGATAACCAGTTTTGCCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTAcg                                                                                          <  2:74879/251‑1 (MQ=255)
                                                                                                                                                                        tCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCAGGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAAAAAGGTAGCCTCCAGACCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAGGGGGACGATTCTAACGACGGTTGGCTTAATTGTCAGTCGTCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTcc                                                                 >  1:241781/1‑251 (MQ=255)
                                                                                                                                                                         aacccctccggCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCTGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCa                                                                <  2:459469/241‑1 (MQ=255)
                                                                                                                                                                          aGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCtt                                                                    >  2:444948/1‑246 (MQ=255)
                                                                                                                                                                          aGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCtt                                                                    <  1:444948/246‑1 (MQ=255)
                                                                                                                                                                                              gtcgcCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTc                                             <  2:447091/247‑1 (MQ=255)
                                                                                                                                                                                                 gcCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAAc                                        <  1:136907/251‑1 (MQ=255)
                                                                                                                                                                                                 gcCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTGCCAGAATACGGTGATAGTTTTCGAAAc                                        <  1:178316/251‑1 (MQ=255)
                                                                                                                                                                                                             aaGGCGAGTAAGCCATAGTTTATGCAGAATGTATTGTAGCGAACGTTTAAATGAATTTACCAAGGTAGCCTACAGGCCATTGTTTAGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAATAAAGTGGACGATTCTAACGACGGTTAGCTTACTTGTCCATCATACGTATCAGCAAAGTTTTTACGTGTTTTTACCTGCGCCATCCTTTCCAGAATACGGTGATAGTTTTCGAACCGGGTTTccacgt                            >  2:432716/2‑247 (MQ=255)
                                                                                                                                                                                                                      aaGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCtt                                                                    >  1:202643/1‑202 (MQ=255)
                                                                                                                                                                                                                      aaGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGGTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCtt                                                                    <  2:202643/202‑1 (MQ=255)
                                                                                                                                                                                                                       aGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTACCTGATCCGGCCTACATGCCGGCTCCTGAAACAAAGGAGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGCGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCTTc                  >  2:415720/1‑251 (MQ=255)
                                                                                                                                                                                                                          cATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGCTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCTTCCTc               <  2:236008/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCTTCCTCAACCGCTTcc     >  1:38114/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCAGGCCTACATGCCGGATCCTGAAAAAAAGGGGCCGATTCTAACGTCGGTTAGCTTAATTGTCAGTCATCCGTATCAGAACAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCCGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCTTCCTCAACCGCTTCCCgg  >  2:475216/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
CGGCAGGCATGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCTTCCTCAACCGCTTCCCGG  >  NC_000913/4390093‑4390574

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