Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 4,740,819 2 bp→CT intergenic (+56/+696) gltA → / ← yeiW citrate synthase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,740,8190TC97.8% 157.9 / ‑4.1 45intergenic (+56/+697)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (23/21);  minor base A (0/1);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01
*AE0154514,740,8200GT97.8% 153.4 / NA 45intergenic (+57/+696)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (23/21);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTCATCGAGCGCTCGCTGTACCCGAACGTCGACTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCGGCTCTTCTTTCAGGTTAC  >  AE015451/4740531‑4741107
                                                                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                                                                               
ccTTCATCGAGCGCTCGCTGTACCCGAACGTCGACTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctct                                                                                                                                                                                                                                                                                     <  1:403213/300‑1 (MQ=255)
 cctCATCGAGCGCTCGCTGTACCCGAACGTCGACTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCATGTGCCGGCCTCtt                                                                                                                                                                                                                                                                                    <  3:55763/299‑1 (MQ=255)
               cgctGTACCCGAACGTCGACTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCTTCGGGCAGCC‑TTTTTTATTTGTCCTGTTCCGGCCTCTTCGCGTCCTTGcccc                                                                                                                                                                                                                                                                      >  2:706309/1‑300 (MQ=255)
                gctGTACCCGAACGTCGACTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcgc                                                                                                                                                                                                                                                                     >  4:701767/1‑301 (MQ=255)
                   gTACCCGAACGTCGACTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                                                                                                                                                  >  1:289979/1‑301 (MQ=255)
                    tgctcttccgatctcTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAag                                                                                                                                                                                                                                                                 <  4:423987/287‑1 (MQ=255)
                      cccGACCGCCGACTTCTACGC‑GGGAACATTCCTCAAGGCGACCGGCATCCCGACCCGCATGTTCACCGTGCTCTTCCCCCTGGCACGTACCGTGGGCTGGATCTCGAACTGAAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                                                                                                                                                                                               <  1:214793/301‑1 (MQ=255)
                      cccGAACGTCGACTTCTACTCGGGGATCA‑TCCTCAAGGCGCTCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCCGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTTCCGGCCTCTTCGCGGCCTTGccctcgcagcg                                                                                                                                                                                                                                                               >  2:556560/1‑293 (MQ=255)
                      cccGAACGTCGACTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGGGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCTGCCTCTTCGCCGCCTTcccccccaaagg                                                                                                                                                                                                                                                               >  2:80856/1‑293 (MQ=255)
                                  cTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCCCGGCCTTGCCCGCGAagagatcggaagagcg                                                                                                                                                                                                                                                   >  2:423878/1‑289 (MQ=255)
                                                      cacAAGGCGATCGGCAGTCCGACAAGCATGTTCACCGTGATCTGCGCCCCGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGtttt                                                                                                                                                                                                                                <  4:556730/299‑1 (MQ=255)
                                                                  ggCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTACCCTCGGGCAGCC‑TTTTTGATTTTTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCCCGCAAAGGCCCGCCCAGTCCCCACACAACCCCAGTTTTtccccccacccc                                                                                                                                                                                                                    >  2:468386/1‑291 (MQ=255)
                                                                         cGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGTCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTCTCTCCGGCCCCCCACGCAg                                                                                                                                                                                                             >  2:677089/1‑301 (MQ=255)
                                                                           cccAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGcc                                                                                                                                                                                                           <  3:289978/300‑1 (MQ=255)
                                                                            ccAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGGGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGccc                                                                                                                                                                                                          <  4:80875/301‑1 (MQ=255)
                                                                            ccAGCATGTTCACCGTGATCTTCGCCCTGGCACCTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGTGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGccc                                                                                                                                                                                                          <  4:706554/301‑1 (MQ=255)
                                                                                     gacgccgttccgatgtccTGGCACGTACCGGGGGCTGGATCTCGAACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGt                                                                                                                                                                                                 <  2:605685/285‑1 (MQ=255)
                                                                                               cTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCAAAGCCACCACAAAACCTCAGTTTTTCTCCGCCCGCTAACGCAGCCCTTTCAGGGTATTGAACGgc                                                                                                                                                                                       >  3:485874/1‑301 (MQ=255)
                                                                                                     ccTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTagatcggaagagcaca                                                                                                                                                                                 >  4:605882/1‑286 (MQ=255)
                                                                                                      tcggCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcca                                                                                                                                                                                <  4:651469/299‑1 (MQ=255)
                                                                                                         gCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGACCGGCCCATcgcaga                                                                                                                                                                             >  2:651258/1‑296 (MQ=255)
                                                                                                                        gggATCTCGCACTGGAAAGAAATGCTCTCCCGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAAAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTACCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGcca                                                                                                                                                              <  4:468515/300‑1 (MQ=255)
                                                                                                                         ggaacacGCACGGGAAAGAAATGCACTCCCGCCCGTACAAGATCGACCGTCCG‑CGCCAGCTGTAAACCGGCGAACAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAACTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGccac                                                                                                                                                             <  1:485872/295‑1 (MQ=255)
                                                                                                                            tctcGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCAt                                                                                                                                                          <  4:677311/301‑1 (MQ=255)
                                                                                                                                 cagtggaaagaaatctccccccgCCCGTACCACATCGGCCGTCCGACG‑CAGCTGTACACCGGCGAGCAAGAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCg                                                                                                                                                     <  2:701522/279‑1 (MQ=255)
                                                                                                                                   cTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTGGCCGGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGATTTTCTCCGCCAGCGCACGCAGCCCTTCCAGGGTATAGAACGTCGAATACACCACAAACTTGTGGGGCACCATCGCCggg                                                                                                                                                   >  2:43347/1‑300 (MQ=255)
                                                                                                                                        aagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCt                                                                                                                                              <  4:43364/301‑1 (MQ=255)
                                                                                                                                                              cAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCAcc                                                                                                                        >  4:230518/1‑301 (MQ=255)
                                                                                                                                                                     ggCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAgg                                                                                                                 >  3:127148/1‑301 (MQ=255)
                                                                                                                                                                                       gTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCa                                                                                                <  1:127146/301‑1 (MQ=255)
                                                                                                                                                                                                                  cacgacgctcttccgctctaaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCg                                                                     <  2:26278/282‑1 (MQ=255)
                                                                                                                                                                                                                           gAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACtttt                                                            >  4:438773/1‑301 (MQ=255)
                                                                                                                                                                                                                            aaataCCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTg                                                           <  2:230464/297‑1 (MQ=255)
                                                                                                                                                                                                                              aGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTTAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCCCCCTTGGGCACCAGCTTCCAGAAGCCGTCGCCGTTTGcc                                                         >  1:367061/1‑301 (MQ=255)
                                                                                                                                                                                                                                     aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTTGGCCCCAGCTTCCAGAAGCCgagatcggaagagcacacgt                                                  >  4:26288/1‑282 (MQ=255)
                                                                                                                                                                                                                                       ttccgatctagatgCCGAAGGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCAGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACa                                                <  2:189321/288‑1 (MQ=255)
                                                                                                                                                                                                                                                aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCAACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGACCaagatcggaa                                       >  4:189367/1‑293 (MQ=255)
                                                                                                                                                                                                                                                 ggagttcagacgtgtgctcttccgatctggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTtct                                      <  4:720769/273‑1 (MQ=255)
                                                                                                                                                                                                                                                    gCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTTCTCtt                                   >  4:604396/1‑301 (MQ=255)
                                                                                                                                                                                                                                                     ccgatctcAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTTCTCTTc                                  <  4:318582/294‑1 (MQ=255)
                                                                                                                                                                                                                                                     ccGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGTCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTACGTTCTTTCACCCTTGGGCACCAGCTTCCAGACGCCGTCGCCTTTTGCCCCCCGATCCAAAACCTTCTATTc                                  >  2:561859/1‑301 (MQ=255)
                                                                                                                                                                                                                                                            cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCCCCCGAACAAACCCCTTCTCTTcccaccgc                           >  2:318509/1‑295 (MQ=255)
                                                                                                                                                                                                                                                                             ggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTtctagatcggaagagcgtggtgtcggggagg          >  2:720524/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                                      tttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAACCCCTTCTCTTCCGGCAGGTAAGTCGGCTCTTCTTTCAGGTTAc  >  1:53290/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                      tttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCGGCTCTTCTTTCAGGTTAc  <  3:367061/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                                                                               
ACTTCATCGAGCGCTCGCTGTACCCGAACGTCGACTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCGGCTCTTCTTTCAGGTTAC  >  AE015451/4740531‑4741107

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑