Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 71,877 A→G 20.8% E176G (GAG→GGG)  yabI → DedA family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091371,8770AG20.8% 169.3 / 4.9 87E176G (GAG→GGG) yabIDedA family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (8/57);  new base G (19/2);  total (28/59)
Fisher's exact test for biased strand distribution p-value = 6.70e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGT‑CGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTG  >  NC_000913/71650‑72107
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        
aaCTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGT‑CGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGTTCGATATTCCTGCCGGAATG‑CAGACG‑GGGCAGTTTAAATGGt                                                                                                                                                                                                                               >  1:1147127/1‑239 (MQ=255)
aaCTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGT‑CGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGGGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTTCTTCCTGCCAGGGGTTCTGGCGGGCGCGGGGGCCGAATTTCCTTCCGCAATG‑CAGGGC‑GGGGGGTTTAAAAGGt                                                                                                                                                                                                                               >  2:941939/1‑239 (MQ=255)
aaCGGAACATGGTTTGCAACAACCCAGCGTGTTC‑CCCTTTTGGTGGGGTCGTTTTTTTGGCCCGACGCGGCCGCTGGTGCAAA‑TGGTGGCGGGATTGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGt                                                                                                                                                                                                                               <  2:1048896/239‑1 (MQ=255)
     aaCATGCGTTGCATCAACACAGCATGTTCACCATTCTGGT‑CGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGGTCGATATTCCTGCCGGAATG‑CCGAGC‑GGGGGGTTTAAATTGTtgct                                                                                                                                                                                                                           >  1:1048896/1‑238 (MQ=255)
     aaCATGCGTTGCATCAACACAGCATGTTCACCATTCTGGT‑CGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGGGGCTAAATTTATTACGCCGCAAATTATCGGGTGCCTTCTGTGGCCCCCGTTTTACTTCCTGCCCGGGATTCTGGCGGGGGCGGGGATCGGTTTTCCTGCCGGAATG‑CCGAGG‑GGGGGGTTTTAATGGtggggg                                                                                                                                                                                                                          >  2:1374884/1‑234 (MQ=255)
               gCATCAACACAGCATGTTCACCATTCTGGT‑CGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTAAGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGGTTTACTTCCTGCCCGGGATTCTGCCGGGCGCGGGGGGTGTTATTCTTGCCGGAATG‑CCACGC‑GGTGGGGTTATATGGTtggctcggggcacaa                                                                                                                                                                                                                >  2:272294/1‑226 (MQ=255)
               cgcacAACACAGCTGGGTCACCATTCGGGT‑CGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑GGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAg                                                                                                                                                                                                               <  1:1284882/235‑1 (MQ=255)
                          gCATGTTCACCATTCTGGT‑CGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGGGGGTTTAAATGGTTGCTG‑CTGGCAACACC‑GGGGTTTTtg                                                                                                                                                                                                   >  1:943462/1‑239 (MQ=255)
                           cATGTTCACCATTCTGGT‑CGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGCATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGGCCCCGTTTTACTTTCTGCCAGGGATTCTGGCGGGCGCCGGGATGCATATTCCTTCCCGCATG‑CAGAGCTGGTGAGTTT‑AATGGTTGCTT‑CTGGCACCCCC‑GGGGttttttt                                                                                                                                                                                                  >  2:801691/1‑237 (MQ=255)
                           cATGTTCACCATTCTGGT‑CGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGGTCGATATTTCTGCCGGAATG‑CAGGGC‑GGGGGGTTTAAATGGGTGCTG‑CTGGCAACAGG‑GGTGGTTTtgt                                                                                                                                                                                                  >  1:1286554/1‑239 (MQ=255)
                           cATGTTCACCATTCTGGT‑CGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGGGAGTTTAAATGGTTGCTG‑CTGGCAAAAGC‑GGGGgtttttt                                                                                                                                                                                                  >  1:290828/1‑237 (MQ=255)
                             tGTTCACCATTCTGGT‑CGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGTTTCTTGCGGGCGCGGGGGTCGGTATTCCTGCCGGGATG‑CGGAGC‑GGTGGGTTTAAATGGTTTGTG‑CTGGCAAACGC‑GGGGTTTTTgggg                                                                                                                                                                                                >  2:147883/1‑239 (MQ=255)
                             atgttctccatgtggtgcggtggTTTTGTGGGCCCGTGTGGCCCGCTGGTGCCAG‑GGGGGGGGGGAATGGGGGTGCTGGCGGGGGCTAAATTTATTCATCCGAATATTACCGGCTGCCTGGGTGGGCCGCAGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGGTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTtgt                                                                                                                                                                                                  <  2:883151/221‑1 (MQ=255)
                                ttctcctttgtgggtgggctttttGTTGGCCCGACTGTTCCGCGGGTGCCAA‑TGGGGGCGGTAAGGCTGGATCTGCCGGTGGCTAAATTTATTGCGCGGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGt                                                                                                                                                                                              <  2:1261614/219‑1 (MQ=255)
                                     aTTCTGGT‑CGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTTCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTtggtgg                                                                                                                                                                                        >  1:1310609/1‑239 (MQ=255)
                                       tCTGGT‑CGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGGTATTCCTGCCGGAATG‑CAGAGC‑GGGTGGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGGGTTTTTGTGGGTtggtggtt                                                                                                                                                                                      >  1:385574/1‑237 (MQ=255)
                                       tCTGGT‑CGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCCGAATG‑CAGAGC‑GGGGGGTTTAAATGGGTGCTG‑CTGGCCACCGG‑GGTGGTTTTGTGGGGTGGtggcg                                                                                                                                                                                      >  1:1069099/1‑238 (MQ=255)
                                       agaagacaggtcgTTTTGTTGGCCCGACGGGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCTGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGtggc                                                                                                                                                                                       <  1:820312/231‑1 (MQ=255)
                                         tGGT‑CGGTCGTTTGGTTGGCCCGACGCGTCCGCTGGTGCCAA‑GGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGtggctgg                                                                                                                                                                                    <  1:941939/239‑1 (MQ=255)
                                          ggt‑cggtcGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGtggctggc                                                                                                                                                                                   <  1:22085/239‑1 (MQ=255)
                                           gt‑cggtcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGGATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTTCCAGGGGTTTCTGCGGGGGCGGCGGTCGGAATTCCCTCCCGAATGCCGGAGG‑GGTG‑GTTTAAAAGGGTGTTGTGTGACAACAGG‑GGGGttttttgggggggggggggggg                                                                                                                                                                                   >  2:187175/1‑222 (MQ=255)
                                           gt‑cggtcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGGATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGGTTCTGGCGGGCGCGGCGGTCGATATTCCTGCCGGAATG‑CGGAGC‑GGGGGGTGTAAATGGTTGCTG‑CTGGCGACACA‑GGGGGTTTTGTGGGTTGgtgggtgggt                                                                                                                                                                                  >  2:563728/1‑237 (MQ=255)
                                           gt‑cggtcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCCGGGGTTCTGGCGGGCGCGGCGGTCGATATTTCTGCCCGGATG‑CAGAGG‑GGGGGGGTTAAATGGTTGTTGTCTGGAAAAACG‑GGGGTTTTTGGGGGTggggggtgggg                                                                                                                                                                                   >  2:17716/1‑234 (MQ=255)
                                           gt‑cggtcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCCGGGATTCTGGCGGGCGCGGCGATCGGTATTCCTGCCGGAATG‑CAGAGC‑GGTGGGGTTAAATGGTTGCTG‑CTGGCCACAGC‑GGTGTTTTTGTGGGGTGGtggctg                                                                                                                                                                                     >  1:258108/1‑236 (MQ=255)
                                           gt‑cggtcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGGTTTACTTCCTGCCAGGGATTCTGGCGGGCGGGGGGGCCGGTATTTCTGCCCGAATG‑CAGGGG‑CGTGGGTTTAAATGGTTGTGT‑CTGGGCACAGC‑GGGGGTTTTGTGGGTGGGggggtgggg                                                                                                                                                                                  >  2:184497/1‑237 (MQ=255)
                                           gt‑cggtcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGCCTGCCTGCCGTGGCCGCCGGTTTACTTCCTGCCCGGGGTTCTGGTGGGCGCGGCGAGCGATATTCCCGCCGGAATG‑CACGCC‑GGTGGGTTTAAATGGTAGCTG‑CGGGCAAACGCGGGGGTTTTTGggggggggggcgggag                                                                                                                                                                                   >  2:508254/1‑227 (MQ=255)
                                            gtgggttgTTTTGGTGGCCCGACGCGCCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCGGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGtggctggc                                                                                                                                                                                   <  1:850959/235‑1 (MQ=255)
                                               gttttttttGTTGGCCCGGCG‑GGCCGGTTGGGCCAAGTGGGGGGGGGAATTGTGGATCGGCGGGTGGCTAAATTTATTCCGCCGAATATTAGCGGCTGCCTGCTGGGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCtgtg                                                                                                                                                                               <  2:1165383/234‑1 (MQ=255)
                                               ggtcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGGGGTTAAATGGTTGCTG‑CTGGCAACAGC‑GGGGTTTTTGTGGGTTGGTGGCTGGCtgtg                                                                                                                                                                               >  1:774223/1‑239 (MQ=255)
                                                 tcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGCATATTATCGGCTGCCTGCTGCGGCCGCCGTTTTACTTCCTGCCAGGGATTCTTGCGGGCGCGGCGATCCATATTTCTGCCGGAAGG‑CACAGG‑CGTGGGCTTCAATGGGTGCTG‑CGGCGAACAGC‑GGGGGTTTTTTGGGTtggggggagggcgggcg                                                                                                                                                                             >  2:1029839/1‑228 (MQ=255)
                                                 tcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGGGTTTAAATGGTTGCTG‑CTGGCAACACC‑GGGGTTTTTGTGGGTTGGGGGTTGGCtgtgtt                                                                                                                                                                             >  1:851803/1‑237 (MQ=255)
                                                 tcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGGTTTACTTCCCGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGGGC‑GGTGGGTTTAAATGGTTGTCG‑CGTGCAACAGC‑GGGGTGTTTTGGGGTGGGGGGGTGGGtgtggg                                                                                                                                                                             >  2:1313883/1‑237 (MQ=255)
                                                  ctgtttGTTGGCCCGACGCGTCCGCTGGTGCAAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCt                                                                                                                                                                             <  1:547731/235‑1 (MQ=255)
                                                    ttttttTGCCCCGGCGGGTCCTCGGGTGGCGG‑GGGGGGGTGGAAGTCTGGTTCGGCCGGTGGCTAAATTTATTCCGCCGAATTTTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGc                                                                                                                                                                          <  2:113052/239‑1 (MQ=255)
                                                    ttttGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGGGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGc                                                                                                                                                                          <  1:428044/239‑1 (MQ=255)
                                                    tttgGTTAGCCCGGCGGTTCCGCTGGTGCCAA‑TGGGGGCGGGAATGCGGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCGGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGc                                                                                                                                                                          <  1:681124/239‑1 (MQ=255)
                                                           ggCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGGCTGCTGTGGCCGCCGTTTTACTTCCTTCCAGGGATTCTGGCGGGCGCGGGGGTCGATATTCCTTCCGGGAGT‑CAGTGG‑GGTGGGGTTTAAAGGGTGCTG‑GTGGCGACAAC‑GGGGGTTTTGTGGGTGGGGGgctggttggggggggggctttt                                                                                                                                                                   >  2:820312/1‑225 (MQ=255)
                                                              ccccgggtcccGTTGGTGCCAT‑TGTTGGGGGAATTGTTGGATCTGCCGGTGGCTAAATTTATTACCCCGAATATTATCGGCTGCCTGCGGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGgcg                                                                                                                                                                <  2:774223/230‑1 (MQ=255)
                                                              ccGACGGGTCCGCTGGTGCGAG‑TGGTGGCGGAAATGCTGGTTCTGCCGGTGGCTAAATTTTTTACGCCGAATTTTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGgcg                                                                                                                                                                <  1:272294/239‑1 (MQ=255)
                                                                 aCGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGGATG‑CAGAGC‑GGGGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGGGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCgggttttgggcag                                                                                                                                                             >  1:715866/1‑228 (MQ=255)
                                                                      tCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CCGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGG‑GGGGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGGGGGTAttggggggggggg                                                                                                                                                        >  1:1344721/1‑227 (MQ=255)
                                                                        cgtgggggacAA‑GGGTGGGGGGAAGGGTTGATCGGCCGGTGGCTAAATTTATGACGCGGAATATTATCGGCTGCCTGCGGGGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTaaa                                                                                                                                                      <  2:943462/230‑1 (MQ=255)
                                                                               gCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTGTTTACCCCGAATATTATCGGCTGCCTGCTGGGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAGAGCGACt                                                                                                                                                <  1:560291/238‑1 (MQ=255)
                                                                               gCCAA‑GGGTGGAGGGAATGCTGGATTTGCCGGTGGCTAAATTTATTACGCCGAATTTTACCGGCTGCCTGCTGGGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACt                                                                                                                                                <  1:52778/238‑1 (MQ=255)
                                                                                     tggtggCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGGGTTTTTGGGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGGGcggcggggaagaggcgggcc                                                                                                                                            >  1:63455/1‑223 (MQ=255)
                                                                                        tggCGGGAAGTCTGGATCTGCCGGTGGCTAAATTTATTACCCCGAATATTATCGGCTGCCTGCGGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGttt                                                                                                                                        <  1:147883/238‑1 (MQ=255)
                                                                                        tggCGGAAATGCGGGCTCTGCCGGGGGCTAAATTTAGTCCGCCGAATTTTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGttt                                                                                                                                        <  1:59856/238‑1 (MQ=255)
                                                                                          gggggagtggtggtcttGCCGGGGGGTAAATTTTTTACGCCAAATTTTATCGGTGCCCTGCGGGGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGa                                                                                                                                      <  2:1099979/222‑1 (MQ=255)
                                                                                               aagtcTGGATCTGCCGGGGGATAACTTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCat                                                                                                                                 <  1:40326/234‑1 (MQ=255)
                                                                                                          ggCCGGTGGCAAAATTGTTTACTCCGAATATTAGCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCcg                                                                                                                      <  1:1029839/237‑1 (MQ=255)
                                                                                                                      tggtgtataagagacagattatCGGCTGCCGGCTGTGGCCGCCGTTTTACTTCCGGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGATGGCTGTGCTGGCGGTTATGGCGAAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTttgtt                                                                                                          <  1:534499/221‑1 (MQ=255)
                                                                                                                                  tgatttttttgggCTGCCTGCTGGGGCCCGCGTTTTTCTTCCTGCCAGGGATTCTGGCGGGGGCGGCGATCGATATTCCGGCCGGAATG‑CAGAGG‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCg                                                                                             <  2:918233/228‑1 (MQ=255)
                                                                                                                                  cgagtttttccggggtccggtggtgGCCGCGTTTTTACTTCCTGCCAGGGATTCGGGCGGGCGCGGCGATCGATCTTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCg                                                                                             <  2:1286554/217‑1 (MQ=255)
                                                                                                                                   gattagtcgcgcgtgcctgcTGTGGCGGCCGTTTTACTTCCTGCCAGGGATGCTGGCGGGCGCGGCGATCGATATTCCTGTCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGt                                                                                            <  1:677716/226‑1 (MQ=255)
                                                                                                                                    aaTTTTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTa                                                                                                                                                        <  1:313935/178‑1 (MQ=255)
                                                                                                                                             ggCTGCCGGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGgcg                                                                                                                                                                <  1:791000/161‑1 (MQ=255)
                                                                                                                                              gtgtggataagagccagccgTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGc                                                                                  <  1:1251781/222‑1 (MQ=255)
                                                                                                                                               tgtctgtgttggggcgcggTTTTTCTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGcc                                                                                 <  2:1276237/224‑1 (MQ=255)
                                                                                                                                                ttctcgctgGGGCCCCCTTTTTACTTCCTGCCAGGGTTTCTGGGGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGGGC‑GGGGAGTTTAAATGGTTGCTG‑CTGGCAACGGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCAt                                                                               <  2:980135/234‑1 (MQ=255)
                                                                                                                                                tgcctgcTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAAGG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCAt                                                                               <  1:681382/239‑1 (MQ=255)
                                                                                                                                                gtcctgtTTTGCCCGGCGTTTTGCTTCCTGGCCGGGATGCGGGCGGGCGGGGCGATCGGTATTCCGGCCGGAATG‑CAGGGC‑GGTGAGTTTAAAGGGTTGCTG‑CGGGCAACGGC‑GGTGTGTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCAt                                                                               <  2:1069099/237‑1 (MQ=255)
                                                                                                                                                  cctgcTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTGGCTG‑CTGGCAACAGC‑GGGGTTTTTGTGGGTTGGTGGCTGGCTGGGCTGGGGGGTTTGGCGGGCGGGGGAAACGACGGGCCGTTTTGGGCATTATTTTTTCCGCGGGCGGTTTTTGGGGGTgggggccggtgttttttttcagc                                                                             >  2:1051612/1‑217 (MQ=255)
                                                                                                                                                              gccgGTTTACGTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGgcgc                                                                 <  1:1051612/239‑1 (MQ=255)
                                                                                                                                                                 gTTTTACTTCCTGCCAGGGATTTCGGCGGGCCGGGCGATCGATATGCCGGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTT‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCtgg                                                              <  1:508254/239‑1 (MQ=255)
                                                                                                                                                                  ttttttttcttggcgggggTTTTGGCGTGCGCAGCGCTCGGTATTGCTGCCGGAATG‑AAGGGC‑GGTTGGTTTAAATGGGTGGTG‑CTGTCAACAGG‑GGTGTTTTTGGGGGTTGGGGGTGGGCTGTGGTGGCGGTTATGGCGCAGAGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCtggt                                                             <  2:468790/224‑1 (MQ=255)
                                                                                                                                                                  ttttGCTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGtt                                                                                                                                         <  1:774671/163‑1 (MQ=255)
                                                                                                                                                                    ttACTTCCTGCCAGGGATTCTGGCGGGGGCGGCGAGCGTAATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCtggtgg                                                           <  1:1282216/239‑1 (MQ=255)
                                                                                                                                                                                            gggcgCGGGGATCGATTTTCCTGCCGGAATG‑CAGGGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATg                                    <  2:725515/238‑1 (MQ=255)
                                                                                                                                                                                            gggcgCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATg                                    <  1:500308/238‑1 (MQ=255)
                                                                                                                                                                                            gggcgCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATg                                    <  1:1313883/238‑1 (MQ=255)
                                                                                                                                                                                            gggcgCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATg                                    <  1:1272449/238‑1 (MQ=255)
                                                                                                                                                                                            gggcgCCGCGATCGATATGCCTGCCGGAATT‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATg                                    <  1:801691/238‑1 (MQ=255)
                                                                                                                                                                                            gggccgggggttcgatTTTCCTGCCGGAATG‑CAGGGG‑GGTGAGTTTAAATGGTTGCTG‑CTGGAAACAGC‑GGTGTTTTTGGGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGc                                   <  2:715866/228‑1 (MQ=255)
                                                                                                                                                                                              gcgcGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGcc                                  <  1:17716/238‑1 (MQ=255)
                                                                                                                                                                                               cgcgGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCg                                 <  1:187175/238‑1 (MQ=255)
                                                                                                                                                                                                 cggggttcgtttttCCCGCCGGAATT‑GAGGGG‑GGTGGGTTTAAATGGTTGTGG‑CTGGCAACAGG‑GGTGTTTTTGGGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGgtg                              <  2:1147127/227‑1 (MQ=255)
                                                                                                                                                                                                    cgatcgatATTCCTGCCGGAATG‑CAGGGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGta                            <  1:563728/238‑1 (MQ=255)
                                                                                                                                                                                                                       actG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGCCCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTaaa         <  1:873966/236‑1 (MQ=255)
                                                                                                                                                                                                                       aaTG‑CAGAGC‑GGTGAGTTTAAATGTTTTCTG‑CTGGCAACAGC‑GGTGTTTTTGGGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTaaa         <  1:218261/238‑1 (MQ=255)
                                                                                                                                                                                                                            cAGAGG‑GGTTAGTTTAAATGTTGTCGG‑GTGGCAACGGG‑GGTGTTTTTGGGGGTTGGTGGCTGGCTGGGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGt    <  2:258108/239‑1 (MQ=255)
                                                                                                                                                                                                                             agagC‑GGTGAGTTTAAATGGTTGCTG‑GTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGt    <  1:1278485/238‑1 (MQ=255)
                                                                                                                                                                                                                             agagC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCTGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGt    <  1:425466/238‑1 (MQ=255)
                                                                                                                                                                                                                             agagC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGt    <  1:1374884/238‑1 (MQ=255)
                                                                                                                                                                                                                             agagC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGt    <  1:729813/238‑1 (MQ=255)
                                                                                                                                                                                                                               acagggTGAGTTTAAATGGTTGCTG‑CTGGCAACCGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGtt   <  1:1131831/234‑1 (MQ=255)
                                                                                                                                                                                                                                gggggTGATTTTTAATGGTTGCTG‑CTGGGCACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTTTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTg  <  2:290828/235‑1 (MQ=255)
                                                                                                                                                                                                                                gC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTg  <  1:84812/237‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        
AACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGT‑CGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAA‑TGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATG‑CAGAGC‑GGTGAGTTTAAATGGTTGCTG‑CTGGCAACAGC‑GGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTG  >  NC_000913/71650‑72107

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