Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 698,183 C→T A98A (GCG→GCA)  unknown ← Aspartate‑semialdehyde dehydrogenase (EC 1.2.1.11)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1698,1830CT98.7% 224.2 / ‑5.0 79A98A (GCG→GCA) unknownAspartate‑semialdehyde dehydrogenase (EC 1.2.1.11)
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (35/43);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

CACCAATAATTCGCGCATATTTTTCGCACCGGCACCACTGGCGGCTTGGTAAGTGGCAACGGAGATCCATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCCGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCTTGTCCTGCTTGAGAAGTGGTAAAGAAGATCGGATTGATATTCGTAAAATCATTTTCCTGCTGCATACGAT  >  1/697943‑698421
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              
cACCAATAATTCGCGCATATTTTTCGCACCGGCACCACTGGCGGCTTGGTAAGTGGCAACGGAGATCCATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCAcc                                                                                                                                                                                                                                       <  1:798981/250‑1 (MQ=255)
cACCAATAATTCGCGCATATTTTTCGCACCGGCACCACTGGCGGCTTGGTAAGTGGCAACGGAGATCCATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCAcc                                                                                                                                                                                                                                       >  1:385690/1‑250 (MQ=255)
  ccAATAATTCGCGCATATTTTTCGCACCGGCACCACTGGCGGCTTGGTAAGTGGCAACGGAGATCCATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAg                                                                                                                                                                                                                                    >  1:489095/1‑251 (MQ=255)
      tacttCGCGCATATTTTTCGCACCGGCACCACTGGCGGCTTGGTAAGTGGCAACGGAGATCCATTCCACTAAAACTTTTTCAAATAAGCCGCCTAACGATATGAGCATCACGCTGACAGTACGGTTACCGCATACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAAc                                                                                                                                                                                                                                <  1:422750/248‑1 (MQ=255)
      taatTCGCGCATATTTTTCGCACCGGCACCACTGGCGGCTTGGTAAGTGGCAACGGAGATCCATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAAc                                                                                                                                                                                                                                <  1:247491/251‑1 (MQ=255)
                atatTTTTCGCACCGGCACCACTGGCGGCTTGGTAAGTGGCAACGGAGATCCATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGcc                                                                                                                                                                                                                      <  1:669092/251‑1 (MQ=255)
                      ttCGCACCGGCACCACTGGCGGCTTGGTAAGTGGCAACGGAGATCCATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGc                                                                                                                                                                                                                <  1:851682/251‑1 (MQ=255)
                            ccGGCACCACTGGCGGCTTGGTAAGTGGCAACGGAGATCCATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATc                                                                                                                                                                                                                                           <  1:528473/218‑1 (MQ=255)
                            ccGGCACCACTGGCGGCTTGGTAAGTGGCAACGGAGATCCATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTaa                                                                                                                                                                                                          >  1:516781/1‑251 (MQ=255)
                                    acTGGCGGCTTGGTAAGTGGCAACGGAGATCCATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGAt                                                                                                                                                                                                  >  1:754695/1‑251 (MQ=255)
                                           gcTTGGTAAGTGGCAACGGAGATCCATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAAc                                                                                                                                                                                                                                <  1:741839/214‑1 (MQ=255)
                                           gcTTGGTAAGTGGCAACGGAGATCCATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACtt                                                                                                                                                                                            >  1:292002/1‑250 (MQ=255)
                                                        cAACGGAGATCCATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATcg                                                                                                                                                                              <  1:374991/251‑1 (MQ=255)
                                                        cAACGGAGATCCATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATc                                                                                                                                                                               <  1:56517/250‑1 (MQ=255)
                                                        cAACGGAGATCCATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATc                                                                                                                                                                               <  1:298502/250‑1 (MQ=255)
                                                               gaTCCATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCtt                                                                                                                                                                       >  1:77882/1‑251 (MQ=255)
                                                                aTCCATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAg                                                                                                                                                                                                                                    <  1:349673/189‑1 (MQ=255)
                                                                  ccATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTgg                                                                                                                                                                     >  1:317730/1‑250 (MQ=255)
                                                                  ccATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTgg                                                                                                                                                                     >  1:114375/1‑250 (MQ=255)
                                                                     ttCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAgg                                                                                                                                                                 >  1:828773/1‑251 (MQ=255)
                                                                            aaaTCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGAt                                                                                                                                                          <  1:777413/251‑1 (MQ=255)
                                                                              aTCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATaa                                                                                                                                                        <  1:189185/251‑1 (MQ=255)
                                                                                   tttCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTcc                                                                                                                                                   <  1:716936/251‑1 (MQ=255)
                                                                                      cAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGa                                                                                                                                                                                                   <  1:445860/200‑1 (MQ=255)
                                                                                           aaaCCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGtttttt                                                                                                                                           >  1:709318/1‑251 (MQ=255)
                                                                                                gCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCAcc                                                                                                                                                                                                                                       <  1:31485/154‑1 (MQ=255)
                                                                                                gCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAAc                                                                                                                                                                                                                                <  1:289349/161‑1 (MQ=255)
                                                                                                gCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTcc                                                                                                                                                   >  1:792536/1‑238 (MQ=255)
                                                                                                gCCTAACGCTATGAGCATAAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCGCCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAAt                                                                                                                                                                                <  1:397545/209‑1 (MQ=255)
                                                                                                      cGCTATGAGCATCAAGCTGACGGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGt                                                                                                                                       <  1:873914/244‑1 (MQ=255)
                                                                                                      cGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGAt                                                                                                                                                          <  1:344296/225‑1 (MQ=255)
                                                                                                      cGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCaa                                                                                                                                <  1:124090/251‑1 (MQ=255)
                                                                                                              gCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGaa                                                                                                                         >  1:219839/1‑250 (MQ=255)
                                                                                                                  cAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCAt                                                                                                                    >  1:559141/1‑251 (MQ=255)
                                                                                                                       tGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTc                                                                                                               <  1:7597/251‑1 (MQ=255)
                                                                                                                         acagtacagtTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAg                                                                                                                                        <  1:480437/224‑1 (MQ=255)
                                                                                                                         acagtacagtTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCaa                                                                                                             >  1:434809/1‑251 (MQ=255)
                                                                                                                         acagtacagtTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCaa                                                                                                             >  1:351155/1‑251 (MQ=255)
                                                                                                                         acagtacagtTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCaa                                                                                                             >  1:804265/1‑251 (MQ=255)
                                                                                                                            gtacagtTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAAc                                                                                                                                                                                                                                >  1:267976/1‑133 (MQ=255)
                                                                                                                            gtacagtTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGtt                                                                                                                                               <  1:319015/214‑1 (MQ=255)
                                                                                                                            gtacagtTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGc                                                                                                                      <  1:618395/239‑1 (MQ=255)
                                                                                                                            gtacagtTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGt                                                                                                          >  1:393380/1‑251 (MQ=255)
                                                                                                                            gtacagtTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGCCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCaatgg                                                                                                          >  1:457024/1‑248 (MQ=255)
                                                                                                                             tacagtTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAg                                                                                                            <  1:279664/248‑1 (MQ=255)
                                                                                                                                 gtTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTa                                                                                                                                                                                                                                  >  1:340955/1‑126 (MQ=255)
                                                                                                                                 gtTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCgg                                                                                                     >  1:490618/1‑251 (MQ=255)
                                                                                                                                 gtTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACG‑CTCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCgg                                                                                                     >  1:805093/1‑250 (MQ=255)
                                                                                                                                         cTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATc                                                                                                                                                                               <  1:657863/169‑1 (MQ=255)
                                                                                                                                          tACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTAc                                                                                            <  1:109242/251‑1 (MQ=255)
                                                                                                                                           aCATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATcgc                                                                                                                                                                             <  1:683715/169‑1 (MQ=255)
                                                                                                                                                            gTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCAtt                                                                                                                   <  1:57721/210‑1 (MQ=255)
                                                                                                                                                               tttCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGc                                                                                                                      >  1:701218/1‑204 (MQ=255)
                                                                                                                                                               tttCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTAccgccg                                                                                       >  1:623615/1‑235 (MQ=255)
                                                                                                                                                               tttCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCtt                                                                       >  1:171859/1‑251 (MQ=255)
                                                                                                                                                               tttCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCtt                                                                       >  1:252841/1‑251 (MQ=255)
                                                                                                                                                                    aGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCTTGTcc                                                                   <  1:390103/250‑1 (MQ=255)
                                                                                                                                                                      ccTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCTTGTCCTg                                                                 >  1:365824/1‑250 (MQ=255)
                                                                                                                                                                         tCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAAc                                                                                                                                                                                                                                >  1:265476/1‑88 (MQ=255)
                                                                                                                                                                            cTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGa                                                                                                                                                              <  1:604309/151‑1 (MQ=255)
                                                                                                                                                                            cTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCTTGTCCTGCTTgaga                                                          >  1:822686/1‑251 (MQ=255)
                                                                                                                                                                             tGATCACGTGTTGATTCACCGGGTCCAACACGTTAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCTTGTCCTGCTTGAGaa                                                         <  1:361426/251‑1 (MQ=255)
                                                                                                                                                                              gATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGa                                                                                                                                                              >  1:221284/1‑149 (MQ=255)
                                                                                                                                                                              gATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGttt                                                                                                                                              >  1:690001/1‑165 (MQ=255)
                                                                                                                                                                                       ttGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCTTGTCCTGCTTGAGAAGTGGTAaag                                                <  1:709234/250‑1 (MQ=255)
                                                                                                                                                                                        tGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCTTGTCCTGCTTGAGAAg                                                        <  1:336516/241‑1 (MQ=255)
                                                                                                                                                                                                  ggTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGa                                                                                                                                                              <  1:898007/129‑1 (MQ=255)
                                                                                                                                                                                                  ggTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCTTGTCCTGCTTGAGAAGTGGTAAAGAAGATCGGATTg                                    >  1:499778/1‑251 (MQ=255)
                                                                                                                                                                                                   gTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAg                                                                          >  1:590692/1‑212 (MQ=255)
                                                                                                                                                                                                   gTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCTTGTCCTGCTTGAGAAg                                                        <  1:343639/230‑1 (MQ=255)
                                                                                                                                                                                                   gTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCTTGTCCTGCTTGAGAAg                                                        >  1:827323/1‑230 (MQ=255)
                                                                                                                                                                                                      cAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCTTGTCCTGCTTGAGAAGTGGTAAAGAAGATCGGATTGatat                                <  1:565010/251‑1 (MQ=255)
                                                                                                                                                                                                           cGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGtt                                                                                                                                               <  1:568474/135‑1 (MQ=255)
                                                                                                                                                                                                           cGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCTTGTcc                                                                   <  1:417114/211‑1 (MQ=255)
                                                                                                                                                                                                                          cttcCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCTTGTCCTGCTTGAGAAGTGGTAaag                                                >  1:421826/1‑215 (MQ=255)
                                                                                                                                                                                                                            tcCATCCGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCTTGTCCTGCTTGAGAAATGGTAAAGAAGATCGGATTGATATTCGTAAAATCATTTTc                <  1:897577/245‑1 (MQ=255)
                                                                                                                                                                                                                            tcCATACGCAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCTTGTCCTGCTTGAGAAGTGGTAAAGAAGATCGGATTGATATTCGTAAAATCATTTTCctg             >  1:766969/1‑248 (MQ=255)
                                                                                                                                                                                                                                    cAATGCGGACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCTTGTCCTGCTTGAGAAGTGGTAAAGAAGATCGGATTGATATTCGTAAAATCATTTTCCTGCTGCATACGAt  <  1:188916/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gACGCTGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCTTGTCCTGCTTGAGAAGTGGTAAAGAAGATCGGATTGata                                 <  1:140180/213‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              
CACCAATAATTCGCGCATATTTTTCGCACCGGCACCACTGGCGGCTTGGTAAGTGGCAACGGAGATCCATTCCACTAAATCTTTTTCAAATAAACCGCCTAACGCTATGAGCATCAAGCTGACAGTACAGTTACCGCCTACATAGGTTTTCACGCCGTTTTTCAAGCCTTCGCTGATCACGTGTTGATTCACCGGGTCCAACACGATAATCGCGTCTTCTTCCATACGCAATGCGGACGCCGCATCCACCCAGTAACCGCTCCAGCCGGTGGCTTTTAATTTCGGATAAACTTCGTTGGTGTAATCGCCGCCTTGGCAGGTGACGATAATGTCCAGTTTTTTCAGTTCTTCAATGTCGAACGCATTTTTCAAGGTGCCGGCGTCTTTACCGCCGAATGTCGGTGCAGCTTGTCCTGCTTGAGAAGTGGTAAAGAAGATCGGATTGATATTCGTAAAATCATTTTCCTGCTGCATACGAT  >  1/697943‑698421

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