Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 677,006 C→G 8.0% A99G (GCG→GGG)  ybeR → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913677,0060CG 8.0% 201.3 / 2.4 86A99G (GCG→GGG) ybeRuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (42/35);  new base G (3/5);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

TTACCTTCGCATTCATCACCCAAAATGTATTTAGGTACATGCTTTTTTGACCTTTCTTCCTCATGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACGGGCTCTACGTTATTATAGTAACTGGCGGCAATACTTTGCCGGTTATACCTTTGGTCGGCAATATTGGCA  >  NC_000913/676780‑677228
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                               
ttACCTTCGCATTCATCACCCAAAATGTATTTAGGTACATGCTTTTTTGACCTTTCTTCCTCATGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGTTGAAC‑TAAATGAGCAAGGTCAGGCTTATgcg                                                                                                                                                                                                                                >  2:172955/1‑228 (MQ=255)
ttACCTTCGCATTCATCACCCAAAATGTATTTAGGTACATGCTTTTTTGACCTTTCTTCCTCATGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATgc                                                                                                                                                                                                                                 >  2:105681/1‑227 (MQ=255)
      tcccattgctcccccAAAATGTATTTAAGTACGTGTTTTTTTGCCTTTGCTTCCTCATGGGGAATTGATGACCGGTGTGACCTGCTGCGGACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAAttt                                                                                                                                                                                                                         <  1:169166/216‑1 (MQ=255)
         cATTCATCACCCAAAATGTATTTAGGTACATGCTTTTTTGACCTTTCTTCCTCATGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGGG‑CAATTTGTc                                                                                                                                                                                                                      >  1:396833/1‑228 (MQ=255)
              atcaCCCAAAATGTATTTAGGTACATGCTTTTTTGACCTTTCTTCCTCATGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGGGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTAGGCG‑AAATTTGTTCCCAg                                                                                                                                                                                                                 >  1:144913/1‑228 (MQ=255)
                            aTTTAGGTACATGCTGTTTGGACCGTTCTTCCTCATGGGGAATTGATGACCGGGATGACCTGCGGCGGACAATTCACCGAATGATCGATACAGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCt                                                                                                                                                                                                   <  2:538781/228‑1 (MQ=255)
                             tttAGGTACATGCTTTTTTGACCTTTCTTCCTCCTGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCCC‑ACCTGGTTTCGTTATTCGCCATGTGTATGGCGTGCTTATCTTCCTGAAT‑TACTTGAGCACGGTCAGGCGTA‑GCGTCACTTTGTCGCCATCACAGCCAAATCCtt                                                                                                                                                                                                  >  1:191784/1‑227 (MQ=255)
                             tttAGGTACATGCTTTTTTGACCTTTCTTCCTCATGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCTTGTGAATGGCGTGACTATCTTACTGAAT‑TAAATGATCAAGGCCGGGCTTTTGCG‑CAGTGTGCCGCCAGCAccgccacatgtgg                                                                                                                                                                                                  >  1:587309/1‑220 (MQ=255)
                             tttAGGTACATGCTTTTTTGACCTTTCTTCCTCATGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGATC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACAGCCGAATGCtg                                                                                                                                                                                                  >  2:585974/1‑228 (MQ=255)
                             tttAGGTACATGCTTTTTTGACCTTTCTTCCTCATGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCtg                                                                                                                                                                                                  >  2:397981/1‑228 (MQ=255)
                             tttAGGTACATCCTTTTTTGACCTTTCTTCCCCCTGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGCTAACGGTCATGCCGCTCG‑GTTGGCCGGGTTTTATCAC‑CGCTGGTTTCGTTATTCCCCTTGTGATTGGGGGGACTATCTTGCTGAAC‑TAAATGTGCCTGGACAGGGTTATGCG‑CACTTTGTTGCCAGCAccccccagcgctg                                                                                                                                                                                                  >  1:273263/1‑220 (MQ=255)
                                  gTACCTGCTTTTTTGACCTTTCTTCCTCATGGGGAATTGATGACCGTGTTGGCCTGCTGCGCCCAATTCCCCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCg                                                                                                                                                                                             <  2:587309/228‑1 (MQ=255)
                                   tACATGCTTTTTTGCCCTTTCTTCCTCAGGGGGAATTGATGACCGTGATGACCTGCTGCGCCCAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCg                                                                                                                                                                                             <  2:79498/227‑1 (MQ=255)
                                     cATGCTTTTTTGACCTTTCTTCCTCATGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGCCTATCTTGCTGAAC‑TAAATGAGCACGGTCAGGCTTATGCG‑AAATTTTTCGCCAGCACCGCCGAATGCGGTGGCGAAg                                                                                                                                                                                          >  2:507966/1‑228 (MQ=255)
                                     cATGCTTTTTTGACCTTTCTTCCTCATGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAGTGAGCAAAGTCAGGCTTATGCG‑CAATTTGTGGCCAGCACCGCCGATTGCTGTGGCGaa                                                                                                                                                                                           >  2:8825/1‑227 (MQ=255)
                                     cATGCTTTTTTGACCTTTCTTCCTCATGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CACTTTGTCGCCAGCACCGCCGAATGCTGCGGCGaa                                                                                                                                                                                           >  2:454736/1‑227 (MQ=255)
                                     cATGCTTTTTTGACCTTTCTTCCTCATGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAg                                                                                                                                                                                          >  2:341704/1‑228 (MQ=255)
                                     cATGCTTTTTTGACCTTTCTTCCTCATGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTCCTTAAC‑TAAGTGAGCAAGGTCAGGCTTATGCG‑CAATTTTTCGCCAGCACAGCCGaaacctgggtggaag                                                                                                                                                                                          >  1:201058/1‑215 (MQ=255)
                                     cATGCTTTTTTGACCTTTCTTCCTCATGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGGGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATATGTTCCCAGCACCGCCGAAAGCTGTAGGGAAg                                                                                                                                                                                          >  1:79498/1‑228 (MQ=255)
                                     cATGCGTTTCTGACCTTTCTTCCTCATGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGGTTTATCAC‑CGCCGGGTTCGTTATTCGCCATGTGGGTGGAGTGCCTATCCTGCTGAAC‑TAAATGATCAAGGTCAGGGTTATGGC‑CAATTTGTCGCGAGAACCGCCGAAGGCGGTGGCgcgg                                                                                                                                                                                          >  1:701300/1‑225 (MQ=255)
                                      aTGCTTTTTTGACCTTTCTTCCTCATGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGACC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAAGTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAg                                                                                                                                                                                          >  2:649969/1‑227 (MQ=255)
                                                           cTCATGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGAttg                                                                                                                                                                    >  2:621400/1‑227 (MQ=255)
                                                                ggggAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAGTGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAGTGCTGTGGCGAAGGGGGAATCAAAGCCGGGGCGTATGTcc                                                                                                                                                               >  1:361986/1‑228 (MQ=255)
                                                                              cGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGGGAAGGGGGAATCAAAGCCTGGGGTTATGCCCGGTTGGGGTTTTTa                                                                                                                                                 >  1:117605/1‑228 (MQ=255)
                                                                              cGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAAAGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTa                                                                                                                                                 >  2:359951/1‑228 (MQ=255)
                                                                                tgattaCCTGGTGCGCACAATTCCCCTAATTACCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTCCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAg                                                                                                                                               <  2:320100/228‑1 (MQ=255)
                                                                                tgatgaCCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTTTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTaa                                                                                                                                                >  2:104053/1‑227 (MQ=255)
                                                                                gtgtgaCCTTCTGGTCACATTTCCCCGAATGGTCGATAACGGTCATGCCGCTGG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTAGTCGCCATGGGAATGGCGTGACTATCTTGCTGAAC‑TAAAGGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCTCCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAg                                                                                                                                               <  1:189241/225‑1 (MQ=255)
                                                                                          ctgcGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑AAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAg                                                                                                                                                                              >  2:253958/1‑187 (MQ=255)
                                                                                           gtcccccatttCGCCGAATGAACGTTTATGGTCTTGCCCCTCG‑GTTGGCGGGTTTTTATCAC‑GGCTGGTTTCGTTAGTCGCCATGTGAATTGCGTGACTTTCTTGCTGAAC‑TAAATGATCAAGGTCAGGCTTATGGG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCg                                                                                                                                    <  1:172955/219‑1 (MQ=255)
                                                                                           ggcgcgcATTTCCCCGAATGATCTATAACTGTCGTGCCGCTGG‑GGTGGCGGGTTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCGGGCTTATGCG‑CAATTTGTCGCCAGTTCCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCg                                                                                                                                    <  1:585974/227‑1 (MQ=255)
                                                                                            ggagagggttataagatacagcGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCg                                                                                                                                    <  2:197590/206‑1 (MQ=255)
                                                                                            gcgccaaaTTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCg                                                                                                                                    <  2:117605/221‑1 (MQ=255)
                                                                                            gcgcACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCg                                                                                                                                    <  2:201058/227‑1 (MQ=255)
                                                                                                       aCCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCgg                                                                                                                                                             >  2:494002/1‑191 (MQ=255)
                                                                                                         cGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCGGAATGGGCGTCCTtataaattgt                                                                                                                      >  2:90939/1‑220 (MQ=255)
                                                                                                         cGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGGGACCTTCATAACTgg                                                                                                                      >  2:485149/1‑228 (MQ=255)
                                                                                                                tCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGGGTACTTGATAACTGGTTggggg                                                                                                               >  1:624171/1‑224 (MQ=255)
                                                                                                                gcGATACCTGGCATGTCGGTAG‑GTGGGCAGGGTTTTATCAC‑GGGTGGTTCCGTTAGTCGCCAGGTGAATGGCGTGACTATCTTGCTGAAC‑TAAAGGGGCAAGGGCAGGCTTAGGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTGTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCg                                                                                                               <  2:191784/227‑1 (MQ=255)
                                                                                                                 cgaacaCGGTAATGCCGTTCG‑GTTGGCAGGGTTTAAGCAC‑CGGTGGTTTCGTAATTCGCCATTTGAATGGCGTGACTATCGTGCTGAAC‑TAAATGAGCAAGGTCATGCTTATGCG‑CAATTGGTCGCCAGCACCGCCTAATGCTGTGTCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCga                                                                                                              <  2:273263/223‑1 (MQ=255)
                                                                                                                     aaCGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGCAGGGGGGATCAAAGCCTGGGTTTATGCCCGGATGGGGTTTTTAAGCGGAATGGGCGCCCTTAATAAATGGTTGAGGgagaaa                                                                                                          >  1:634673/1‑224 (MQ=255)
                                                                                                                     aaCGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATTTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCgaagaa                                                                                                          >  2:680539/1‑228 (MQ=255)
                                                                                                                       tgtgcgtcccctgtg‑gTTGGGCGGGTTTTAACAC‑CGTTGGTTTCGTTATGCGGCAGGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCgaagaaga                                                                                                        <  1:359951/214‑1 (MQ=255)
                                                                                                                       tggTCGTGCCGTCGG‑GTTGGCAGGGTTTTACCAC‑CGCTGGTTTCGTTATCCGCCATGTGAATGGCGTGACTATCTGGCTGAAA‑TAAACGAGCAAGGTCAGGCTTATGCG‑AAATTTGTCGCCAGCCCCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCgaagaaga                                                                                                        <  1:622822/227‑1 (MQ=255)
                                                                                                                       cGGTCATGCCGCTGG‑GTTGGCAGGGTTTTATCGC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCgaagaaga                                                                                                        <  2:445935/228‑1 (MQ=255)
                                                                                                                       cGGTCATGCCGCTCG‑TTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCgaagaaga                                                                                                        <  2:409550/228‑1 (MQ=255)
                                                                                                                        gttcctGCCGCTCG‑GTTTGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCgaagaaga                                                                                                        <  2:105690/222‑1 (MQ=255)
                                                                                                                            aTGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGGG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGCATCAAAGCCTGGGATTATGTCCGGTTGGGTTTTTTAAGGCGAATGGGCGTACTTAATAAccttgcttttatacaattttga                                                                                                    >  1:197590/1‑206 (MQ=255)
                                                                                                                             tGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGGG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCt                                                                                                  <  2:396833/228‑1 (MQ=255)
                                                                                                                              ggcggttggttgggCAGGTTTTTATCAC‑GGCTGGTTTCTTTATTGGCCATGGGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGGCGCGAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCt                                                                                                  <  1:454736/216‑1 (MQ=255)
                                                                                                                              gCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGAACTTAATAACTGGTTGAGCGGAGAGGAAAGCCTa                                                                                                 >  2:208138/1‑228 (MQ=255)
                                                                                                                                ttggccttgttggggggTTTTTTCCCTCGCGGGTTGCGTGAGTCGCCAGTGGAATGGCGGTATTTTCCTCCTGAAC‑TAAATGAGCAGTGTCGGGCTTAGGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTa                                                                                                 <  1:244652/214‑1 (MQ=255)
                                                                                                                                cgcctg‑gTTGGCAGGGTTTTTCTAC‑CGCTTGTTTCGTTATTCGCCATTTGAATGGCGTTACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGCCGCCAGCACCGCCGAATGCTGTTGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTAt                                                                                                <  2:634673/222‑1 (MQ=255)
                                                                                                                                cGCTCGTGTGGGCGGGTTTTTTTCCC‑CCCTGGTTAGGTTATTCCGGATTGGAACGGCGTGACTATCTTGCTAAAG‑AAACTGGGAAAGGTGGGGTTTATGGG‑CAATTTGTCGCCAGAACCGCCGATGGCTGTGGCGAAGGGGGAATCAAAGTCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTAt                                                                                                <  1:649969/228‑1 (MQ=255)
                                                                                                                                             aGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCCGGGCTTATGTCCGGATGGGTTTTTTAAGCCGACTGGTCGTACTTAATACTTGGTTGAGCGAAGAGGAAAGTCTATGGATCCAgctcgg                                                                                   >  2:2524/1‑223 (MQ=255)
                                                                                                                                                gTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGGTTTTTAAGCCGAAAGGGc                                                                                                                                     >  1:625448/1‑175 (MQ=255)
                                                                                                                                                 ttttATCCC‑CGCTGTTTTGGTTATTCGCCATGTGAAGGGGGTTACTATCTTGCGGAAC‑AAAATGAGGAAGGTCAGGCTTATGCG‑CAATTGGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCAtt                                                                               <  1:105681/228‑1 (MQ=255)
                                                                                                                                                      tCAC‑CGCTGGTTTCGTTATCCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGGGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCt                                                                          <  2:46362/228‑1 (MQ=255)
                                                                                                                                                         c‑cGGTGGTTTTGTGGTTCGCCGTGTGAAGGGGGGGAGTTTCGGGGTGAAC‑GAAATGAGCAAGGGCAGGCTGAGGCG‑AAAGTTGGCGGCGGCTGCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTGAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGCATGCCTATGGATCCAGTCGCGCATTTATATAc                                                                        <  1:680539/227‑1 (MQ=255)
                                                                                                                                                           cGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACgg                                                                      >  2:214242/1‑228 (MQ=255)
                                                                                                                                                                 ttttttttttCGCGATGGGAATGGGGTGACTTGCTTGCTGAAT‑CAAATGAGCAAGGTCAGGTTTATGCG‑AAATTTGTCGCCAGCACCGCCGAATGCTGTGTCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACGGGCTCTa                                                                <  1:407004/220‑1 (MQ=255)
                                                                                                                                                                 ttttgttgtgcgcgcggtacaTGGCGTGAATTTGTTG‑TGAACAAAAATGGGCGAGGTTAGGGTTACGGG‑AAATTGGTCGGCAGAACCGCCGGAGGCTGTTGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGAGGGGTTGTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCAAAGAAGAAAGTCTATGGATCCAGTCGCGCATACATCTGCGGGCTCTa                                                                <  1:268226/208‑1 (MQ=255)
                                                                                                                                                                 tgtagtgtttaggcgtgggAAGGGCGTGATTATCGTGCTGAAC‑TAAAGGAGCAAGGTCAGGATTATGAG‑CAATGTGTCGCGAGCAACGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACGGGCTCTa                                                                <  1:545510/210‑1 (MQ=255)
                                                                                                                                                                      ttatgccgcttgggagttgCGTGACTTGTTTGCTGAAA‑TAAGTGAGCAAGGTCGGGCGTATGCT‑GACTTGGTCGCCCACTCTGCCTAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTGCGGATGGGTTTTTTAAGCCGAATGGGTGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCAGCTACGGGCTGTACGtta                                                           <  1:191688/210‑1 (MQ=255)
                                                                                                                                                                          tCGCCATGTGAATCGCGTGACTATCTTGCTGAAC‑TAAATGAGCACGGTCAGGCTTATGCG‑CACTTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGCTGGGTTTTTTAAGCAGAATGGGCGCACTTAATAACTGGTTGAGCGAAGAAGAATGTCTATGGATCCAGTCGCGCATTCATCTCCGGGCTCTAAGTtattat                                                       >  1:168970/1‑228 (MQ=255)
                                                                                                                                                                                 ggggttGGGGTGACTTTTTTGCTTCAC‑GAAGGGAGGAAGGTCGGGCTTATGCA‑CAATTTGGCGACAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACGGGCTCTACGTTATTATAGTAACt                                                <  1:397981/223‑1 (MQ=255)
                                                                                                                                                                                  tgAATGGCGGTACTTTCTTGCTGAAC‑TAAAGTAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACGGGCTCTACGTTATTATAGTAACt                                                <  2:168970/227‑1 (MQ=255)
                                                                                                                                                                                         cGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACGGGCTCTACGTTATTATAGTAACTGGCGGCaa                                        <  2:624171/228‑1 (MQ=255)
                                                                                                                                                                                           tGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGTGCGCAGAAGAAAGGCTATGGATCCAGTCGCGCATTCATCTACGGGCTCTGCGTTATTATAGTAACTGGCGGCACTa                                      >  2:329247/1‑228 (MQ=255)
                                                                                                                                                                                              cTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATC‑AAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGCATGGGCGGACTTAATAACTGGTTGAGCGAAGAAGAAAGGCTATGGTTCCATTCGCGCATTCATTTACGGGCTCTCCGTTATTTTAGTAACTGGCGGCAATACtt                                   >  1:43094/1‑227 (MQ=255)
                                                                                                                                                                                                ttcgtggtGAAC‑TAAAGATGCAAGGTCGGGTTTAGGCG‑CAATTTTGCGCGAGCACCGGCGAATGGTGTTGCGGAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACGGGCTCTACGTTATTATAGTAACTGGCGGCAATACTTTg                                 <  1:341704/221‑1 (MQ=255)
                                                                                                                                                                                                 tCTTGCTTAAC‑TAAATGGAGAAGGTCAGGCTTATGCG‑CAATTTGGCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACGGGCTCTACGTTATTATAGTAACTGGCGGCAATACTTTg                                 <  2:144913/227‑1 (MQ=255)
                                                                                                                                                                                                 tCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGGTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGTTCCAGTTGCGCATTCATCTACGGGCTCTACGTTATTATTGTAAATGGCGGCAATACTTTGc                                >  2:35243/1‑228 (MQ=255)
                                                                                                                                                                                                      cTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCGGAATGGGCGTACTTAATAACTGGTTGAGCGGAGAAGAAAGTCTATGGAGCCAGTCGCGCATTCATCTACGGGCTCAAAGTACTTAAAGTAACCggcggccaacattgtccgggg                           >  1:16016/1‑213 (MQ=255)
                                                                                                                                                                                                         aaC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGTAAGTCTATGGATCCAGTCGCGCATTCATCTACGGGCTCTa                                                                >  1:557504/1‑188 (MQ=255)
                                                                                                                                                                                                         aaC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACGGGCTCTa                                                                <  2:557504/188‑1 (MQ=255)
                                                                                                                                                                                                         aaC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACGGGCTCTACGGTATTATAGTAACTGGCGGCAATACTTTGCCGGTtata                        >  2:571469/1‑228 (MQ=255)
                                                                                                                                                                                                         aaC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGGCGCGCATTCATCTACGGGCTCTACGTTATTAAAGTAACTGGGGGCAATACttttcacggttta                        >  1:662546/1‑218 (MQ=255)
                                                                                                                                                                                                                    ggaaGGTTAGGCTTTTGCA‑AAATTTGTCGCAAGACCCGACGAATTCTGTGGCGAAGGGGTAATCAAAGGCTGGGATTATGGCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACGGGCTCTACGTTATTATAGTAACTGGCGGCAATACTTTGCCGGTTATACCTTTGGTCg              <  1:208138/226‑1 (MQ=255)
                                                                                                                                                                                                                         gTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGCATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCgaagaag                                                                                                         >  2:298855/1‑132 (MQ=255)
                                                                                                                                                                                                                         gTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACGGGCTCTACGTTATTTTAGTAACTGGCGGCAATACTTTTCCGGTTATACCTTTgg                 >  2:49800/1‑220 (MQ=255)
                                                                                                                                                                                                                         gTCAGGCTTATGCG‑CAATTTGTCACCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCATAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGTTCTAAGCAGAAAGTCTATGGATCCAGTAGCGCATTCATCTACGGGCTCAACGTTATTATAGTAACTGGCGGCAATACTTTGCCGGTTTATCCATTGGTCGCCaa          >  1:343348/1‑227 (MQ=255)
                                                                                                                                                                                                                          tCAGGGTTTGGCG‑CAATTTGTCGCGAGGACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACGGGCTCTACGTTATTATAGTAACTGGCGGCAATACTTTGCCGGTTATACCTTTGGt                <  1:49800/220‑1 (MQ=255)
                                                                                                                                                                                                                             ggCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATACCTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTAc                                                                        >  1:60614/1‑161 (MQ=255)
                                                                                                                                                                                                                              gCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACg                                                                       <  2:60614/161‑1 (MQ=255)
                                                                                                                                                                                                                               cTTATGGG‑CAATTTGCCGCCAGCCCCGCCGAATGCTGTGGCGAAGGGGGAATC‑AAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACGGGCTCTACGTTATTATAGTAACTGGCGGCAATACTTTGCCGGTTATACCTTTGGTCGGCAATATTgg    <  2:43094/226‑1 (MQ=255)
                                                                                                                                                                                                                                 tttGCC‑CATTTTGCCGCCAGCTCCGCCGAATGGTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGGCGAATGGGCGTGATTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACGGGATCTACGTCATTATAGTAACTGGCGGCAATACTTTGCCGGTTATTCCTTTGGTCGGCAATATTGGCa  <  1:632708/225‑1 (MQ=255)
                                                                                                                                                                                                                                   tGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACGGGCTCTACGTTATTATAGTAACTGGCGGCAATACTTTGCCGGTTATACCTTTGGTCGGCAATATTg     >  2:444368/1‑222 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                               
TTACCTTCGCATTCATCACCCAAAATGTATTTAGGTACATGCTTTTTTGACCTTTCTTCCTCATGGGGAATTGATGACCGTGATGACCTGCTGCGCACAATTCACCGAATGATCGATAACGGTCATGCCGCTCG‑GTTGGCAGGGTTTTATCAC‑CGCTGGTTTCGTTATTCGCCATGTGAATGGCGTGACTATCTTGCTGAAC‑TAAATGAGCAAGGTCAGGCTTATGCG‑CAATTTGTCGCCAGCACCGCCGAATGCTGTGGCGAAGGGGGAATCAAAGCCTGGGATTATGTCCGGATGGGTTTTTTAAGCCGAATGGGCGTACTTAATAACTGGTTGAGCGAAGAAGAAAGTCTATGGATCCAGTCGCGCATTCATCTACGGGCTCTACGTTATTATAGTAACTGGCGGCAATACTTTGCCGGTTATACCTTTGGTCGGCAATATTGGCA  >  NC_000913/676780‑677228

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑