Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,935,391 T→C V269A (GTA→GCA)  unknown → Phosphoserine aminotransferase (EC 2.6.1.52)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,935,3910TC98.5% 191.0 / ‑2.7 66coding (806/1083 nt)unknownPhosphoserine aminotransferase (EC 2.6.1.52)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (42/23);  minor base . (1/0);  total (43/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

ATTTATGCCGGTGCACAAAAAAACCTCGGTCCGGCAGGGATTACCATCGTGATCGTACGTGACGATCTCATCGGTAAGGCACGTAAAGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGTAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACAAAGTGCTTGGCGGTATGCGCGCGTCCATTTATAACGCCATGCCGATTGAAGGCGTACAAGCATTGATTGC  >  1/1935153‑1935634
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
aTTTATGCCGGTGCACAAAAAAACCTCGGTCCGGCAGGGATTACCATCGTGATCGTACGTGACGATCTCATCGGTAAGGCACGTAAAGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGaaa                                                                                                                                                                                                                                           >  1:161656/1‑249 (MQ=255)
          gTGCACAAAAAAACCTCGGTCCGGCAGGGATTACCATCGTGATCGTACGTGACGATCTCATCGGTAAGGCACGTAAAGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGat                                                                                                                                                                                                                               >  1:799891/1‑251 (MQ=255)
          gTGCACAAAAAAACCTCGGTCCGGCAGGGATTACCATCGTGATCGTACGTGACGATCTCATCGGTAAGGCACGTAAAGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGa                                                                                                                                                                                                                                <  1:681812/250‑1 (MQ=255)
          gTGCACAAAAAAACCTCGGTCCGGCAGGGATTACCATCGTGATCGTACGTGACGATCTCATCGGTAAGGCACGTAAAGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGa                                                                                                                                                                                                                                <  1:627435/250‑1 (MQ=255)
                   aaaaCCTCGGTCCGGCAGGGATTACCATCGTGATCGTACGTGACGATCTCATCGGTAAGGCACGTAAAGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGAc                                                                                                                                                                                                                      <  1:697893/251‑1 (MQ=255)
                                             aTCGTGATCGTACGTGACGATCTCATCGGTAAGGCACGTAAAGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAAtgt                                                                                                                                                                                            >  1:334460/1‑251 (MQ=255)
                                              tCGTGATCGTACGTGACGATCTCATCGGTAAGGCACGTAAAGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAAtgtg                                                                                                                                                                                           >  1:536954/1‑251 (MQ=255)
                                              tCGTGATCGTACGTGACGATCTCATCGGTAAGGCACGTAAAGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCAGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCGATAAtgtg                                                                                                                                                                                           >  1:447267/1‑251 (MQ=255)
                                                 tGATCGTACGTGACGATCTCATCGGTAAGGCACGTAAAGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGG‑CAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGa                                                                                                                                                                                       >  1:286679/1‑251 (MQ=255)
                                                  gATCGTACGTGACGATCTCATCGGTAAGGCACGTAAAGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGat                                                                                                                                                                                                                               >  1:701342/1‑211 (MQ=255)
                                                     cGTACGTGACGATCTCATCGGTAAGGCACGTAAAGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGa                                                                                                                                                                                                                                <  1:507345/207‑1 (MQ=255)
                                                      gTACGTGACGATCTCATCGGTAAGGCACGTAAAGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAg                                                                                                                                                                                    >  1:153785/1‑250 (MQ=255)
                                                                      tCGGTAAGGCACGTAAAGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCGAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTaatg                                                                                                                                                                   >  1:436348/1‑251 (MQ=255)
                                                                                  gTAAAGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCttt                                                                                                                                                       >  1:725245/1‑251 (MQ=255)
                                                                                     aaGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCt                                                                                                                                                                                                         >  1:498110/1‑198 (MQ=255)
                                                                                     aaGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTAc                                                                                                                                                                                                       >  1:742809/1‑200 (MQ=255)
                                                                                     aaGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCt                                                                                                                                                         >  1:823645/1‑246 (MQ=255)
                                                                                     aaGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTac                                                                                                                                                     >  1:97742/1‑250 (MQ=255)
                                                                                     aaGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTac                                                                                                                                                     >  1:213213/1‑250 (MQ=255)
                                                                                          acacCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTAttt                                                                                                                                                                                                                          >  1:881/1‑176 (MQ=255)
                                                                                           cacCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTaccacccg                                                                                                                                               >  1:515237/1‑248 (MQ=255)
                                                                                               gTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCgatg                                                                                                                                          >  1:195965/1‑251 (MQ=255)
                                                                                                        gAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTaccac                                                                                                                                                  >  1:292132/1‑234 (MQ=255)
                                                                                                        gAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGAtgaactg                                                                                                                                  <  1:185333/250‑1 (MQ=255)
                                                                                                                                    tCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTaccac                                                                                                                                                  >  1:850125/1‑206 (MQ=255)
                                                                                                                                    tCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCAccg                                                                                                      >  1:108638/1‑250 (MQ=255)
                                                                                                                                          aTCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTAttt                                                                                                                                                                                                                          <  1:849570/128‑1 (MQ=255)
                                                                                                                                          aTCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAgg                                                                                               >  1:138204/1‑251 (MQ=255)
                                                                                                                                          aTCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAgg                                                                                               >  1:706341/1‑251 (MQ=255)
                                                                                                                                           tCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGt                                                                                              >  1:512446/1‑251 (MQ=255)
                                                                                                                                              acacCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTTTTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTa                                                                                           >  1:500136/1‑251 (MQ=255)
                                                                                                                                                cccccGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACa                                                                                         <  1:135764/250‑1 (MQ=255)
                                                                                                                                                  cccGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAgg                                                                                       <  1:58917/251‑1 (MQ=255)
                                                                                                                                                  cccGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAgg                                                                                       <  1:391374/251‑1 (MQ=255)
                                                                                                                                                        gACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTaa                                                                                 >  1:230717/1‑251 (MQ=255)
                                                                                                                                                         aCTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTaaa                                                                                <  1:619422/251‑1 (MQ=255)
                                                                                                                                                           tttCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCGCCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAg                                                                              <  1:397438/251‑1 (MQ=255)
                                                                                                                                                                  tGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTaccac                                                                                                                                                  <  1:60167/176‑1 (MQ=255)
                                                                                                                                                                  tGGTACTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACaa                                                                       >  1:56943/1‑251 (MQ=255)
                                                                                                                                                                    gTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTgg                                                                                                                                                                                          >  1:688689/1‑134 (MQ=255)
                                                                                                                                                                    gTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTCCCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACAAAg                                                                     >  1:453879/1‑251 (MQ=255)
                                                                                                                                                                    gTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCggtggtaaag                                                                                                      >  1:220067/1‑209 (MQ=255)
                                                                                                                                                                    gTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAg                                                                                        <  1:622327/232‑1 (MQ=255)
                                                                                                                                                                    gTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACAAAg                                                                     >  1:556337/1‑251 (MQ=255)
                                                                                                                                                                     tATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGACACAAAg                                                                     <  1:645780/250‑1 (MQ=255)
                                                                                                                                                                       tttATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACAAAGTg                                                                   >  1:406017/1‑250 (MQ=255)
                                                                                                                                                                        ttATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACAAAGTg                                                                   >  1:102395/1‑249 (MQ=255)
                                                                                                                                                                            gtgCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACAAAGTGCtt                                                                >  1:693100/1‑248 (MQ=255)
                                                                                                                                                                            gtgCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACAAAGTGCTTgg                                                              >  1:588854/1‑250 (MQ=255)
                                                                                                                                                                              gCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAg                                                                                        <  1:37918/222‑1 (MQ=255)
                                                                                                                                                                                   ggTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGAtgaactg                                                                                                                                  >  1:413170/1‑175 (MQ=255)
                                                                                                                                                                                      aTTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTaat                                                                                                                                                                    >  1:85266/1‑138 (MQ=255)
                                                                                                                                                                                      aTTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACAAAGTGCTTGGCGGTATgcgcg                                                   >  1:385689/1‑251 (MQ=255)
                                                                                                                                                                                      aTTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACAAAGTGCTTGGCGGTATgcgcg                                                   <  1:837156/251‑1 (MQ=255)
                                                                                                                                                                                      aTTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAAACGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAg                                                                              >  1:868229/1‑224 (MQ=255)
                                                                                                                                                                                                 tatcggctatcggcGGCTTGCCGGAAATCGAAAAACGCACCCCGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTCCCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTATAAGGGCACAAAGTGATTGGCGGTATGCGCGCGTCCATTtat                                        <  1:496374/251‑1 (MQ=255)
                                                                                                                                                                                                 tatcggctatcggcGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACAAAGTGCTTGGCGGTATGCGCGCGTCCATTtat                                        <  1:884727/251‑1 (MQ=255)
                                                                                                                                                                                                   tcggctatcggcGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACAAAGTGCTTGGCGGTATGCGCGCGTCCATTTATaa                                      <  1:316587/251‑1 (MQ=255)
                                                                                                                                                                                                    cggctatcggcGGCTTGCCGGAAATCGAAAATCGCATCGGGGCAGAAGCGAAGCGGCTTTGCGATTATTTGGAAGGCAGCGATTTCTACCGCGATAGTGTGGCGAAAGCCACCCGCAGTTTAATGAATGTATCCTTTTCCACCGGCAATGATGAAGTGAACGCCAAATTCATCGCCGCTGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACAAAGTGCTTGGCGGTATGCGCGCGTCCATTTATAAc                                     <  1:686750/251‑1 (MQ=255)
                                                                                                                                                                                                    cggctatcggcGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACAAAGTGCTTGGCGGTATGCGCGCGTCCATTTATAAc                                     <  1:413963/251‑1 (MQ=255)
                                                                                                                                                                                                      gctatcggcGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAAtgt                                                                                                                                                                                            >  1:892157/1‑98 (MQ=255)
                                                                                                                                                                                                      gctatcggcGGCTTGCCGGAAATCGAAAAACGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACAAAGTGCTTGGCGGTATGCGCGCGTCCATTTATAACg                                    >  1:235154/1‑250 (MQ=255)
                                                                                                                                                                                                                               aaaaaCGCAACCTGGCAAAAGCGAAACTGCTTTACGATGATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCGCCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCGTCGCGGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACAAAGTGCTTGGCGGTATGCGCGCGTCCATTTATAACGCCATGCCGATTGAAGGCGTACAAGCa          <  1:490543/251‑1 (MQ=255)
                                                                                                                                                                                                                                aaaaCGCAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACAAAGTGCTTGGCGGTATGCGCGCGTCCATTTATAACGCCATGCCGATTGAAGGCGTACAAGCat         <  1:755448/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cAACCTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTacc                                                                                                                                                    >  1:540519/1‑106 (MQ=255)
                                                                                                                                                                                                                                         ccTGGCAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACAAAGTGCTTGGCGGTATGCGCGCGTCCATTTATAACGCCATGCCGATTGAAGGCGTACAAGCATTGATTGc  >  1:728689/1‑249 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
ATTTATGCCGGTGCACAAAAAAACCTCGGTCCGGCAGGGATTACCATCGTGATCGTACGTGACGATCTCATCGGTAAGGCACGTAAAGCCACACCGTCCATTTGGAATTACGAAATTCAGAAAAACAGCGATTCCATGATCAACACCCCGCCGACTTTCGCATGGTATTTATGTGCCTTGGTATTCAAACATTTATCGGCTATCGGCGGCTTGCCGGAAATCGAAAAACGCAACCTGGTAAAAGCGAAACTGCTTTACGATTATTTGGACGGCAGCGATTTCTACCGCAATAATGTGGCGAAAGCCAACCGCTCTTTAATGAATGTCACCTTTACCACCGGCGATGATGAACTGAACGCCAAATTCATCGCCGAAGCCACCGCCGCAGGTTTACAGGCGTTAAAAGGGCACAAAGTGCTTGGCGGTATGCGCGCGTCCATTTATAACGCCATGCCGATTGAAGGCGTACAAGCATTGATTGC  >  1/1935153‑1935634

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