Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 2,837,071 Δ6 bp intergenic (‑23/+155) sad‑I ← / ← PP_2489 NAD+‑dependent succinate semialdehyde dehydrogenase/putative Xenobiotic reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154512,837,0710G.100.0% 14.8 / NA 5intergenic (‑23/+160)sad‑I/PP_2489NAD+‑dependent succinate semialdehyde dehydrogenase/putative Xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (1/4);  total (1/4)
*AE0154512,837,0720G.100.0% 14.8 / NA 5intergenic (‑24/+159)sad‑I/PP_2489NAD+‑dependent succinate semialdehyde dehydrogenase/putative Xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (1/4);  total (1/4)
*AE0154512,837,0730T.100.0% 14.8 / NA 5intergenic (‑25/+158)sad‑I/PP_2489NAD+‑dependent succinate semialdehyde dehydrogenase/putative Xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (1/4);  total (1/4)
*AE0154512,837,0740G.100.0% 14.8 / NA 5intergenic (‑26/+157)sad‑I/PP_2489NAD+‑dependent succinate semialdehyde dehydrogenase/putative Xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (1/4);  total (1/4)
*AE0154512,837,0750G.100.0% 14.8 / NA 5intergenic (‑27/+156)sad‑I/PP_2489NAD+‑dependent succinate semialdehyde dehydrogenase/putative Xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (1/4);  total (1/4)
*AE0154512,837,0760C.80.0% 10.0 / ‑1.9 5intergenic (‑28/+155)sad‑I/PP_2489NAD+‑dependent succinate semialdehyde dehydrogenase/putative Xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (1/3);  total (1/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCGCGCTGAAATCTCGACTTCGCCCAAAGCTTCAGCGGTGATCTTGCCCATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGAATGGTGGCATGAAAAGCTCCGGTCTTTCCGTACGTGGTGACAGGGGCTGAAACTCGGATTACTCCACCGGCGTATACCGTTCGAGCCAGTGGGCGTAAG  >  AE015451/2836780‑2837271
                                                                                                                                                                                                                                                                                                   ||||||                                                                                                                                                                                                   
aTTCGCGCTGAAATCTCGACTTCGCCCAAAGCTTCAGCGGTGATCTTGCCCATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgttgcatggtggctggg                                                                                                                                                                                                 >  3:575493/1‑297 (MQ=255)
aTTCGCGCTGAAATCTCGACTTCGCCCAAAGCTTCAGCGGTGATCTTGCCCATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatgg                                                                                                                                                                                                 >  3:124596/1‑301 (MQ=255)
    gcgcTGAAATCTCGACTTCGCCCAAAGCTTCAGCGGTGATCTTGCCCATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggg                                                                                                                                                                                             >  1:56328/1‑300 (MQ=255)
       cTGAAATCTCGACTTCGCCCAAAGCTTCAGCGGTGATCTTGCCCATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGTCGTAAGCCAGGTGAGTCTCCTACCCgtggcatggtgggatggtggcatg                                                                                                                                                                                          >  2:358696/1‑301 (MQ=255)
         aaaaTCTCGACTTCGCCCAAAGCTTCAGCGGTGATCTTGCCCATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggcatggt                                                                                                                                                                                        <  1:579080/300‑1 (MQ=255)
            aTCTCGACTTCGCCCAAAGCTTCAGCGGTGATCTTGCCCATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggcatggtggg                                                                                                                                                                                     >  1:579211/1‑300 (MQ=255)
            aTCTCGACTTCGCCCAAAGCTTCAGCGGTGATCTTGCCCATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtgacatggtgggagggtggg                                                                                                                                                                                     >  1:640623/1‑300 (MQ=255)
                 gACTTCGCCCAAAGCTTCAGCGGTGATCTTGCCCATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggcatggtggcatggt                                                                                                                                                                                <  4:224804/301‑1 (MQ=255)
                         ccAAAGCTTCAGCGGTGATCTTGCCCATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggcatggtggcatggtggcatggg                                                                                                                                                                        >  1:101177/1‑300 (MQ=255)
                            ttctttCAGCGGTGATCTTGCCCATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggcatggtggcatggtggcatggtggc                                                                                                                                                                     <  3:114737/297‑1 (MQ=255)
                             aGCTTCAGCGGTGATCTTGCCCATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatgtgggcatggggcaatgtgggcttggggtactgtggctg                                                                                                                                                                    >  2:503053/1‑272 (MQ=255)
                              ccTTCAGCGGTGATCTTGCCCATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggcatggtggcatggtggcatggtggcat                                                                                                                                                                   <  1:183626/300‑1 (MQ=255)
                                ttCAGCGGTGATCTTGCCCATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggcatggtggcatggtggcatggtggcatgg                                                                                                                                                                 >  3:56303/1‑301 (MQ=255)
                                ttCAGCGGTGATCTTGCCCATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggcatggtggcatggtggcatggtggcagat                                                                                                                                                                 >  1:114734/1‑298 (MQ=255)
                                              tGCCCATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcat                                                                                                                                                   >  3:636535/1‑301 (MQ=255)
                                               gCCCATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatg                                                                                                                                                  <  4:719887/301‑1 (MQ=255)
                                                  cATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggttgcatggtggcatggtg                                                                                                                                               >  3:208324/1‑301 (MQ=255)
                                                      tctAATGTTTGCAGCTGGGCAGAGTGTCTGGCCTGGTGGCGTAGCAGCTGAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTGTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcat                                                                                                                                           <  1:56302/301‑1 (MQ=255)
                                                           ttccgatctgTGGGCATAGTCGCTGGCCTGGGCGTGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtgg                                                                                                                                      <  4:335130/292‑1 (MQ=255)
                                                            gTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggc                                                                                                                                     <  3:101179/301‑1 (MQ=255)
                                                                    gTGGGCATAGTCGCTGGCCTGGGCGTGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggcatggtggcatggtggcctggtggcaaggtggcagggtggcatggtggcatggtggagaggggag                                                                                                                             >  2:335066/1‑292 (MQ=255)
                                                                          atagtccttGGCCTGGGCGTGTAGCAGCTCAGCTGCCCGCTGGCGAATAGCACTGCCTTCGGTAAGGCTGCGTTCTTTCCAGTCGAGGAACGCGCTGTGTGCGTCGTCCAGCGCGTGATTGACCTGGTGGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggcatggtggaatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtg                                                                                                                       <  3:640637/293‑1 (MQ=255)
                                                                                       tGGGCGCGTAGCAGCTCAGCTGCCCGCTGCAGAATGGCACTGCGTTCGGTAAAGCTGCGTGCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGGCTGATTGGCCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatgatggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggaatggtggcatg                                                                                                          <  4:525179/301‑1 (MQ=255)
                                                                                                      tcagctGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtgccgtggtagcatggtggcatgttggcatggttgaatggtggcatgagatcggaagagcac                                                                                           >  3:75490/1‑286 (MQ=255)
                                                                                                        cgctGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCgtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggcatggtggaatggtggcatgaaaagctccggtctttc                                                                                         <  3:56329/300‑18 (MQ=255)
                                                                                                                                                cGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGAATGGTGGCATGAAAAGCTCCGGTCTTTCCGTACGTGGTGACAGGGGCTGAAACTCGGATTACTCCAcc                                 <  4:516511/301‑1 (MQ=255)
                                                                                                                                                 gTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGAATGGTGGCATGAAAAGCTCCGGTCTTTCCGTACGTGGTGACAGGGGCTGAAACTCGGGTTACTCCACCg                                >  4:466698/1‑301 (MQ=255)
                                                                                                                                                         cAGTCGAGGAGCTCGGTATGTGCGTCGTCGAGCGCCTTATGGGCCTGGTCGTCGGTGGTTTCTGGGAAGGGCTTGAGCAATTGCCCGTTCTAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCC‑‑‑‑‑‑‑‑GAT‑‑‑‑‑‑‑TCGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGAATGGTGGCATGAAAAGCTCCGGTCTTTCCGTACGTGGTGACAGGGGCTGAAACTCGGATTACTCCACCGGCGTATAc                        <  4:336208/301‑1 (MQ=255)
                                                                                                                                                          aGTCGAGGAACGCTCTATGTGCGGCTTCCAACGTCCGATTGGCCTGGTCGGCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATCGGGGTTGATAGTGGCGTAAGCCATGTGATGCTCCTATTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGAATGGTGGCATGAAAAGCTCCGGTCTTTCCGTACGTGGTGACAGGGGCTGAAACTCGGATTACTCCACCGGCGTATAcc                       <  3:579226/301‑1 (MQ=255)
                                                                                                                                                                               cggcggccAGCGGCTGAGTGACCTGGACGTCTGTGGTTTCTGGAAAGGTATTGCGCAATTGCCCGTTATAGGGGTTGATAGTTGCGTAAGCCATTGGATTCTCCTACT‑‑‑‑‑‑‑‑‑‑‑‑‑CA‑‑‑TGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGTTGGCATGGTGGCATGGTGGCATGGTGGAATGGTGGCATGAAAAGCTCCGGTCTTTCCGTACGTGGTGACAGGGGCTGAAACTCGGATTACTCCACCGGCGTATACCGTTCGAGCCAGTGGGCGTAAg  <  4:358769/295‑1 (MQ=21)
                                                                                                                                                                                                                                                                                                   ||||||                                                                                                                                                                                                   
ATTCGCGCTGAAATCTCGACTTCGCCCAAAGCTTCAGCGGTGATCTTGCCCATCTCTAAGGTTAGCAGGTGGGCATAGTCGCTGGCCTGGGCGCGTAGCAGCTCAGCTGCCCGCTGGAGAATGGCACTGCGTTCGGTAAAGCTGCGTTCTTTCCAGTCGAGGAACGCGCTATGTGCGTCGTCCAGCGCCTGATTGACCTGGTCGTCGGTGGTTTCTGGGAAGGTCTTGAGCAATTCCCCGTTATAGGGGTTGATAGTGGCGTAAGCCATGTGATTCTCCTACTCGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGCATGGTGGAATGGTGGCATGAAAAGCTCCGGTCTTTCCGTACGTGGTGACAGGGGCTGAAACTCGGATTACTCCACCGGCGTATACCGTTCGAGCCAGTGGGCGTAAG  >  AE015451/2836780‑2837271

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: