Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 607,933 G→A L792F (CTC→TTC)  unknown ← Ribonuclease E (EC 3.1.26.12)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1607,9330GA98.7% 262.6 / ‑3.4 78L792F (CTC→TTC) unknownRibonuclease E (EC 3.1.26.12)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (39/38);  minor base T (1/0);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

CACTTCAGCTGTCGCTTCCATTGGAGCGGATACCACATTGTTATTTAAGCGTTGTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAGGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATCGCGCCCTTTATCGACTTCAGCATTATCATCAAGATCAATCACGTCCGTTACTACAATCTCTTCTT  >  1/607688‑608175
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
cACTTCAGCTGTCGCTTCCATTGGAGCGGATACCACATTGTTATTTAAGCGTTGTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAATGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTAt                                                                                                                                                                                                                                               >  1:441179/1‑251 (MQ=255)
        cTGTCGCTTCCATTGGAGCGGATACCACATTGTTATTTAAGCGTTGTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCttt                                                                                                                                                                                                                                        <  1:71436/250‑1 (MQ=255)
              cTTCCATTGGAGCGGATACCACATTGTTATTTAAGCGTTGTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTttgattga                                                                                                                                                                                                                                 >  1:490568/1‑251 (MQ=255)
                 ccATTGGAGCGGATACCACATTGTTATTTAAGCGTTGTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACg                                                                                                                                                                                                                              >  1:638681/1‑251 (MQ=255)
                           ggATACCACATTGTTATTTAAGCGTTGTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGccc                                                                                                                                                                                                                      <  1:127258/249‑1 (MQ=255)
                            gATACCACATTGTTATTTAAGCGTTGTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCAt                                                                                                                                                                                                                   >  1:490039/1‑251 (MQ=255)
                              tACCACATTGTTATTTAAGCGTTGTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGt                                                                                                                                                                                                                 <  1:751712/251‑1 (MQ=255)
                              tACCACATTGTTATTTAAGCGTTGTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGt                                                                                                                                                                                                                 <  1:80524/251‑1 (MQ=255)
                              tACCACATTGTTATTTAAGCGTTGTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGt                                                                                                                                                                                                                 <  1:866548/251‑1 (MQ=255)
                              tACCACATTGTTATTTAAGCGTTGTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGt                                                                                                                                                                                                                 <  1:826722/251‑1 (MQ=255)
                                ccACATTGTTATTTAAGCGTTGTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCg                                                                                                                                                                                                               >  1:29092/1‑251 (MQ=255)
                                   caTTGTTATTTAAGCGTTGAAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATc                                                                                                                                                                                                            <  1:597409/251‑1 (MQ=255)
                                        ttATTTAAGCGTTGTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCacac                                                                                                                                                                                                       <  1:175668/251‑1 (MQ=255)
                                        ttATTTAAGCGTTGTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCacac                                                                                                                                                                                                       <  1:95887/251‑1 (MQ=255)
                                           tttAAGCGTTGTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGt                                                                                                                                                                                                                 >  1:228351/1‑238 (MQ=255)
                                                   ttGTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGccc                                                                                                                                                                                                                      >  1:514222/1‑225 (MQ=255)
                                                   ttGTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACcgctcg                                                                                                                                                                                            >  1:901069/1‑251 (MQ=255)
                                                     gTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGt                                                                                                                                                                                                                 >  1:720418/1‑228 (MQ=255)
                                                     gTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACcgctcg                                                                                                                                                                                            >  1:496612/1‑249 (MQ=255)
                                                     gTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACcgctcg                                                                                                                                                                                            >  1:627316/1‑249 (MQ=255)
                                                           aTGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTc                                                                                                                                                                                     <  1:97340/250‑1 (MQ=255)
                                                           aTGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTc                                                                                                                                                                                     >  1:356781/1‑250 (MQ=255)
                                                             gCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTc                                                                                                                                                                                     >  1:767058/1‑248 (MQ=255)
                                                                    cTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGcc                                                                                                                                                                           <  1:763976/251‑1 (MQ=255)
                                                                             tAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCa                                                                                                                                                                  >  1:717637/1‑251 (MQ=255)
                                                                                cGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAaat                                                                                                                                                               <  1:495097/251‑1 (MQ=255)
                                                                                  gAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGt                                                                                                                                                                                                                 <  1:707424/199‑1 (MQ=255)
                                                                                      cATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTttgat                                                                                                                                                                                                                                    >  1:772620/1‑176 (MQ=255)
                                                                                                      aataaATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCa                                                                                                                                                                                                                    >  1:525698/1‑176 (MQ=255)
                                                                                                      aataaATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAaat                                                                                                                                                               >  1:182624/1‑229 (MQ=255)
                                                                                                              ccGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTg                                                                                                                                 >  1:336893/1‑251 (MQ=255)
                                                                                                              ccGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTg                                                                                                                                 >  1:872091/1‑251 (MQ=255)
                                                                                                                gTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGcg                                                                                                                               >  1:447430/1‑251 (MQ=255)
                                                                                                                 ttAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGcga                                                                                                                              <  1:47459/251‑1 (MQ=255)
                                                                                                                               tcatTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGatt                                                                                                                <  1:296960/251‑1 (MQ=255)
                                                                                                                               tcatTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTCCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGatt                                                                                                                <  1:610014/251‑1 (MQ=255)
                                                                                                                                  tttttACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAaat                                                                                                                                                               >  1:66885/1‑201 (MQ=255)
                                                                                                                                    tttACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGa                                                                                                           <  1:900097/251‑1 (MQ=255)
                                                                                                                                    tttACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGa                                                                                                           <  1:801079/251‑1 (MQ=255)
                                                                                                                                    tttACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGa                                                                                                           <  1:97766/251‑1 (MQ=255)
                                                                                                                                    tttACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGa                                                                                                           <  1:465802/251‑1 (MQ=255)
                                                                                                                                    tttACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGa                                                                                                           <  1:472479/251‑1 (MQ=255)
                                                                                                                                                    aTCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGccc                                                                                                                                                                                                                      <  1:474423/128‑1 (MQ=255)
                                                                                                                                                    aTCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTAc                                                                                                                                                                                                  <  1:406092/148‑1 (MQ=255)
                                                                                                                                                    aTCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTgg                                                                                           >  1:768093/1‑251 (MQ=255)
                                                                                                                                                       atatCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACg                                                                                              >  1:653890/1‑245 (MQ=255)
                                                                                                                                                         atCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCaaa                                                                                       >  1:460016/1‑250 (MQ=255)
                                                                                                                                                               cGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGCTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGtt                                                                                >  1:379709/1‑251 (MQ=255)
                                                                                                                                                                 ccGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGa                                                                              >  1:211686/1‑251 (MQ=255)
                                                                                                                                                                         gTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAAtt                                                                      >  1:525463/1‑251 (MQ=255)
                                                                                                                                                                          tGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAAt                                                                       <  1:658363/249‑1 (MQ=255)
                                                                                                                                                                             cAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGa                                                                                                <  1:378928/221‑1 (MQ=255)
                                                                                                                                                                               gTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTc                                                                                                                                                                                     >  1:821819/1‑134 (MQ=255)
                                                                                                                                                                               gTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGa                                                                                                           >  1:221020/1‑208 (MQ=255)
                                                                                                                                                                               gTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATcgcg                                                                >  1:83768/1‑251 (MQ=255)
                                                                                                                                                                               gTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATcgcg                                                                >  1:569353/1‑251 (MQ=255)
                                                                                                                                                                                tGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTAt                                                                    <  1:674206/246‑1 (MQ=255)
                                                                                                                                                                                    aCTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGa                                                                                                                                                                                                              >  1:688094/1‑104 (MQ=255)
                                                                                                                                                                                       ttGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATCGCGCCCTTTAt                                                        <  1:68839/251‑1 (MQ=255)
                                                                                                                                                                                         gCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGATGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATCGCGCCCTTTATCg                                                      >  1:856212/1‑251 (MQ=255)
                                                                                                                                                                                          cTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATCGCGCCCTTTATCGa                                                     >  1:38674/1‑251 (MQ=255)
                                                                                                                                                                                                             gattgattTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATCGCGCCCTTTATCGACTTCAGCATTATCATCAAg                                  >  1:757623/1‑251 (MQ=255)
                                                                                                                                                                                                               ttgattTTGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATCGCGCCCTTTATCGACTTCAGCATTATCATCAAGat                                <  1:122499/251‑1 (MQ=255)
                                                                                                                                                                                                                    tttGTTGCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATCGCGCCCTTTATCGACTTCAGCATTATCATCAAGatcaatc                           <  1:752578/251‑1 (MQ=255)
                                                                                                                                                                                                                     ttgttgCTCAAGAAGTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACAACGTTGACGTGAATTATCGCGCCCTTTATCGACTTCAGCATTATCATCAAGatcaatca                          <  1:74241/251‑1 (MQ=255)
                                                                                                                                                                                                                              aagaagTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATCGCGCCCTTTATCGACTTCAGCATTATCATCAAGATCAATCACGTCCGTta                 <  1:374850/251‑1 (MQ=255)
                                                                                                                                                                                                                                gaagTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATCGCGCCCTTTATCGAc                                                    >  1:265038/1‑214 (MQ=255)
                                                                                                                                                                                                                                   gTTCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGACTTATCGCGCCCTTTATCGACTTCAGCATTATCATCAAGTTCAATCACGTCCTTTACTACa            >  1:388837/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ttCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATCGCGCCCTTTATCGACTTCAGCATTATCATCAAGATCAATCACGTCCGTtactac             <  1:267507/249‑1 (MQ=255)
                                                                                                                                                                                                                                    ttCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATCGCGCCCTTTATCGACTTCAGCATTATCATCAAGATCAATCACGTCCGTtactac             <  1:228711/249‑1 (MQ=255)
                                                                                                                                                                                                                                    ttCATCAACAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATCGCGCCCTTTATCGACTTCAGCATTATCATCAAGATCAATCACGTCCGTtactac             <  1:541334/249‑1 (MQ=255)
                                                                                                                                                                                                                                           aCAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATCGCGCCCTTTATCGACTTCAGCATTATCATCAAGat                                <  1:654528/223‑1 (MQ=255)
                                                                                                                                                                                                                                           aCAGATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATCGCGCCCTTTATCGACTTCAGCATTATCATCAAGATCAATCACGTCCGTta                 >  1:532184/1‑238 (MQ=255)
                                                                                                                                                                                                                                              gATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGa                                                                                                           >  1:220215/1‑145 (MQ=255)
                                                                                                                                                                                                                                              gATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAAt                                                                       <  1:747540/181‑1 (MQ=255)
                                                                                                                                                                                                                                              gATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATCGCGCCCTTTATCGACTTCAGCATTATCATCAAGa                                 <  1:129580/219‑1 (MQ=255)
                                                                                                                                                                                                                                              gATAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATCGCGCCCTTTATCGACTTCAGCATTATCATCAAGATCAATCACGTCCGTTACTACAATCTcttctt  >  1:784171/1‑250 (MQ=255)
                                                                                                                                                                                                                                               aTAAGAAGCTATTTTCTTTTGATTGAACGGACTGCCCTATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATCGCGCCCTTTATCGACTTCAGCAt                                            <  1:445598/207‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
CACTTCAGCTGTCGCTTCCATTGGAGCGGATACCACATTGTTATTTAAGCGTTGTAACGATGCTTGTACTTTTTTCTTAACGGAATCATTTGCCGGTTGGATAATAAATTCCGTTAATGGCTGAACATCATTTTTACTTTTTTCAACAATCATATCCAACGCCGGCATAGTGACAGTGATACTTTGCTCCGGTTCGCCTGCATCAGATTGATTTTGTTGCTCAAGAAGTTCATCAACAGATAAGAGGCTATTTTCTTTTGATTGAACGGACTGCCCCATGTCGATCCACACTTTACCGCTCGCTAATTCCGGTGAAGCCACTGCCGCAAATAATGGCATCGGCGATTTTTCTTGCGTATTGCGACGACGGCGCTGATTATTGACGCGCAAATGACGTGGCAAACGACGTTGACGTGAATTATCGCGCCCTTTATCGACTTCAGCATTATCATCAAGATCAATCACGTCCGTTACTACAATCTCTTCTT  >  1/607688‑608175

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