Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 3,845,705 +G intergenic (‑72/+26) PP_3394 ← / ← PP_3395 putative 3‑hydroxy‑3‑methylglutaryl‑CoA lyase/Transcriptional regulator, LysR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154513,845,7031.G94.1% 106.3 / ‑0.6 34intergenic (‑70/+28)PP_3394/PP_3395putative 3‑hydroxy‑3‑methylglutaryl‑CoA lyase/Transcriptional regulator, LysR family
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (15/17);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CGGTTACTTGCAGGTCCGGGTACGTGAGGGCGTGAGCCACCACTTCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑C‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCC  >  AE015451/3845408‑3845995
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                    
gggTTACTTGCAGGTCCGGGTACGTGAGGGCGTGAGCCACCACTTCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑C‑gggct                                                                                                                                                                                                                                                                                                 <  1:468884/300‑4 (MQ=255)
    tACTTGCAGGTCCGGGTACGTGAGGGCGTGAGCCACCACTTCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATc                                                                                                                                                                                                                                                                                              <  1:285660/301‑1 (MQ=255)
          gcgggccGGGGACGGGAGAGCGGGAGCGACCACTTCGCGGGCATCGGCCAGATGCGGCAGCAAGCGTGCGGGCACGAACGACGCGACTGCCAGATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGAATGATCACTTGCAGGCGCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGt                                                                                                                                                                                                                                                                                        <  2:581628/296‑1 (MQ=255)
          cAGGTCCGGGTACGTGAGGGCGTGAGCCACCACTTCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GAG‑GGCTGATCCATCgc                                                                                                                                                                                                                                                                                        >  1:512211/1‑300 (MQ=255)
             gTCCGGGTACGTGAGGGCGTGAGCCACCACTTCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAg                                                                                                                                                                                                                                                                                     <  1:825283/301‑1 (MQ=255)
                   gTACGTGAGGGCGTGAGCCACCACTTCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATCATCACTTGCCGGCTCTGCAAGCCATCGCGCCGGCCGACCTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGCTTC‑G‑CGGGCTGATCCCTCGTCAGCCACGa                                                                                                                                                                                                                                                                               >  2:629374/1‑301 (MQ=255)
                                                      cGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCACGCGCCCTGCCAGGCCGGTGCACgc                                                                                                                                                                                                                                            >  1:395669/1‑301 (MQ=255)
                                                       ggCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACgcc                                                                                                                                                                                                                                           >  1:369658/1‑300 (MQ=255)
                                                        gCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGTCGGTGCACgcag                                                                                                                                                                                                                                          >  1:268062/1‑301 (MQ=255)
                                                                      caAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCCCCCTGCACGGCCGGTGccagccaccaccggtggtaaa                                                                                                                                                                                                                            >  2:183496/1‑281 (MQ=255)
                                                                           gTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCAGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGAGCAGCCGGCGGCATACCAGGCCGGTGGCGGCAGCAACAGGTGGGCAAGCAgg                                                                                                                                                                                                                       >  1:672796/1‑301 (MQ=255)
                                                                                  aacgaacgaCGCGACGTCCCTAGGGCGCACGCCGGAGCCGGAAGCGGCGTTGCTCCGGCGGCGTTTGGCTGCTGTCGGGATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTGTGTCGGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCAAGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgag                                                                                                                                                                                                                <  2:714777/300‑1 (MQ=255)
                                                                                            gcgACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCga                                                                                                                                                                                                      <  1:870187/301‑1 (MQ=255)
                                                                                             cgacttccattggggtcccGCCGGGGCCGTGAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCgag                                                                                                                                                                                                     <  2:395542/283‑1 (MQ=255)
                                                                                                     atatGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGtgtg                                                                                                                                                                                             <  1:867700/301‑1 (MQ=255)
                                                                                                                  gtgcgcCGTAAGCGGCGTCGCTCAAGCGGCTTTTGGCTGCTGTCGGCCTGATCCCTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATg                                                                                                                                                                                <  1:489381/299‑1 (MQ=255)
                                                                                                                  gggCGCCGTAAGCGGCGTCGATACAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCTATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATg                                                                                                                                                                                <  1:33522/300‑1 (MQ=255)
                                                                                                                       ccGTAGGCGGCGTGGATCCAGCGGCGTTTGGCTGCTGTAGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCGGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGAt                                                                                                                                                                           <  1:629597/301‑1 (MQ=255)
                                                                                                                                         cAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAgcgc                                                                                                                                                         >  1:135576/1‑301 (MQ=255)
                                                                                                                                               ggTTTGGTTGCTGTCGTCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAgtc                                                                                                                                                   <  2:369545/300‑1 (MQ=255)
                                                                                                                                                        tctttCGGCATGATCACTTGCAGGCTCTGCAAGCCATCTCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGtcagtcag                                                                                                                                          <  2:512036/297‑1 (MQ=255)
                                                                                                                                                           gTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGGCGCGCATGCAGCTGACGAATTCACCAAGCGCCCAGTCGGCAGTCAgcgg                                                                                                                                       >  1:856767/1‑300 (MQ=255)
                                                                                                                                                                                      aaGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGGGCGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCCATCGTCATGCAGCGCAATGACCCGCAAGCCAAACAAcaccgc                                                                                                            >  2:140380/1‑301 (MQ=255)
                                                                                                                                                                                         ccATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGAGTACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCTCAATGACCCGCAGGTCGAACAACGCTGCATa                                                                                                         >  1:261611/1‑301 (MQ=255)
                                                                                                                                                                                                gcacgCCGAGTTCATTGATTTTTATTGGTGTCATTCG‑TGTCTCTCTGTCTGAAGCGTGGCGCCAAGCTAAGAGCTGGATGACGGGCAGAAAAGCCGTGTTTC‑G‑CGGGCTGATCCATCGGCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGac                                                                                                  <  2:322651/297‑1 (MQ=255)
                                                                                                                                                                                                          ttCATTGATTGTTAGTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCTCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCtt                                                                                        <  2:315645/301‑1 (MQ=255)
                                                                                                                                                                                                                     ttattgttGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCACTGACCCGCAGGTCGAACAACGCCGCATAGCGGGaccccggctctcgcgggcacccc                                                                             >  2:102652/1‑280 (MQ=255)
                                                                                                                                                                                                                           tgggcA‑TCGGGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCagg                                                                       <  2:135536/297‑1 (MQ=255)
                                                                                                                                                                                                                                                  tgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTAGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGAATAGCGGGACACGGCCTCTTGCGGGAACACCGACAGGCCAAGCCCAGACTGCACAAGAc                                                 >  2:305667/1‑301 (MQ=255)
                                                                                                                                                                                                                                                   gagtcGTTGGGCCAGACTAAGAGCTGGATGCCGGGCAGAAAAGACGGGTTTCGG‑G‑GGCTGATCCATAGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACg                                                <  2:787631/297‑1 (MQ=255)
                                                                                                                                                                                                                                                         ttGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGTCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGc                                          <  2:267972/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                      agagCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACGGCACCAAACGGTAGTCCGCGTcgaccgac                             >  2:774230/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                          cTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACCACGCCGCATAGCGGGACACGGCCTCCTGCGGGAACGCCGCCAGGCCAAGGCCAAACTGCACCAGACGGGATGCCGCGTCCAACGAGGaca                         >  2:416411/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                           aaggcgggtTTC‑G‑CGGGCTGATCCATCGTGAGCAACGAGGGAGCCGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGc        <  1:140429/294‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                 tgtTTC‑G‑CGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACACGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCTAAAATGCACCAGCCGGTATGCCGCGTCGAACGACGACACGTTcgaacgaagatcggaagag  >  1:799793/1‑290 (MQ=255)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                    
CGGTTACTTGCAGGTCCGGGTACGTGAGGGCGTGAGCCACCACTTCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCA‑TCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑C‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCC  >  AE015451/3845408‑3845995

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