Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 5,674,751 +G intergenic (+113/+27) rhlE‑II → / ← yceI ATP‑dependent RNA helicase/polyisoprene‑binding protein, acidic stress response factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154515,674,7501.G96.0% 168.7 / ‑0.3 51intergenic (+112/+28)rhlE‑II/yceIATP‑dependent RNA helicase/polyisoprene‑binding protein, acidic stress response factor
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (30/18);  total (33/18)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CCCTGCGCGTGACGATGCTGCCCGGCAAGAGCCGGCCGTGCGCAACCCGCGTCAGCCGGAAAAACAGCCAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGC‑G‑CGCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGT‑CGCCAGTCACGTCAGCCGTCTTGTCGCCAGTAGACTTGACGCTGGTGGAGACGAATTTGGCG  >  AE015451/5674455‑5675042
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                        
ccatGCGAGTGACGAGGCGGCCCGGCAAGAGCGGGCCGTGCGGAACCCGCGTCAGCCGGAAAAACAGCCAGTGATCATTCGAAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCGGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAgc‑ggccc                                                                                                                                                                                                                                                                                                     <  2:294535/298‑5 (MQ=255)
            cctcttccgatctcAAGAGCCGGCCGTGCGCAACCCGCGTCAGCCGGAAAAACAGCCAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCa                                                                                                                                                                                                                                                                                       <  1:228739/288‑1 (MQ=255)
                 gtgcCCGGCAAGAGACGGCCGTGCGCAACCCGCGTCAGCCGGAAAAACAGCCAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTgg                                                                                                                                                                                                                                                                                  <  1:76614/300‑1 (MQ=255)
                       ggCAAGAGCCGGCCGTGCGCAACCCGCGTCAGCCGGAAAAACAGCCAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGaaa                                                                                                                                                                                                                                                                            >  1:524532/1‑301 (MQ=255)
                         cAAGAGCCGGCCGTGCGCAACCCGCGTCAGCCGGAAAAACAGCCAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGCAAtt                                                                                                                                                                                                                                                                          >  2:332940/1‑301 (MQ=255)
                         cAAGAGCCGGCCGTGCGCAACCCGCGTCAGCCGGAAAAACAGCCAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTCTTCaagatcggaaaggg                                                                                                                                                                                                                                                                          >  2:228649/1‑288 (MQ=255)
                          aaGAGCCGGCCGTGCGCAACCCGCGTCAGCCGGAAAAACAGCCAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTa                                                                                                                                                                                                                                                                         <  1:187123/301‑1 (MQ=255)
                                                       ccGGAAAAACAGCCAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGC‑AAGCGG‑CCCCC‑‑CCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTAGACACTGATAtc                                                                                                                                                                                                                                           >  1:701507/1‑301 (MQ=255)
                                                        cGGAAAAACAGCCAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtc                                                                                                                                                                                                                                           >  1:60518/1‑301 (MQ=255)
                                                              aaCAGCCAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGACTTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtc                                                                                                                                                                                                                                     >  1:381275/1‑301 (MQ=255)
                                                                       tGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCCAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCCGACATTCCTTCTGACTCCCACCTTCggccacgcattccagggtccgggtctt                                                                                                                                                                                                                            >  2:101344/1‑275 (MQ=255)
                                                                             ttCGCAAGGAGGCAAAGCTGGACCGTTGCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGt                                                                                                                                                                                                                      <  1:569435/301‑1 (MQ=255)
                                                                                          cacgCTGGACCGTTTCCCGCCGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGC‑AAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCtt                                                                                                                                                                                                        <  2:701208/298‑1 (MQ=255)
                                                                                             gcgggACCGTTTCCCGACGGCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGGACCGAGGGGGCGGCTTCGCAGCACTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTg                                                                                                                                                                                                      <  2:588965/298‑1 (MQ=255)
                                                                                                                 cTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGCGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGATCGGCCTGGTGTCGTGATAGCGCCGCAAGGCAG‑CCCCC‑‑ACGGTTTTCACCTGACAGAAATTACTTCTCTTTCACACCTTCGACACTGTTATCCAGGCTCAGGCTCTCCGAAGTAACACCCGGGCCTTTGATTCCGAGATCgttcaggtt                                                                                                                                                                                  >  2:148406/1‑301 (MQ=255)
                                                                                                                        gcTGGATGACGTGCCAAGCCGCCCACGCTGTGAGCGCCCAGCGCTGATGACCCGCAAGGGCTGATAGAAGGACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtg                                                                                                                                                                           <  2:235270/300‑1 (MQ=255)
                                                                                                                               gaccggcaaagccgacccagctgtaaGCGCACAGCGCTGCTGGCCAGGAAGGGCTGATAGAAGTACCTAGGGGGCTGCTTAGCAGACCGATCACGAAACAAGGCCGATCCTGCAGAGAATTTCTTGAGCTCGACAGGAGCGGGCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCGTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGtggtggtgg                                                                                                                                                                    <  2:346332/276‑1 (MQ=255)
                                                                                                                                     ccAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGa                                                                                                                                                              <  2:524311/301‑1 (MQ=255)
                                                                                                                                              ccACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCTTTCAGGTTCAGCGTTGTGGTTGCGTTGAAGCCAgcgc                                                                                                                                                     >  2:374039/1‑301 (MQ=255)
                                                                                                                                                   cgGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTACAGAGAATTTCTTGATCTCTACAGGAGCGGACTTGTGTTGTGATAGGGCCGCAAAGC‑GCCCCCC‑‑TCGGCTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATTCAGCTCCAGGGTCTGCCAAGTAGCACCCGGCCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAgcacacccc                                                                                                                                                >  2:152134/1‑294 (MQ=255)
                                                                                                                                                                    ggctgAACCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCtt                                                                                                                               <  2:622336/300‑1 (MQ=255)
                                                                                                                                                                                    cTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCCGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGTCCGCCTGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAg                                                                                                               >  2:10775/1‑301 (MQ=255)
                                                                                                                                                                                     tGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAgg                                                                                                              >  2:344542/1‑301 (MQ=255)
                                                                                                                                                                                              tACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAATGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGCCCTTGTGTCGTGATAGGGCCGCCAAGCGG‑CCCCC‑‑TCGGTTCTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGTCCTTTGATGCCCAAGCCGTTCAGGTTCAGCGTGGTTGTGTCGTTGAAGCCAGCGCGCTCGCCACCCCACGGGCCCTTTCCTTCACCGTTgcatgtcgcctcga                                                                                                     >  1:541754/1‑288 (MQ=255)
                                                                                                                                                                                                ccGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCCGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACGGTTGAAGGTAGCCTTTAAc                                                                                                   >  2:156072/1‑301 (MQ=255)
                                                                                                                                                                                                                  aGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTCGTGGCGTTGAAGCCATCGCGCTCGCCACCCCCTGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTTAACGTCACCGGCTTGGTGaca                                                                                 >  1:514933/1‑301 (MQ=255)
                                                                                                                                                                                                                   gCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTTGTGGTGGCGTTGAAGCCAGTGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGGCACCGGCTTGGTGacac                                                                                >  1:459018/1‑301 (MQ=255)
                                                                                                                                                                                                                   gCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGacac                                                                                >  1:327820/1‑301 (MQ=255)
                                                                                                                                                                                                                   ccccTATCACGACACAAGGCCGCTCCTGCATAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGacac                                                                                <  1:362140/300‑1 (MQ=255)
                                                                                                                                                                                                                            cGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGCACTTGTGTCGTGATACCGCCGCAAAGCGG‑CTCCC‑‑TCGGTTTTCAGCGGGCGGCAACTACTTCTGCTTCACACCTCTTACACTGATATTCAGGTCCAGGGTCTGCGAAGTAGTCCCCGGGCCTTTGATGCCGAAGTCGTTCACGTTCAGCGTGGTGGTGGCGTCGAAGCCAGCGCGCTTGCGACCCCATGGCTCCTTGCCTCCACCTTTGAAGGTAGCCTTGAACGCCACGGGCTTGGTGACACtgtgaaggg                                                                       >  1:37910/1‑296 (MQ=255)
                                                                                                                                                                                                                              acacAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGACCGTCACCGGCTTGGTGACACCGTGCATGGTc                                                                     >  1:138123/1‑301 (MQ=255)
                                                                                                                                                                                                                                  aaGGCCGGTCCTGCAGAGAATTTCTTGATCACTACAGGGGCGGCCTTTGGTCGTGATAGGGCCGCAAAGCGG‑CCGCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGt                                                                 <  2:60487/301‑1 (MQ=255)
                                                                                                                                                                                                                                       cctacacgacgcccttccgatcttctcTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGT‑CGCCa                                                           <  2:249701/278‑1 (MQ=255)
                                                                                                                                                                                                                                        gCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGTTAGCCTTGAACTTCACCGGCTTGGTGACACCGTGCCTGGTCAGGT‑CGCCAg                                                          >  2:587718/1‑301 (MQ=255)
                                                                                                                                                                                                                                         cTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTTAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGT‑CGCCAGt                                                         >  1:11959/1‑301 (MQ=255)
                                                                                                                                                                                                                                            cTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGT‑CGCCAGTCAc                                                      >  2:315160/1‑301 (MQ=255)
                                                                                                                                                                                                                                             tGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGT‑CGCCAGTCACg                                                     <  1:156142/301‑1 (MQ=255)
                                                                                                                                                                                                                                                  agaATTTCTTGATCTCTACAGGAGCGGCCTTGTGGCGTGATATGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGT‑CGCCAGTCACGTCAGc                                                <  2:176938/301‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaTTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTGCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGT‑CGCAAGTCACGTCAGCCg                                              >  2:329128/1‑301 (MQ=255)
                                                                                                                                                                                                                                                             aTCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTAGTGACACCGTGCATGGTCAGGT‑CGCCAGTCACGTCACCCGTCTTCTCGc                                     >  2:744275/1‑301 (MQ=255)
                                                                                                                                                                                                                                                              tctcTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGT‑CGCCaagatcggaagagcacacgtctga                                    >  1:249801/1‑278 (MQ=255)
                                                                                                                                                                                                                                                                tctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGT‑CGCCAGTCACGTCAGCCGTCTTGTCGCCAg                                  <  1:333067/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGT‑CGCCAGTCACGTCAGCCGTCTTGTCGCCAGTaa                               >  2:189357/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                       gAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGACGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGT‑CGCCAGTCACGTCAGCCGTCTTGTCGCCAGTAGACtt                           >  1:238471/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                           ggCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTTCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGACAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGT‑CGCCAGTCACGTCAGCCGTCTTGTCGCCAGTAGACTTGACg                       >  1:557067/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                           ggCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGTTGACACCGTCCATGGTCACGT‑CGCCAGTCACGTCAGCCGTCTTTTCGCCAGTAGACTTGAcc                       >  2:682630/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                             ccTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGT‑CGCCAGTCACGTCAGCCGTCTTGTCGCCAGTAGACTTGACGCt                     <  1:315279/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                               ttGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGT‑CGCCAGTCACGTCAGCCGTTTTGTCGCCAGTAGACTTTCCGCtgg                   >  2:75626/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                       tGATAGGGCCGCAAAGC‑G‑‑GCCCTTTCGGTTTTCAGCTTGCAGAAATTACTTCTGCTTCACACCCTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGGTCAGGTTCAGCGTGGTGGTGGCGTTGACGCCATCGCGCTCGCCACCCCATGCGTCCTTGCCTTCACCGCTGAAGGCAGCTTTGCACGTCACCGTCTTGCTCACACACTGCACGGTCAGGTACACTGGTCA‑GTCAGCCGTCTTGTCGTCAGTGGACtggtgcatggtggagacg           >  2:714846/1‑284 (MQ=255)
                                                                                                                                                                                                                                                                                           aGGGCCGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGT‑CGCCAGTCACGTCAGCCGTCTTGTCGCCAGTAGACTTGACGCTGGTGGAGACGAAtt       <  1:189437/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            gagCCGCAAAGCGG‑CCACC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTGCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTGGAAGCCAGCGCGATCTCCACCCCATGGGTCCTTGCCTTCACAGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGT‑TGCCAGTTCCGTCATCCGTCTTTTCGCCAATAGAATTACCGCTGGTGAAGACGAAttt      >  2:737730/3‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                cGCAAAGCGG‑CCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGT‑CGCCAGTCACGTCAGCCGTCTTGTCGCCAGTAGACTTGACGCTGGTGGAGACGAATTTGGCg  >  1:732397/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                        
CCCTGCGCGTGACGATGCTGCCCGGCAAGAGCCGGCCGTGCGCAACCCGCGTCAGCCGGAAAAACAGCCAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGC‑G‑CGCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGTCAGGT‑CGCCAGTCACGTCAGCCGTCTTGTCGCCAGTAGACTTGACGCTGGTGGAGACGAATTTGGCG  >  AE015451/5674455‑5675042

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