Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451196,3720AG24.0% 56.5 / 24.6 51K626K (AAA→AAG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (16/22);  new base G (5/7);  total (21/30)
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.44e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAC‑CCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA  >  AE015451/196076‑196670
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                             
cAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACTGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTCACGGCCAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCTTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAccc                                                                                                                                                                                                                                                                                                           <  1:420450/301‑1 (MQ=255)
            gATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACGACGATAAACGACTCGGTCGACACCACTACCGCGACCCGGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACACAGGGCAGCgccc                                                                                                                                                                                                                                                                                               >  1:497727/1‑298 (MQ=255)
               ttCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGAtt                                                                                                                                                                                                                                                                                            <  1:459628/301‑1 (MQ=255)
                 ccAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTc                                                                                                                                                                                                                                                                                          >  1:450324/1‑301 (MQ=255)
                   aGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGAGTCGGTAGACACCACTACCGCGACTCTGACCGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑Ca                                                                                                                                                                                                                                                                                       >  1:146909/1‑301 (MQ=255)
                   aGACACCGGCCACTGAAGTCTACGACCACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGTAACTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTC‑Ca                                                                                                                                                                                                                                                                                       <  2:527069/301‑1 (MQ=11)
                      ccccGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAACTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGAc                                                                                                                                                                                                                                                                                    <  1:726071/301‑1 (MQ=255)
                            ccAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCCGTAACCGAAGGCTGCGTGATCCCTTACACCGCGACCCTGAGCACTCCTGCCTAGACGACGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCTAGGGCAGCGTCGATTTC‑CAGACCCCGGc                                                                                                                                                                                                                                                                              >  1:216328/1‑301 (MQ=14)
                               aTGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGaa                                                                                                                                                                                                                                                                           >  2:90904/1‑301 (MQ=255)
                                tGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAAt                                                                                                                                                                                                                                                                          <  1:90942/301‑1 (MQ=255)
                                         ctacaacGGTTCGACTGTCAGAGTCACCAACGAGAACGCCACAGGCGGTAATTGCGAGCAACTGACCCCGATCCCGACCCCGGCTCACACCACGATCAACGACTCGGTCGATACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTTACTCTGTTTAATGCTACCCAGACGACCGAGATC‑GACACTCTGTCCCACCGCGACGCTCTTCCGGTCGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCCAATGACGTCTac                                                                                                                                                                                                                                                                 <  2:362886/299‑1 (MQ=2)
                                               ttccgatctggTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACaacaac                                                                                                                                                                                                                                                           <  1:510101/292‑1 (MQ=255)
                                                 gtgcGACCGTCAGCGTCACCATCGAAAACGCCACAGGTGGGAATTGCGAGCAACGGACCCCGAACCGGACCACGGATCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCGTGCCCAGACGCGGGTG‑CACGTGACTCCGTCCAACGGCCAAGCCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTC‑CAGACCCCGGCGAATGACGTCTACAACAACgg                                                                                                                                                                                                                                                         <  1:368938/298‑1 (MQ=11)
                                                    cGACGGTCAGCCTCACCATCGAGAACGCCACGGGCGGTAATTCCAAGCAACTGACCCCTATCCAGACGCCAGCTCAGACACCGATCAACGACCCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTc                                                                                                                                                                                                                                                      <  1:260355/301‑1 (MQ=14)
                                                     gACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑g                                                                                                                                                                                                                                                    >  1:667813/1‑301 (MQ=255)
                                                        ggTCAGCGTCACCATCGAGGACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACTGCGAGCCCGTCGGTCACCGAAGGCGTCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACt                                                                                                                                                                                                                                                 <  1:179078/301‑1 (MQ=255)
                                                        ggTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAACCCCAGGGCCGCGTCGATTTC‑CAGACCCCGGCCAATGACGTCTacacaaacagataggaa                                                                                                                                                                                                                                                  >  2:509893/1‑287 (MQ=255)
                                                         gTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GAccc                                                                                                                                                                                                                                                >  1:353110/1‑299 (MQ=38)
                                                            aGCGGCACCATCGAGACCGCCACGGGCTGCAATTTCGAAAAACTGACCCCCAACCCGACGCCGCCACAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCCGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACCGTCa                                                                                                                                                                                                                                             <  1:422873/301‑1 (MQ=11)
                                                             gcgccACCATAGAGAGCGCGACGGGCGGCAAGTTCGGACAACGGGACCCGAACCCGACGCCGGGTCAGACCAGGAGCAAAGAATCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATAACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACCGTCAg                                                                                                                                                                                                                                            <  2:84259/297‑1 (MQ=11)
                                                              ggttaccctggggaACGCCACGGGCGGGAAAGTAGAGCAAATGACCCCGATCCCGTCGCCGGGTCAGAGACCGAACAACGGCTCGGGCGCCACCACTGGCGTTAACCCGACGGCTATCCCGTCGTTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAgt                                                                                                                                                                                                                                           <  1:493892/289‑1 (MQ=14)
                                                               gTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATCACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGACTAACGTCTacagtcacggctcgatcgtcagcg                                                                                                                                                                                                                                           >  1:267070/1‑280 (MQ=11)
                                                                                    ggcagTAAGTTCGTGCAACTGACCCAGAACCCGACGGCGGCTGAGAACAAGATCAACGACTCGGTCGACACCACTGCCGCGACCCGGACGGCGAGCCCGTCGGTCACGGAAGGCGACGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGAACCTGTCCAGCGGCCAAACCATTACCGTTGAAGCGGGCAAAAGCAATGACAGCGTCGATTTC‑CCGACCCCGGCCAATGACGTCTACAACAACGGTTC‑GACTGTCAGTGTCACCATCGAGAACGCCACcg                                                                                                                                                                                                                     <  1:104150/297‑1 (MQ=14)
                                                                                    ggCGGTAATTTCGAGCAACTGACCCCGAACCCGCCGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACCGTCAGCGTCACCATCGAGAACGCCAcag                                                                                                                                                                                                                     <  2:352975/301‑3 (MQ=12)
                                                                                          aaTTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAggagaccggacgaggatcctacaggcaagg                                                                                                                                                                                                               >  2:418806/1‑272 (MQ=11)
                                                                                                      cTGACCCCGAACCCGACGGCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCCGTGCCCTAGAAGCCGGCACAACCCAGCGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAAc                                                                                                                                                                                                   <  2:450134/301‑1 (MQ=14)
                                                                                                      cTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCCAACCATTACCGTTTAAGACGTCAAAACCCCGGGCATAGTCGATTAC‑GCGACCCCGTCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGTATCACCATCGAGAACGCCACCGGCGTTAATTTCGAGCAAc                                                                                                                                                                                                   >  1:568276/1‑301 (MQ=14)
                                                                                                           cccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTACAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGATCATCTGAc‑c                                                                                                                                                                                             >  1:50662/1‑301 (MQ=255)
                                                                                                                 cccGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACTGCGAGGCCGTCGGTCACCGACGGCGTCGAGATCACTTACACCGTGACCCTCAGCAACCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAC‑CCCGAAt                                                                                                                                                                                       <  1:632180/301‑1 (MQ=255)
                                                                                                                               aGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAC‑CCCGAATCCGACGCCGGcccc                                                                                                                                                                         >  1:331123/1‑298 (MQ=255)
                                                                                                                                          aacggtgcgggcGACACCACTACCGCGGCCCGGACGGCGAGCCCGTCGGTGACCGAAGGCGGCGGGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTTACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCTGGCTAATCACACGGCAGCCACCGTGAC‑CAAACCCCTGCCCATAACGGCCAC‑TCTACCGATCTGACGGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAC‑CCCGAATCCGACGCCGGCTCAGACCACGATCa                                                                                                                                                              <  2:667530/290‑1 (MQ=2)
                                                                                                                                                     cGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCACGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAC‑CCCGACTCCGACGCCGGCTCAGACAACGATCAACGACTCgttc                                                                                                                                                   >  1:455417/1‑298 (MQ=255)
                                                                                                                                                           cacTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAC‑CCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGcc                                                                                                                                             <  1:681047/301‑1 (MQ=40)
                                                                                                                                                            acTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACAGCCCGACAATTCGCCTTGTAGCGGGC‑AAGTCTAGGTCCCGGTGGGCTTCGCAGACTCCGACGAATGACGACTACAACAACGGTTC‑GACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAACTTCAAGCAACTGAC‑CCCGAATCCGACCACGGCTCCAACCACGATCAACGACTCGGTCGATGccc                                                                                                                                            >  2:725523/1‑300 (MQ=2)
                                                                                                                                                                                             aCCGAAGGCGTCGTGATCAATTACACCGTGACCCTAAGCAATCCTGCCCAGACGCCGGTGACA‑GTTACACTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTTAGTGTCACCATCGAGACCGCCACCGCCGATAACTACAATCTACCGACACCCG‑ACCCGCTGCCGATCCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCa                                                                                                           <  2:330996/301‑1 (MQ=11)
                                                                                                                                                                                                               aCTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTTCAGCGGCCAAACCATGTCCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCCAATGACGTCTACAACAACGGTTC‑GACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAC‑CCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCa                                                                                         <  1:156692/301‑1 (MQ=255)
                                                                                                                                                                                                                             cctaagcaagcccttccgatcGCCGGTAACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAC‑CCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGAcc                                                                           <  1:58732/281‑1 (MQ=255)
                                                                                                                                                                                                                                   cAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAGCTGAC‑CCCGGATCCGACGCCCGCACATACGACGATCAACGACTCGGCCTATGCTACCATCCCGACCCTGACGGCGAACACGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTCACCCTgcgc                                                                     >  1:405609/1‑298 (MQ=14)
                                                                                                                                                                                                                                      gccTGCCCCGACGCCGGAGTCA‑AGGACCCTGACCAAAGGTCAAACCAGTACCGTTGAAGACGGCAAAACCCAGGCCAGCTTCGATGGC‑CAGAGCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGTGTCACCATCGAACACGCCACCGGTTTTAAATGCGACGAACTCAC‑CCCAGATCCGCAGCCGGCTCCGACCCGGATCTCCTCCGCGCTCGATGCCACAACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAAt                                                                  <  2:455227/300‑1 (MQ=9)
                                                                                                                                                                                                                                             ccggacgcacggaacaatgccccTGTCAAAAGGCCAAACCATTAGCGTTCAGACCGGCAAAACCCCGGGCTGCGTCGATTTC‑CAGACCCCGGCAAATGACGTCTACAACAACGGTTC‑GACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAC‑CCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGcc                                                            <  1:426967/281‑1 (MQ=14)
                                                                                                                                                                                                                                                gacgtgtgctcttccgatcttGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAC‑CCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAg                                                         <  1:338955/281‑1 (MQ=255)
                                                                                                                                                                                                                                                     cccgacgcgtttccgtgctcAACGGCCAAACGATTACCGTTGAAGCCGGTAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCCAATGACGTCTACAACAACGGTTC‑GACGGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAC‑CCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                    <  2:299174/282‑1 (MQ=11)
                                                                                                                                                                                                                                                                    tGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAC‑CCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCGGCCCAGagatcggaagagcgtcgtgt                                     >  2:338828/1‑282 (MQ=255)
                                                                                                                                                                                                                                                                           cGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAC‑CCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAAc                              <  1:725832/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                               cAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAATTTCGAGCAACTGAC‑CCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGCCGAGCCCGTCGGTCACCGAAGGAGGCGGGAACACCTACACCGTGACCCTAGGCAATCCTGCCCCGAcacccgtagacacggccgggtctgacgtag                          >  2:498626/1‑272 (MQ=11)
                                                                                                                                                                                                                                                                                    gaTTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAC‑CCCGCACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTCCACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGCGACCCTGTCCACCGGCCAAAcc                     >  1:328883/2‑301 (MQ=11)
                                                                                                                                                                                                                                                                                    ccatACCGTTGAAGCCTGCAATACCCAGGTCACCGTCGATTTC‑CACACCCCGGCAAATGACGTCGACATCCTCGGTTC‑AACCGTCAGCGTCACCCTCCAGAACGCCACAGGCGGTAATTTCGAGCGACTGAC‑CCCGTATCCGACTCCTGCTCAGACCACGATCCACGACTCGGTCCATGCTACCGCCGCTACCCTTCCGGCGAGCCCGTCTGTCACCGAATGCGGCGTAATCACCTACACCTTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAAcc                     <  1:80403/298‑1 (MQ=2)
                                                                                                                                                                                                                                                                                       cgctcttcggatctgCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAC‑CCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCAtt                  <  2:528221/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         ccgttaaaggaggCAAAACCCAGGGCAGACTCGATTTC‑CAGGCCCGTGGGAATGACCTCTACCACAACGGTTC‑GACTGTCAGAGTCAACATCGAGAACGCCACAGGCGGTAATTTCGAACAACTGTC‑CCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCGCCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTAAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTAc                <  2:360696/290‑1 (MQ=11)
                                                                                                                                                                                                                                                                                               aaGCCGGCAAAACCCAGGGCAGCGTCGACTTC‑CAGACCCCGGCGTATGACGTCTACAACAACGGTTC‑GACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAC‑CCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGa          >  1:363022/1‑301 (MQ=11)
                                                                                                                                                                                                                                                                                                     gCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAC‑CCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTagatcggaagagcc    >  1:528442/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                                       aaaaCCCAGGGCAGCGTCGACTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAC‑CCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCCACACCGTCACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCCTTGAAGCCGGCa  >  1:722378/1‑301 (MQ=11)
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                             
CAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACA‑GTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTC‑CAGACCCCGGCGAATGACGTCTACAACAACGGTTC‑GACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAC‑CCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA  >  AE015451/196076‑196670

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