Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 5,300,963 T→G 5.3% D937A (GAT→GCT)  recB ← Chi activated ATP‑dependent DNA helicase and dsDNA/ssDNA exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154515,300,9630TG 5.3% 257.7 / 4.0 71D937A (GAT→GCT) recBChi activated ATP‑dependent DNA helicase and dsDNA/ssDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base T (35/32);  new base G (2/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

ATTGCAGCGCTGGCCAACGGTACGCTCTATCAGTTGAGGGTTAGCGCTGACCTCGCTGAAGCCCTCACGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGCCTGGGCAAGTGATGTGCGGGCAACTTGCCGCCAAAGCCTGCAACCAGTCCGTCAGCTGCGCCGAGCCAGCCAAGGCAAGGCCGCCCCC  >  AE015451/5300728‑5301183
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                            
aTTGCAGCGCTGGCCAACGGTACGCTCTATCAGTTGAGGGTTAGCGCTGACCTCGCTGAAGCCCTCACGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGAt                                                                                                                                                                                                                              >  1:1062255/1‑236 (MQ=255)
    cAGCGCTGGCCAACGGTACGCTCTATCAGTTGAGGGTTAGCGCTGACCTCGCTGAAGCCCTCACGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGc                                                                                                                                                                                                                           >  1:687754/1‑235 (MQ=255)
            gCCAACGGTACGCTCTATCAGTTGAGGGTTAGCGCTGACCTCGCTGAAGCCCTCACGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGa                                                                                                                                                                                                                   <  1:1176450/235‑1 (MQ=255)
             ccAACGGTACGCTCTATCAGTTGAGGGTTAGCGCTGACCTCGCTGAAGCCCTCACGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAg                                                                                                                                                                                                                 <  1:1479867/236‑1 (MQ=255)
              cAACGGTACGCTCTATCAGTTGAGGGTTAGCGCTGACCTCGCTGAAGCCCTCACGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAgg                                                                                                                                                                                                                <  1:1326621/236‑1 (MQ=255)
               aaCGGTACGCTCTATCAGTTGAGGGTTAGCGCTGACCTCGCTGAAGCCCTCACGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAg                                                                                                                                                                                                                 <  2:763104/234‑1 (MQ=255)
                         tctATCAGTTGAGGGTTAGCGCTGACCTCGCTGAAGCCCTCACGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGAt                                                                                                                                                                                                                              >  2:1088340/1‑211 (MQ=255)
                               aGTTGAGGGTTAGCGCTGACCTCGCTGAAGCCCTCACGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAgg                                                                                                                                                                                                                <  2:1337184/219‑1 (MQ=255)
                                      ggTTAGCGCTGACCTCGCTGAAGCCCTCACGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAgg                                                                                                                                                                                                                >  1:1127782/1‑212 (MQ=255)
                                      ggTTAGCGCTGACCTCGCTGAAGCCCTCACGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAgg                                                                                                                                                                                                                <  2:1127779/212‑1 (MQ=255)
                                      ggTTAGCGCTGACCTCGCTGAAGCCCTCACGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATcc                                                                                                                                                                                             >  1:710272/1‑231 (MQ=255)
                                        ttAGCGCTGACCTCGCTGAAGCCCTCACGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCCGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACCCGCAGGGCa                                                                                                                                                                                                             >  2:270103/1‑213 (MQ=255)
                                                 aCCTCGCTGAAGCCCTCACGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGc                                                                                                                                                                             <  1:16748/236‑1 (MQ=255)
                                                    tCGCTGAAGCCCTCACGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGATCAGACCCTCGGCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGcc                                                                                                                                                                            <  2:635518/234‑1 (MQ=255)
                                                    tCGCTGAAGCCCGCACGGCCTGCCCACGCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGa                                                                                                                                                                                                                   <  2:124367/195‑1 (MQ=255)
                                                       cTGAAGCCCTCACGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTa                                                                                                                                                                                                        <  1:921329/203‑1 (MQ=255)
                                                                   cGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATg                                                                                                                                                                       >  1:1322156/1‑224 (MQ=255)
                                                                   cGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATg                                                                                                                                                                       <  2:1322153/224‑1 (MQ=255)
                                                                   cGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATg                                                                                                                                                                       >  2:1219087/1‑224 (MQ=255)
                                                                     gCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTg                                                                                                                                                                                      >  2:1170608/1‑207 (MQ=255)
                                                                            cACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCgcg                                                                                                                                                         >  1:286308/1‑227 (MQ=255)
                                                                            cACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTg                                                                                                                                                         >  2:900526/1‑229 (MQ=255)
                                                                                    cAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAgg                                                                                                                                                                                                                <  1:1430447/166‑1 (MQ=255)
                                                                                       aCCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCg                                                                                                                                                                                                                            <  2:747404/151‑1 (MQ=255)
                                                                                         cATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAgg                                                                                                                                                                                                      <  2:65537/171‑1 (MQ=255)
                                                                                          aTGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGtt                                                                                                                                       >  2:965642/1‑233 (MQ=255)
                                                                                          aTGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCgg                                                                                                                                    <  1:913704/236‑1 (MQ=255)
                                                                                            ggaggaAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGACggg                                                                                                                                                                                                               >  2:51513/1‑159 (MQ=255)
                                                                                              aggaAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATg                                                                                                                                                                       >  1:30277/1‑197 (MQ=255)
                                                                                                      gCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCTATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGc                                                                                                                            >  1:721718/1‑232 (MQ=255)
                                                                                                       ccGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATg                                                                                                                         <  2:1054777/234‑1 (MQ=255)
                                                                                                       ccGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTATGATCGCCGACACGAAg                                                                                                                                                                                                                 <  1:542948/146‑1 (MQ=255)
                                                                                                          gggTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATa                                                                                                                    >  1:903836/1‑236 (MQ=255)
                                                                                                           ggTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGa                                                                                                                      >  2:766334/1‑233 (MQ=255)
                                                                                                           ggTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGAGCGCCGACACg                                                                                                                                                                                                                    >  2:1345133/1‑139 (MQ=255)
                                                                                                                  cccccGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGc                                                                                                                                                                                                  >  2:843838/1‑150 (MQ=255)
                                                                                                                  cccccGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATa                                                                                                                    >  2:667610/1‑228 (MQ=255)
                                                                                                                      cgcgGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATc                                                                                                        <  1:691053/236‑1 (MQ=255)
                                                                                                                       gcgGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGtgtg                                                                                                                              <  2:1275528/213‑1 (MQ=255)
                                                                                                                        cgGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCg                                                                                                      >  1:1420391/1‑236 (MQ=255)
                                                                                                                         gggAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGc                                                                                                     >  1:818349/1‑236 (MQ=255)
                                                                                                                                  tGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCt                                                                                                                                                                                    <  2:1371752/148‑1 (MQ=255)
                                                                                                                                    gATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGacac                                                                                                                                                                                                                     >  2:491937/1‑113 (MQ=255)
                                                                                                                                    gATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGt                                                                                                                                                                   <  2:1247201/163‑1 (MQ=255)
                                                                                                                                    gATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGc                                                                                          >  1:1366557/1‑236 (MQ=255)
                                                                                                                                            cgcTGTCGGCCGGGCCCTCGCGCAACGCCTGGGCATCGGGCCCCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCt                                                                                                            <  1:239244/210‑1 (MQ=255)
                                                                                                                                                     ccGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGacac                                                                                                                                                                                                                     <  2:88281/96‑1 (MQ=255)
                                                                                                                                                                 gcAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCGGCGGCAGGCCTGGGCAAGTGATGTGCGGGCAACTTgcc                                                             >  1:258799/1‑236 (MQ=255)
                                                                                                                                                                  cAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATcc                                                                                                       <  1:1514808/193‑1 (MQ=255)
                                                                                                                                                                      acCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCAc                                                                                                                                          <  2:710270/154‑1 (MQ=255)
                                                                                                                                                                      acCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATCGCGTGGTCTGCGTGCTGGCAGCAGCTCACGt                                                                                                                                        >  2:335928/1‑156 (MQ=255)
                                                                                                                                                                        cTGGGCGTCGGGCACGTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGTTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGc                                                                                                                                                                             <  2:93423/117‑1 (MQ=255)
                                                                                                                                                                        cTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTTGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGCCTGGGCAAGTGATGTGCGGGCAACTTGCCGCCaaa                                                       <  1:799139/235‑1 (MQ=255)
                                                                                                                                                                                   gCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATGGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTACTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGc                                                                                                     >  2:799137/1‑178 (MQ=255)
                                                                                                                                                                                   gCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGCCTGGGCAAGTGATGTGCGGGCAACTTGCCGCCAAAGCCTGCAACCAg                                           >  1:1483351/1‑236 (MQ=255)
                                                                                                                                                                                   gCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGATCAGACCCTCGGCGCCGACACGAAGGGCACTGTAGGAAGCGATc                                                                                                                                                                                              <  2:1506677/89‑1 (MQ=16)
                                                                                                                                                                                   gCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGATCAGACCCTCGGCGCCGACACGAAGGGCACTGTAGGAAGCGATc                                                                                                                                                                                              >  1:1506680/1‑89 (MQ=21)
                                                                                                                                                                                         cgtcgtCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGCCTGGGCAAGTGATGTGCGGGCAACTTGCCGCCAAAGCCTGCAACCAGTcc                                        >  2:628349/1‑233 (MQ=255)
                                                                                                                                                                                          gtcgtcCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCAGAGGGTCTGATCGCCGACTCGCAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGCCTGGGCaa                                                                                  >  2:528991/1‑190 (MQ=255)
                                                                                                                                                                                            cgtcCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGCCTGGGCAAGTGATGTGCGGGCAACTTGCCGCCAAAGCCTGCAACCAGTCCGt                                      <  1:1284022/232‑1 (MQ=255)
                                                                                                                                                                                                    ggTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCGGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGCCTGGGCAAGttgcccg                                                                          <  2:1188897/187‑7 (MQ=255)
                                                                                                                                                                                                          tgtgCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGCCTGGGCAAGTGATGTGCGGGCAACTTGCCGCCAAAGCCt                                                   <  1:1027277/205‑1 (MQ=255)
                                                                                                                                                                                                           gtgCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGCCTGGGCAAGTGATGTGCGGGCAACTTGCCGCCAAAGCCTGCAACCAGGCCGTCAGCTGCGCCGagccagcc                   >  1:398519/1‑236 (MQ=255)
                                                                                                                                                                                                              ccTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACa                                                                                                                  <  1:1007625/138‑1 (MQ=255)
                                                                                                                                                                                                              ccTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGCCTGGGCAAGTGATGTGCGGGCAACTTGCCGCCAAAGCCTGCAACCAGTCCGTCAGCTGCGCCGAGCCAGCcaag                >  1:1507091/1‑236 (MQ=255)
                                                                                                                                                                                                               cTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGCCTGGGCAAGTGATGTGCGGGCAACTTGCCGCCAAAGCCTGCAACCAGTCCGTCAGCTGCGCCGAGCCAGCcaagg               >  1:748556/1‑236 (MQ=255)
                                                                                                                                                                                                                  cttcgtGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGCCTGGGCAAGTgc                                                                              >  2:1168697/6‑169 (MQ=255)
                                                                                                                                                                                                                     gCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGCCTGCGCAAGTGATGTGCGGGCAACTTgcc                                                             >  1:1216305/1‑184 (MQ=255)
                                                                                                                                                                                                                           ggCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGCCTGGGCAAGTGATGTGCGGGCAACTTGCCGCCAAAGCCTGCAACCAg                                           >  2:858578/1‑196 (MQ=255)
                                                                                                                                                                                                                           ggCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGCCTGGGCAAGTGATGTGCGGGCAACTTGCCGCCAAAGCCTGCAACCAg                                           <  1:858580/196‑1 (MQ=255)
                                                                                                                                                                                                                             cgcCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGCCTGGGCAAGTGATGTGCGGGCAACTTGCCGCCAAAGCCTGCAACCAGTCCGTCAGCTGCGCCGAGCCAGCCAAGGCAAGGCCGccccc  >  1:1324550/1‑235 (MQ=255)
                                                                                                                                                                                                                               ccGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGCCTGGGCAAGTGATGTGCGGGCAACTTGCCGCCAAAGCCTGCAACCAGTCCGTCAGCTGCGcc                            <  2:1386193/207‑1 (MQ=255)
                                                                                                                                                                                                                                 gAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCGCGGCAGCCGCATGGCGTGGTCGGCGTGCTGGCAGCAGCCCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCg                                                                                                >  1:1130089/1‑137 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                            
ATTGCAGCGCTGGCCAACGGTACGCTCTATCAGTTGAGGGTTAGCGCTGACCTCGCTGAAGCCCTCACGGCCTGCCCACTCCAGCAAACCATGGAGGAAGGTGCCGGGGTTAGGCCCCCGCGGGAAGCGATGGATATCGCCGCTGTCGGCCGGGACCTCGCGCAACACCTGGGCATCGGGCACCTCGTCGTCCAGCAGTTGCTGTGCCTGCGAGCTGTCGGCGCCGAGGGTCTGATCGCCGACACGAAGGGCACTGTAGGAAGCGATCCACCAGTGCTCGGCAGCCGCATGGCGTGGTCTGCGTGCTGGCAGCAGCTCACGTTCGGCGTGTGGCATGCGATACATCTGCTCATCCGCCTGCGGCAGGCCTGGGCAAGTGATGTGCGGGCAACTTGCCGCCAAAGCCTGCAACCAGTCCGTCAGCTGCGCCGAGCCAGCCAAGGCAAGGCCGCCCCC  >  AE015451/5300728‑5301183

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