Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 2,958,523 C→T A428V (GCA→GTA)  PP_2589 → Aldehyde dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154512,958,5230CT98.0% 173.3 / NA 50A428V (GCA→GTA) PP_2589Aldehyde dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (23/26);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATCCGGCGAGCCGCGCCGGGGCCATTGTCGATGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGCAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTG  >  AE015451/2958239‑2958815
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                    
ggATCCGGCGAGCCGCGCCGGGGCCATTGTCGATGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAgg                                                                                                                                                                                                                                                                                      >  3:375121/1‑300 (MQ=255)
        cGAGCCGCGCCGGGGCCATTGTCGATGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGCCtt                                                                                                                                                                                                                                                                              >  4:283260/1‑301 (MQ=255)
           ccacgacgctcttccgatctgATGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGTCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCaa                                                                                                                                                                                                                                                                           <  1:512028/281‑1 (MQ=255)
                  gggggCCATTGTCGATGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTg                                                                                                                                                                                                                                                                    <  1:375121/300‑1 (MQ=255)
                  cGGGGCCATTGTCGATGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTg                                                                                                                                                                                                                                                                    >  3:232596/1‑301 (MQ=255)
                         aTTGTCGATGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCggg                                                                                                                                                                                                                                                             <  1:582916/301‑1 (MQ=255)
                               gATGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGTCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCaaagatcggaagagcacacgtc                                                                                                                                                                                                                                                       >  3:512030/1‑281 (MQ=255)
                                aTGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGcgcg                                                                                                                                                                                                                                                      >  1:95115/1‑301 (MQ=255)
                                aTGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGAGCCTTCGACACGGGAGAAGAAGCCGTGCGCCTGGCCAACGAACGCATCGATGGACTGGCAGCCGCGGTGTGGgacgacgacctaacaccggtgtaccgggtccgccg                                                                                                                                                                                                                                                      >  2:343380/1‑275 (MQ=255)
                                   ccGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGc                                                                                                                                                                                                                                                   <  3:49859/301‑1 (MQ=255)
                                        cAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCCCGGGCGTTgg                                                                                                                                                                                                                                              >  1:207061/1‑300 (MQ=255)
                                          gACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGGCTCGGGGGGGAGCGAAGATCTCAACCAGGGGCACCGGGTAGGGCGGGCGTTGaaa                                                                                                                                                                                                                                            >  1:350428/1‑301 (MQ=255)
                                            tccgatctaTCGAGGCCGGCATTGGGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGGGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGc                                                                                                                                                                                                                                          <  1:526450/293‑1 (MQ=255)
                                                 cGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACAACGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTa                                                                                                                                                                                                                                     <  3:154499/301‑1 (MQ=255)
                                                    aTCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGcagatcgga                                                                                                                                                                                                                                  >  3:526452/1‑293 (MQ=255)
                                                      cttccgatctaTTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTa                                                                                                                                                                                                                                <  4:672169/291‑1 (MQ=255)
                                                       gAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTAt                                                                                                                                                                                                                               <  3:207059/301‑1 (MQ=255)
                                                                aTTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTaagatcggaag                                                                                                                                                                                                                      >  2:671947/1‑291 (MQ=255)
                                                                        tgcCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGAt                                                                                                                                                                                                              >  3:95996/1‑301 (MQ=255)
                                                                                 agagGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGAt                                                                                                                                                                                                     >  4:478955/1‑301 (MQ=255)
                                                                                            gATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGcc                                                                                                                                                                                          <  1:232598/301‑1 (MQ=255)
                                                                                             aTTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTTACGGTGCCg                                                                                                                                                                                         >  1:267045/1‑301 (MQ=255)
                                                                                                ggTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTc                                                                                                                                                                                      <  3:350429/301‑1 (MQ=255)
                                                                                                               aCGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAg                                                                                                                                                                       >  3:465189/1‑301 (MQ=255)
                                                                                                                                 ttCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGAt                                                                                                                                                     <  3:267044/301‑1 (MQ=255)
                                                                                                                                 gttcagacgtgtgctcttccgatctaTTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGAt                                                                                                                                                     <  4:531647/276‑1 (MQ=255)
                                                                                                                                  tCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGAAGGTGCTGTTCGGTGGGGGCCAGCAGGCgggggtggggcgggatg                                                                                                                                                    >  2:64686/1‑288 (MQ=255)
                                                                                                                                   gggaaaagtaaaTCGGGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCt                                                                                                                                                   <  2:478821/290‑1 (MQ=255)
                                                                                                                                       ccttcccctacacgccgctcttccgatcGGTGGCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCa                                                                                                                                               <  2:603959/274‑1 (MQ=255)
                                                                                                                                                cGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCAttcg                                                                                                                                      <  4:485815/301‑1 (MQ=255)
                                                                                                                                                 gAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCAttcgt                                                                                                                                     <  3:95118/301‑1 (MQ=255)
                                                                                                                                                          aTTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGGTGGGGTTTGGCcgctatagagcggaggaggggcgtggtgggg                                                                                                                            >  2:531497/1‑273 (MQ=255)
                                                                                                                                                             gatctcGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTAtt                                                                                                                         <  4:167241/296‑1 (MQ=255)
                                                                                                                                                                  cGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGttttagatc                                                                                                                    >  2:167199/1‑293 (MQ=255)
                                                                                                                                                                    gtgtCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCaagatcggaagagcacacgtctgaactcca                                                                                                                  >  4:604156/1‑272 (MQ=255)
                                                                                                                                                                    gtgtCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAAcc                                                                                                                  >  1:543349/1‑300 (MQ=255)
                                                                                                                                                                                     ggCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCa                                                                                                <  4:343464/300‑1 (MQ=255)
                                                                                                                                                                                            cgcgcgAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTgcg                                                                                          <  4:64705/301‑1 (MQ=255)
                                                                                                                                                                                                           ccGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGa                                                                           <  4:322823/300‑1 (MQ=255)
                                                                                                                                                                                                               cccGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGTGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTAtttgacat                                                                       >  2:727543/1‑295 (MQ=255)
                                                                                                                                                                                                                  cctcttccgatctaGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGa                                                                    <  2:688529/288‑1 (MQ=255)
                                                                                                                                                                                                                              cAGCGCCGTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTg                                                        <  1:95993/301‑1 (MQ=255)
                                                                                                                                                                                                                               aGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGGTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACAACCTGGTACCAGTTGCGCGGCCGATGTATTGATCATTGaagatcggaagagc                                                       >  4:688759/1‑288 (MQ=255)
                                                                                                                                                                                                                                  gcCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTc                                                    <  1:465187/301‑1 (MQ=255)
                                                                                                                                                                                                                                  gcCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTc                                                    <  3:543351/301‑1 (MQ=255)
                                                                                                                                                                                                                                          ccgctttcccgttctctgctccccggCCAACGTCCGCATCTATGACCTGGTGTCCTGTGTGGGGATCGACGACTCCAACCAGGTGCACCGTGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCAt                                            <  2:263394/277‑1 (MQ=255)
                                                                                                                                                                                                                                              gaagaagaaGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGTGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGccc                                        <  4:727803/301‑1 (MQ=255)
                                                                                                                                                                                                                                                         cGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACAGCGTCGTCTGCTTCGCGGGCAtagatcggaagagcac                             >  4:263460/1‑286 (MQ=255)
                                                                                                                                                                                                                                                                              cATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAGGTCTGAGGGtt        >  3:136376/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                    tGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCTCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTg  >  2:386/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                    
GGATCCGGCGAGCCGCGCCGGGGCCATTGTCGATGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGCAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTG  >  AE015451/2958239‑2958815

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: