Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 266,024 A→G 14.3% S22P (TCT→CCT)  ykfB ← CP4‑6 prophage, uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913266,0240AG14.3% 169.4 / 2.1 74S22P (TCT→CCT) ykfBCP4‑6 prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (52/9);  new base G (2/10);  total (55/19)
Fisher's exact test for biased strand distribution p-value = 7.41e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGACGGAATAATGTACTGCTGGTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCG  >  NC_000913/265789‑266252
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                      
ccggACGGAATTATGTACTGCTGGTGATCCATCACTACGGTTATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCaga                                                                                                                                                                                                                                        <  2:132508/234‑1 (MQ=255)
   gACGGAATAATGTACTGCTGGTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTTCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑AAAGAAGGGGCAtt                                                                                                                                                                                                                                  >  1:402871/1‑239 (MQ=255)
                  tgctgcTGATCCATCACTATGGAGATCGTAGTATCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGACACCATTCTGGTCCGTGTTCCCTGTGCAGCTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCTGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGACGAGGCCTA‑AAAGCTCGC‑GAGcct                                                                                                                                                                                                                 >  1:25106/1‑237 (MQ=255)
                  tgctgGTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑CTAGCAGAGGCATT‑AAAGCTCGC‑GAGCAg                                                                                                                                                                                                                 >  2:24552/1‑239 (MQ=255)
                  tgctgGTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GGGCAg                                                                                                                                                                                                                 >  2:8889/1‑239 (MQ=255)
                  tgctgGTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GGGCAg                                                                                                                                                                                                                 >  2:244637/1‑239 (MQ=255)
                  tgctgGTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAg                                                                                                                                                                                                                 >  2:195975/1‑239 (MQ=255)
                  tgctgGTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAg                                                                                                                                                                                                                 >  2:388745/1‑239 (MQ=255)
                  tgctgGTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAg                                                                                                                                                                                                                 >  2:620801/1‑239 (MQ=255)
                  tgctgGTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAAAGGCATT‑AAAGCTCCC‑GAGCAg                                                                                                                                                                                                                 >  2:726321/1‑239 (MQ=255)
                  tgctgGTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGATTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCCG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAg                                                                                                                                                                                                                 >  2:406060/1‑239 (MQ=255)
                  tgctgGTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGGAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAACTACC‑GAGCAg                                                                                                                                                                                                                 >  1:839202/1‑239 (MQ=255)
                  tgctgGTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGCAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGCCCGTAACCCTG‑CTGCCACAACCCCGG‑G‑ATACCAGGGGACTA‑AAAGCtaccaacccg                                                                                                                                                                                                                  >  1:108011/1‑230 (MQ=255)
                  tgctgGCGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCGGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ACAGCAATGGAATT‑AAAGCTCCC‑GAGccg                                                                                                                                                                                                                 >  1:91587/1‑237 (MQ=255)
                       gTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGCCCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCCCCC‑GAGCAGAaca                                                                                                                                                                                                             >  1:4540/1‑238 (MQ=255)
                       gTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGAGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGTCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCGGACCTTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCACAGGCATT‑AAAgtcagcggaaggaaaa                                                                                                                                                                                                              >  1:785261/1‑223 (MQ=255)
                            ccATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAAGGGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GGGCAGTACACCCg                                                                                                                                                                                                         >  2:182252/1‑238 (MQ=255)
                              aTCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTTATACAAACCATTCTGGACCGGGTTGCCTGTGCAGTTGAT‑GGTAAAC‑ATCTTTCCGCTGGCTTCGGT‑CACCGTGC‑ATTCCGTATTGCCCTAGCCCTAGCCCTC‑CTCCCCGCACCCCGG‑G‑ATAGTAAAGGTAAT‑AAAGccacgggccatgaaacccgcccc                                                                                                                                                                                                      >  1:383220/1‑217 (MQ=25)
                                      ggTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAAcc                                                                                                                                                                                             >  2:624975/1‑239 (MQ=255)
                                             gtagtGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAg                                                                                                                                                                                      >  2:154055/1‑239 (MQ=255)
                                                gtgtCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GGGCAGTACACCCGCCAGCATAAACCGACTTAGGGa                                                                                                                                                                                   >  2:114714/1‑239 (MQ=255)
                                                 tgtCGTCATCGTCGGAACTGACCACTTGGTCACGGGCAA‑GGGTAAGAAAGACTGAAGCCTGATAGAACCCATTCTGGTCCGGGTTCCCTGTGCAGTTAAT‑GGTAAAC‑GTCTTCCCGCTGGCTTCGGT‑CACGCTGT‑AGTCCGTATTCCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACCCCCGCCAGCATAAACCGACTTAGTGa                                                                                                                                                                                   <  2:165725/238‑1 (MQ=255)
                                                   tcgtcCTCGTGCGAACTTCCCAGTTTGCCAGCGGCAA‑TGGTAAGAAAGACTTAATGTTGATAGAAACCATTCGGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTCCCGGTGGCTGCGGT‑CACGCTGT‑ATTCCGTATTGCCCTAACCGTACCCCTG‑CTGCCAGACCCCCGG‑G‑ATAGCAGAGGCATT‑AAACCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAg                                                                                                                                                                                <  1:381617/239‑1 (MQ=255)
                                                      tcATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGCCCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCCCCAGCCTAACCCGACTAAGGGaaaaaat                                                                                                                                                                             >  1:818098/1‑235 (MQ=255)
                                                        aTCGAGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTC‑CTGCCAGACCCCCGG‑G‑ATAGCAGAGGCATT‑AAAGC‑CGCAGGGCAGTAAACCCGCAACAGAAAACCGACTAAGGAAAAGCAtt                                                                                                                                                                            >  1:367393/1‑238 (MQ=255)
                                                             gCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GACTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGACCGGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ACAGCAGAGGCATA‑AAGGCTCGC‑ACGCGGTCCACCGCCCGGAAAAAACCGGATTATAGAAAGAAtggcaat                                                                                                                                                                       >  1:517999/1‑232 (MQ=255)
                                                                        cttttGTCTTGGGTAA‑GGGTAAGAAAGACTGAGTGGTGATAGAATACATGGTGGTCCGGGTT‑CCTGTGGTGGTTATAGGTAATC‑GTCTGTCCGCTGTCTGCGAT‑CCCG‑TGTAATTCCGCATGGCCATGCCCTTAGCCCTT‑CTGGCAGACCCCCGG‑A‑TTAGCGGGGCAATT‑AAGGCCCCG‑GAGCAGTAACTCCTCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCtt                                                                                                                                                           <  1:308939/236‑1 (MQ=255)
                                                                               cATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCTTATTGCCCTGACCGTAACCCTC‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCCGTACACCCGCCAGCATAAACCGACTTCGTGAAAGTATTgcacttctgttcccccttgccgtg                                                                                                                                                    >  1:892881/1‑216 (MQ=255)
                                                                               aaaggtacgggggAAAAAAAGAGTGAGGGCTTATAAAAACCATGCTGGGCCGGGTCCCCTGTGGCTTGTAG‑GGAAAAC‑GTCTTCCCGCGGGCTGCGGT‑CACGCTGT‑ATCCCGTATTGCCCGGCCCGAACCCCGG‑CTGCAAGACCCCCGG‑G‑ATAGAAGAGGTATT‑AAAGCGCGC‑GAGCAGTCCACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCgt                                                                                                                                                     <  1:404427/229‑1 (MQ=255)
                                                                                                    aCTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAg                                                                                                                                >  2:657707/1‑239 (MQ=255)
                                                                                                    aCTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGCATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTc                                                                                                                                           >  2:404427/1‑228 (MQ=255)
                                                                                                    aCTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGA‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAACTCCC‑GAACAGTACCACCGCCAGAATAAACCGACTTAGTAAAAGGTGTGTCACTCTTGTCTCCTTTGTAGTTTGTTTTTTCCTGATGGTcgg                                                                                                                                >  1:713515/1‑237 (MQ=255)
                                                                                                    aCTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑AACGCTTT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCACGC‑GAGCAGTCCACCCGCCCGACTAAACCAACTTGGGGAAAGAATTTGAATTCTTACCCCCTTTGTCGGTGTTTTATTCCtaaatggaag                                                                                                                                >  1:389842/1‑230 (MQ=255)
                                                                                                     cTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCCCCAGCATAAACCGCCTTAGTGAAAGAATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCTTGATTGTAAgg                                                                                                                               >  1:107842/1‑239 (MQ=255)
                                                                                                                         ttCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGCCTTAGTGAAAGAATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCCGTacc                                                                                                            >  2:148295/1‑236 (MQ=255)
                                                                                                                         ttCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGCATTGTCATTCTTGTCTCCTTTTTCGTTGTTTTATTCCTGATTGTCAGGGTTCGCGGGCATCGGTAGc                                                                                                            >  2:767922/1‑238 (MQ=255)
                                                                                                                               tCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGCGGGTATCAGTAGCTGcccc                                                                                                      >  2:469119/1‑238 (MQ=255)
                                                                                                                                               gcagtggttgctaacgttttttCCGCTGGCTTCGGC‑CACGCTGT‑ATCCCGTTTGCCCCGGCCCGTACCCCTG‑CGCCCAGAACCCCGG‑A‑ATAGCAGGGCTATA‑AAACCTCCC‑GCGCAGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCa                                                                                       <  1:831860/220‑1 (MQ=255)
                                                                                                                                                 ttGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTCGCTGCCCCATCAGTTTGCCTTCATggg                                                                                   >  2:831860/1‑239 (MQ=255)
                                                                                                                                                 ttGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTTCCCCATCAGTTTGCGGTCATgcg                                                                                   >  2:126336/1‑237 (MQ=255)
                                                                                                                                                  gttt‑GGTTAACGGTCTTCCCGTGGGCTGCGGC‑CACGCTTT‑TTGCCGTTTTGCCCAGCCCG‑AACCCTGTCTGCAAGACCCCCGG‑G‑ATAGCAGAGCCTTT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATggg                                                                                   <  1:726321/236‑1 (MQ=255)
                                                                                                                                                   gAT‑GGAAACC‑GTTTTGCCGTTGGCTTCGGT‑CACGCTGT‑ATTCCGTATGGCCCGGACCGTAACCCTG‑CTGCCAAACCCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCCCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCg                                                                                 <  2:402871/239‑1 (MQ=255)
                                                                                                                                                       tttgtaaccgttttCCCGTGGCCTTGGTCCACTCTGT‑TTCCCGTTTTGCCTAGCCCGGACCCGTG‑CCGCCAGACCCCTGG‑C‑AAAGCGGGGGTAAT‑AAAGCGCGC‑GAGAATTCCCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGt                                                                                <  1:195975/228‑1 (MQ=255)
                                                                                                                                                       ggTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTTTCAGGGTTCGA‑GGTATCATTAGCTGCCCCATCATTTTGCCGTCATGGGGGTACt                                                                             >  1:387981/1‑239 (MQ=255)
                                                                                                                                                       ggTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGt                                                                                >  2:718136/1‑237 (MQ=255)
                                                                                                                                                            cctgTCTTTCCGCTTGCTTGGTT‑CCCGTTGT‑ACTCCGTATCCCCCTGACCTTACCCCTG‑CTGCAAACACCCGGA‑A‑AAAGCAGAGGAATT‑AAACCTCTC‑GAGCAGTACCCCCGCCAGCATAAACCGACTTATTGAAATTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCaa                                                                          <  1:862870/236‑1 (MQ=255)
                                                                                                                                                                         tgggctcgggccACGTTTT‑ATCCCTTATGCCCCTGACCGAACCCCTG‑GCGCAAAACCCCCGG‑G‑ATAGCAGGGGCATT‑AAAGCCCGG‑GAGCAGTCACCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCtt                                                              <  2:818098/228‑1 (MQ=255)
                                                                                                                                                                         tGGCTTCGGT‑CCTGCGTT‑ATTCCGTTTGGCCCGACCCTTACCCCTG‑CTGCCA‑AACCCCGG‑GAATAGCAGGGGCTTT‑AAGGCTGCA‑GAGAAGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCTCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCtt                                                              <  2:25106/238‑1 (MQ=255)
                                                                                                                                                                          ctgttttcggtccctgttttacccGTTTTGCCCTGCCCGTACCCCTG‑CTGCCAGACCCCCGG‑G‑AAAGCAGAGGCATT‑AAAGCCCGC‑GAGCAGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCtt                                                              <  1:718136/217‑1 (MQ=255)
                                                                                                                                                                           gCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCg                                                                                          >  2:626577/1‑208 (MQ=255)
                                                                                                                                                                                ggT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATAc                                                      >  2:740827/1‑239 (MQ=255)
                                                                                                                                                                                         tGT‑ATTCCGTATTGCCCGGACCGTACCCCTG‑TCGCCAGACCCCCGG‑G‑ATAGCAGGGGCATT‑AAAGCCCGC‑GAGCAGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGt                                              <  2:107842/239‑1 (MQ=255)
                                                                                                                                                                                                 cGTTTTGCCCTGCCCGAATCCCGC‑GTGCCAGACCCCCGG‑G‑ATAGCAGGGGCATT‑AAAGCTCGC‑GACCATTCCTCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCAcctcc                                       <  2:892881/239‑1 (MQ=255)
                                                                                                                                                                                                   tATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCCGTAGCTGCCCCATACGTTTTCCGTCCTGGGCGTACTCAAAGAATTTTTTTTTGGTTTACGGGTCCGTCCcctccgg                                     >  1:133884/1‑237 (MQ=255)
                                                                                                                                                                                                   tATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGCGGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCAcccccgg                                     >  1:646205/1‑237 (MQ=255)
                                                                                                                                                                                                      tGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCCTCAGTTTGCCGTCATGGGCGTACTCAAAGAATTTTTCTTTGGTAACCGGGTCGGTCACCTCCTggaa                                  >  1:26214/1‑237 (MQ=255)
                                                                                                                                                                                                      tGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTa                                  >  2:289026/1‑239 (MQ=255)
                                                                                                                                                                                                      tGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGACGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTa                                  >  2:861323/1‑239 (MQ=255)
                                                                                                                                                                                                        cccTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGTAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGCTTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTCCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTAt                                 >  2:457647/1‑238 (MQ=255)
                                                                                                                                                                                                          ccttccctgactctttgtcccAGACCCCGGG‑G‑ATGCCAGA‑GCATTAAAAGCTCGG‑GACGTGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTc                               <  2:270113/220‑1 (MQ=255)
                                                                                                                                                                                                           tGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCa                             >  2:735477/1‑239 (MQ=255)
                                                                                                                                                                                                             aCCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGt                           >  2:728601/1‑239 (MQ=255)
                                                                                                                                                                                                             aCCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGAATTTTTTTTTGGTATACGGGTCCTTCACCCCCTGGTTTTCCAGt                           >  1:43944/1‑239 (MQ=255)
                                                                                                                                                                                                                  aaCCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGCCATTGGCGTACTCAAAGTATTTTTCTTTGGCATACGGGTCCGTCCCCTCTTGGTATTCCCGTTTGAt                      >  2:396538/1‑239 (MQ=255)
                                                                                                                                                                                                                       tg‑ctgcCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGGCCGTCACCCCCTGGTTTTCCAGTTTTATGTTAt                 >  1:835521/1‑239 (MQ=255)
                                                                                                                                                                                                                       tg‑ctgcCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGAATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTa                  >  2:826093/1‑238 (MQ=255)
                                                                                                                                                                                                                        tgccgcacgAACCCGGT‑G‑AAAGCGGGGGAATA‑AAACCTCGC‑AAGCAGTCACCCCGGCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTAt                 <  2:133884/231‑1 (MQ=255)
                                                                                                                                                                                                                        tcggccaggAACCCGGG‑G‑ATAGCAGAGGCATT‑AAAGCGCGC‑GAGCAGTCCCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTAt                 <  2:334419/231‑1 (MQ=255)
                                                                                                                                                                                                                         tcctccccacccccGGTG‑ATGGCGGGGGCATT‑AAAGCCTGG‑GGCCAGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTAt                 <  1:154055/229‑1 (MQ=255)
                                                                                                                                                                                                                          ctgcCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTACCCCCCTGTTATTCCAGTTTGATGTTAtgg               >  1:466381/1‑237 (MQ=255)
                                                                                                                                                                                                                                    cccGG‑G‑ATAGCAGGGGCATA‑AAAGCGCGA‑GAGAAGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATAcaca     <  2:108011/239‑1 (MQ=255)
                                                                                                                                                                                                                                      cGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATTGGCAATACACAg    >  2:713787/1‑238 (MQ=255)
                                                                                                                                                                                                                                       gg‑g‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACAAAgc   >  2:812778/1‑238 (MQ=255)
                                                                                                                                                                                                                                       gg‑g‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATAAACAgc   >  2:343788/1‑238 (MQ=255)
                                                                                                                                                                                                                                          g‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAGGCATTTTTCTTTGGAATACGGGTCCGTCCCCTCCTGGTATTCCAGTTTGATGTTATCGGCATTaaacgagg  >  2:432527/1‑233 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                      
CTGGACGGAATAATGTACTGCTGGTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAA‑GGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGAT‑GGTAAAC‑GTCTTTCCGCTGGCTTCGGT‑CACGCTGT‑ATTCCGTATTGCCCTGACCGTAACCCTG‑CTGCCAGAACCCCGG‑G‑ATAGCAGAGGCATT‑AAAGCTCGC‑GAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCG  >  NC_000913/265789‑266252

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