Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 2,297,033 C→A 9.3% R1030L (CGC→CTC) ‡ sbcC ← Exonuclease SbcC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154512,297,0330CA 9.3% 167.8 / 11.0 54R1030L (CGC→CTC) ‡sbcCExonuclease SbcC
Reads supporting (aligned to +/- strand):  ref base C (26/23);  new base A (2/3);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-02

GCACCAGCAGGCCCAGGGCGCTGCCGCCGCGCTTGAGGCGGTAGCGACGGGCCAGCTGACGTAGCTGGGCGTTGGCGTGGTGCAGCAGCAGGTCAAGGTTGTAGCCTTGGGCGATCTTGCGGAACACATCACCGGAGGCCGAACCGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAA  >  AE015451/2296803‑2297239
                                                                                                                                                                                                                                      |                                                                                                                                                                                                              
gcaccagcaGGCCCAGGGCGCTGCCGCCGCGCTTGAGGCGGTAGCGACGGGCCAGCTGACGTAGCTGGGCGTTGGCGTGGTGCAGCAGCAGGTCAAGGTTGTAGCCTTGGGCGATCTTGCGGAACTCATCACCGGAGGCCGAACCGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGgcggcg                                                                                                                                                                                                            >  1:517318/1‑235 (MQ=255)
         ggCCCAGGGCGCTGCCGCCGCGCTTGAGGCGGTAGCGACGGGCCAGCTGACGTAGCTGGGCGTTGGCGTGGTGCAGCAGCAGGTCAAGGTTGTAGCCTTGGGCGATCTTGCGGAACACATCACCGGAGGCCGAACCGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTggc                                                                                                                                                                                                                <  1:202494/222‑1 (MQ=255)
              aGGGCGCTGCCGCCGCGCTTGAGGCGGTAGCGACGGGCCAGCTGACGTAGCTGGGCGTTGGCGTGGTGCAGCAGCAGGTCAAGGTTGTAGCCTTGGGCGATATTGCGGAACACATCACCGGAGGCCGAACCGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGcc                                                                                                                                                                                               >  2:663945/1‑234 (MQ=255)
                                             gACGGGCCAGCTGACGTAGCTGGGCGTTGGCGTGGTGCAGCAGCAGGTCAAGGTTGTAGCCTTGGGCGATCTTGCGGAACACATCACCGGAGGCCGAACCGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTc                                                                                                                                                               >  1:1234908/1‑235 (MQ=255)
                                                     aGCTGACGTAGCTGGGCGTTGGCGTGGTGCAGCAGCAGGTCAAGGTTGTAGCCTTGGGCGATCTTGCGGAACACATCACCGGAGGCCGAACCGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGc                                                                                                                                                      >  1:1022939/1‑236 (MQ=255)
                                                      gatGACGTAGCTGGGCGTTAGCGTGGTGCAGCAGCAGGTCAAGGTTGTAGCCTTGGGCGATCTTGCGGAACACATCACCGGAGGCCGAACCGATCAGGGCATTGAGGCAGGCCCAGCGTGGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCCGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGc                                                                                                                                                      >  1:793760/3‑235 (MQ=255)
                                                                  gggCGTTGGCGTGGTGCAGCAGCCGGTCAAGGCTGTAGCCTTGGGCGATCTTGCGGAACACATCACCGGAGGCCGAACCGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTgctgct                                                                                                                                                                <  1:566223/213‑1 (MQ=255)
                                                                                          ggTCAAGGTTATGGCCTTGGGCGATCTTGCGGAACACATCACCGGAGGCCGAACCGATCAGGGCATTGTGGCGGGCCCAGCGTTTCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTGGTCGGcc                                                                                                                                                                                               <  2:798587/158‑1 (MQ=255)
                                                                                             cAAGGTTGTAGCCTTGGGCGATCTTGCGGAACACATCACCGGAGGCCGAACCGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGc                                                                                                                                                                                                >  1:103970/1‑154 (MQ=255)
                                                                                             cAAGGTTGTAGCCTTGGGCGATCTTGCGGAACACATCACCGGAGGCCGAACCGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAgg                                                                                                                 >  2:1203392/1‑233 (MQ=255)
                                                                                                 ggtGTCGCCTTGGGCGCTCTTGCGGAACACATCCCCGGAGGCCGAGCCGATCAGGGCATTGCGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGgtcgtc                                                                                                                                                                                                   <  2:15161/145‑1 (MQ=255)
                                                                                                        ccTTGGGCGATCTTGCGGAACACATCACCGGAGGCCGAACCGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCgtggt                                                                                                   >  1:6472/1‑236 (MQ=255)
                                                                                                              gCGATCTTGCGGAACACATCACCGGAGGCCGAACCGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTc                                                                                                >  2:1110946/1‑233 (MQ=255)
                                                                                                               ccttgggcgatcttgcggaacacatcaccggaggccgaaccgatcagggcattgaggcgggcccagcgttgccactgctggtgggcctgttcgatttctgccgccagcaGGCC‑‑CACCAGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGgc                                                                                                                                                                           <  1:715182‑M2/46‑1 (MQ=255)
                                                                                                               ccttgggcgatcttgcggaacacatcaccggaggccgaaccgatcagggcattgaggcgggcccagcgttgccactgctggtgggcctgttcgatttctgccgccagcaGGCC‑‑CACCAGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGgc                                                                                                                                                                           >  2:715177‑M2/110‑155 (MQ=255)
                                                                                                               cGATCTTGCGGAACACATCACCGGAGGCCGAACCGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACg                                                                                            >  1:470795/1‑236 (MQ=255)
                                                                                                                      gCGGAACACATCACCGGAGTCCGAACCGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGc                                                                                                                       <  1:33948/202‑1 (MQ=255)
                                                                                                                                  aCCGGAGGCCGAACCGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAg                                                                                                                             >  2:655417/1‑184 (MQ=255)
                                                                                                                                    cGGAGGCCGAACCGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTc                                                                                                                                             <  2:490639/166‑1 (MQ=255)
                                                                                                                                          ccGAACCGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCgg                                                                                                                                  <  1:1207584/171‑1 (MQ=255)
                                                                                                                                            gAACCGATCAGGGCATTGAGGCGGGCCCAGCGTTGTCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGgc                                                                >  1:1140144/1‑235 (MQ=255)
                                                                                                                                                cGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGc                                                                                                                                                      >  2:476028/1‑145 (MQ=255)
                                                                                                                                                cGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGt                                                           >  1:1217802/1‑236 (MQ=255)
                                                                                                                                                    cAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGGGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGGCGGCGGCCTGTTGTGCACGCAATTCGGCGCATtgc                                                                                                                                                                    >  1:233719/1‑127 (MQ=255)
                                                                                                                                                     aGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAg                                                                  >  2:822891/1‑224 (MQ=255)
                                                                                                                                                            tGAGGCGGGCCCCGCGTTGCCGCTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGcc                                                                                                                                                     <  2:229315/134‑1 (MQ=255)
                                                                                                                                                            tGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTg                                                >  1:88118/1‑235 (MQ=255)
                                                                                                                                                               ggCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAggg                                                                                                                >  2:1063222/1‑168 (MQ=255)
                                                                                                                                                               ggCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAggg                                                                                                                <  1:1063227/168‑1 (MQ=255)
                                                                                                                                                               ggCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGAt                                                                         >  1:1025622/1‑207 (MQ=255)
                                                                                                                                                                 cGGGCCCAGCGTTGCCACTGCTGGTGGGGCTGGTCGATTTCTGCCGCCCGCGCCTGGTGGGCCTGCTggcggc                                                                                                                                                                                                             >  2:138218/1‑73 (MQ=255)
                                                                                                                                                                     aggcgggcccagcgttgccactgctggtgggcctgttcgatttctgccgccagcaGGCC‑‑CACCAGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCATCAcgcg                                    <  1:349259‑M2/181‑1 (MQ=255)
                                                                                                                                                                      ggcgggcccagcgttgccactgctggtgggcctgttcgatttctgccgccagcaGGCC‑‑CACCAGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCgt                                                                                                      <  2:1136656‑M2/115‑1 (MQ=255)
                                                                                                                                                                      ggcgggcccagcgttgccactgctggtgggcctgttcgatttctgccgccagcaGGCC‑‑CACCAGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCgt                                                                                                      >  1:1136661‑M2/55‑169 (MQ=255)
                                                                                                                                                                      ccAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGGGgtgt                                                                                                 >  2:990052/1‑176 (MQ=255)
                                                                                                                                                                       cAGCGTTGCCACTGCTGGGTGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCGGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTca                                                        <  2:1082890/216‑2 (MQ=255)
                                                                                                                                                                          cGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCt                                                                                                                      <  2:585740/151‑1 (MQ=255)
                                                                                                                                                                          cGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGAcc                                                                                                        <  2:384529/165‑1 (MQ=255)
                                                                                                                                                                                 aCTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCt                            >  1:470298/1‑234 (MQ=255)
                                                                                                                                                                                   tgctgGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCgtggt                                                                                                   <  2:316417/161‑1 (MQ=255)
                                                                                                                                                                                     ctgGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAg                                                                  <  2:6472/192‑1 (MQ=255)
                                                                                                                                                                                          gggCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCtttt                 >  1:615940/1‑236 (MQ=255)
                                                                                                                                                                                           ggCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCtt                   >  2:1247845/1‑233 (MQ=255)
                                                                                                                                                                                                gTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCcgcg           >  1:858434/1‑236 (MQ=255)
                                                                                                                                                                                                    gATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGt                                                                                          <  2:1124318/153‑1 (MQ=255)
                                                                                                                                                                                                    gATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTGTTTCCGCGCCtt       >  1:1011250/1‑236 (MQ=255)
                                                                                                                                                                                                     aTTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGcc                                                                                                                                                                                               >  2:444776/1‑51 (MQ=255)
                                                                                                                                                                                                     aTTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGcc                                                                                                                                                                                               <  1:444781/51‑1 (MQ=255)
                                                                                                                                                                                                        tCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCtt                                                                                                             <  2:985653/130‑1 (MQ=255)
                                                                                                                                                                                                          tGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCgcgc                                                     <  2:1228709/184‑1 (MQ=255)
                                                                                                                                                                                                          tGCCGCCAGCGCCTGGTGGGCCGGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCaaa  <  1:893052/235‑1 (MQ=255)
                                                                                                                                                                                                            ccgccAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCGCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAg                                                                                                                             <  2:978917/110‑1 (MQ=255)
                                                                                                                                                                                                              gccAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATgg                      >  1:235869/1‑211 (MQ=255)
                                                                                                                                                                                                                 aGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCcgcg                                                                             <  1:87492/153‑1 (MQ=255)
                                                                                                                                                                                                                   cgcCTGGTGGGCCTGGTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGTCGCATGGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGAtg                                                                        <  2:152854/156‑1 (MQ=255)
                                                                                                                                                                                                                   cgcCTGGTGGGCCTGCTGGCGGCGCTGGTCGGCGGCCTGTTGTGCACGCAATGCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGt                                                           <  2:351468/169‑1 (MQ=255)
                                                                                                                                                                                                                             gCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCaaa  <  2:588705/216‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                              
GCACCAGCAGGCCCAGGGCGCTGCCGCCGCGCTTGAGGCGGTAGCGACGGGCCAGCTGACGTAGCTGGGCGTTGGCGTGGTGCAGCAGCAGGTCAAGGTTGTAGCCTTGGGCGATCTTGCGGAACACATCACCGGAGGCCGAACCGATCAGGGCATTGAGGCGGGCCCAGCGTTGCCACTGCTGGTGGGCCTGTTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAA  >  AE015451/2296803‑2297239

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