Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 752,085 G→A 15.7% L29F (CTC→TTC)  ybgQ ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913752,0850GA15.7% 129.4 / 2.2 59L29F (CTC→TTC) ybgQputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (29/20);  new base A (3/7);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACACCGATCTATCTGTGGCAAGGACTGACGGATATCTGGTTTTAAACCAAATTTATCGACCAGTGAATCATTGATGCATGGAATGGTTTTGTCACCCTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTTAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTG  >  NC_000913/751847‑752303
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                          
attccccGATTTTTCTTGGGCAAGGGCTTAATGATATCTGGTGTTAAACCAAATTTATCGACCAGTGAAACATTGATGCATGGAATGGTTTTGTCGCCCTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCgag                                                                                                                                                                                                                            <  2:243020/234‑1 (MQ=255)
        atctatctGTGGCAAGGACGGACGGATATCTGGTTTTAAACCAAATTTATCGACCAGTGCATCATTTATGCATGGAATGGTTTTGTCACCCTTTTTTT‑GCCAGTTAATTTTTTGCCCGTTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATGTT‑CACCGCGAGCAATGCCTCCTTTTTCCTTT‑CATAAGGAAATGTCAATGCGGTCCGACAGTTATTTATCGGAACCATTCa                                                                                                                                                                                                                    >  1:886523/1‑239 (MQ=255)
          ctatctGGGGCAAGGACTGACGGATATCTGGTTTTAAACCAAATTTATCGACCAGTGAATCATTGATGCATGGAATGGTTTTGTCACCCTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGa                                                                                                                                                                                                                  <  1:879420/239‑1 (MQ=255)
               tgtgGCAAGGACTGACGGATATCTGGTTTTAAACCAAATTTATCGACCAGTGAATCATTGATGCATGGAATGGTTTTGTCACCCTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATTCGGTCGCGCATTGCTTTATCGAGAACATTCAGATTGaa                                                                                                                                                                                                             >  1:429906/1‑239 (MQ=255)
                gtgGCAAGGACTGACGGATATCTGGTTTTAAACCAAATTTATCGACCAGTGAATCATTGATGCATGGAATGGTTTTGTCACCCTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTACCGAGAACATTCAGATTGAAt                                                                                                                                                                                                            >  1:791171/1‑239 (MQ=255)
                  ggCAAGGACTGACGGATATCTGGTTTTAAACCAAATTTATCGACCAGTGAATCATTGATGCATGGAATGGTTTTGTCACCCTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAAAAAATTAAAATTAAATTc                                                                                                                                                                                                          >  2:53658/1‑239 (MQ=255)
                   gCAAGGACTGACGGATATCTGGTTTTAAACCAAATTTATCGACCAGTGAATCATTGATGCATGGAATGGTTTTGTCACCCTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCCCATTGATTTATCTAAAAAATTCAAATTGAAttag                                                                                                                                                                                                         >  2:391064/1‑237 (MQ=255)
                       ggACTGACGGATATCTGGTTTTAAACCAAATTTATCGACCAGTGAATCATTGATGCATGGAATGGTTTTGTCACCCTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGAcc                                                                                                                                                                                                      >  1:709917/1‑238 (MQ=255)
                                      tttttttAAACCAAATTTTTCGACCTGGGAATCTTTGTTGCGTGGAATGGTTTTTGCCCCTTTTTTTT‑CCCAGTTAATTTTTTCCCCATTACTGTTTTTTGTGTTAGTCCCCGC‑AACGTAAACAAAATTTT‑CCCGGGGAGCAATGACTCCTTTTTCCTTT‑AATAATAAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCa                                                                                                                                                                                      <  2:709917/236‑1 (MQ=255)
                                      tttttttAAACCAAATTTATCGGCCAGTTAACTATTGATTCATGGAATGTTTTTGTCCCCCTTTTTTT‑CCCAGTTAATTTTTTGCCCATTCCTGATTTTTTTTTTATTCCCCGC‑AACGCTAACAAAATATT‑CCCCGGGGTCAATGCCTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCa                                                                                                                                                                                      <  2:648067/236‑1 (MQ=255)
                                      tGGTTTTAACCCAAATTTATCGACCAGTGAATCATTGATGCATGGAATGTTTTTGTCACCTTTTTTTG‑GCCATTTAATTTTTTGCCAATTACTGATTTTGTTGTTATTCCCCGA‑ACAGGTAACAAAATATT‑CCCCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCa                                                                                                                                                                                      <  1:65635/239‑1 (MQ=255)
                                           ttAAACCAAATTTATCGACCAGTGAATCATTGATGCATGGAATGGTTTTGTCACCCTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCCCCATTAACCAAGGCATCGcc                                                                                                                                                                                 >  2:558211/1‑239 (MQ=255)
                                             aaaCCAAATTTATCGACCAGTGAATCATTGATGCATGGAATGGTTTTGTCACCCTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAgg                                                                                                                                                                                        >  1:419477/1‑230 (MQ=255)
                                                    aTTTATCGACCAGTGAATCATGGTGGCATGGAATGGTTTTGTCCCCCTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcagg                                                                                                                                                                        <  1:81818/239‑1 (MQ=255)
                                                      ttatgagacAGTGAATCATTGATGCATGGAATGGTTTTGTCCCCCTTTTTTT‑GGCAGTTAATTTTTTGCCATTTACTTATTTTGTTGTTTTGCACCGC‑AACGCTAAAAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AAAAATGAAATATCAATGCGGTCGCGCAGTGATTTGTCGAGAACATGCAAATGGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcaggca                                                                                                                                                                      <  1:215425/231‑1 (MQ=255)
                                                           gACCAGTGAATCATTGTTGCATGGTATGGTTTTGTCACCCTTTTTTT‑GCCATTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACCATAACAAAATATT‑CCCCGGGAGCATTGCCTCGTTTTGCCTTT‑AAGAATGAAATATCAATGCGGCCCCGCATTGATTTATCGAGAAAATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATa                                                                                                                                                                 <  1:14102/239‑1 (MQ=255)
                                                                 tGAATCATTGATGCATGGAATGGTTTTGTCACCCTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCACATTGATTTATCGAGAACATTCAGATTGAATTCCACCGCCATTGCGCAAGGCATCGACATTACCAGGCAGGAAACaaaaa                                                                                                                                                           >  2:807936/1‑238 (MQ=255)
                                                                           tttgaTGGATTGTTTTTGTCTCCTTTTTTTT‑CCCATTTATTTTTTTGCCTTTTATGTGTTTTGTTTTTATTACCCAC‑ACCTCTAAAAAAATATG‑CACCGGGGTCAATGACTCCTTTTTCTTTT‑AATAAGGAAATATCAATGCGGTCGCGTATTGATTTATCGAAAACATTCAGATTAAATCGGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATa                                                                                                                                                 <  2:726363/235‑1 (MQ=255)
                                                                           ctgattGGAATGTTTTTGTCCCCCTTTTTTT‑GCAAGTTATTTTTTTGCCAATTATTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCCCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATa                                                                                                                                                 <  1:229598/234‑1 (MQ=255)
                                                                           aTGCATGGAATGGTTTTGTCACCCTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATa                                                                                                                                                 >  1:69868/1‑239 (MQ=255)
                                                                           aTGCATGGAATGGTTTTGTCACCCTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATa                                                                                                                                                 >  1:809337/1‑239 (MQ=255)
                                                                                    tttgTTTTTTCTCCCTTTTTTT‑GCCATTTAATTTTTTGCCCATTCAGGTTTTTGTTGTTATCCCCCGC‑AACGCAAAAAAAATATT‑CACCGGGAGAAATGACTCTTTTTCCCTAT‑AATAATAAAATAACAATGCGGTCGCGCATTAATTTACCGAGAACATTCAGATTGAATACGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCAc                                                                                                                                        <  2:429906/236‑1 (MQ=255)
                                                                                     tGGTTTTTTCCCCCTTTTTTT‑GCCAGTTATTTTTTTGCCCATTACTGATTTGGTTGTTATTCCCCGA‑AACGTTAAAAAAATATT‑CACCGGGAGCAATGCCTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGAATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACg                                                                                                                                       <  1:391064/239‑1 (MQ=255)
                                                                                     tGGTTTTGTCCCCTTTTTTTT‑GCCAGTTAATTTTTTGCCAATTACTGATTTTGTTGTTATTCACCGA‑AACGCTAACAAAATATT‑CACCGGGATCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACg                                                                                                                                       <  1:803446/239‑1 (MQ=255)
                                                                                     tGGTTTTGTCCCCCTTTTTTT‑GCCAGTTAATTTTTTGCCATTTACTAATTTTGTTGTTATTCACCCA‑AACCCTAAAAAAAAATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAAGGAAATATCAATGCGGTCGCGCATTGATTTATCGAAAACATTAAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACg                                                                                                                                       <  1:558211/239‑1 (MQ=255)
                                                                                     tGGTTTTGTCCCCCTTTTTTT‑GCCAGTTAATTTTTTGCCAATTACTGATTTTTTTGTTATTCCCCGC‑AACGCTAACAAAATATT‑CCCCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACg                                                                                                                                       <  1:598699/239‑1 (MQ=255)
                                                                                      ttataggTCCCCCTTTTTTTGGCCATTTATTTTTTTCCCATTTTCTGATTTGTTTTTTTTTCACACC‑AACGCTAAAAAAATATT‑CACGGGGGACAATGACTCTTTTTTCCTTT‑AAAAAAGAAATATCAAGGCGGTCGCGCATTGATTTATCGAAAACATCCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACg                                                                                                                                       <  2:476835/233‑1 (MQ=255)
                                                                                       gTTTTGTCCCCTTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCCCCGC‑AACGCTAAAAAAATTTT‑CACCGGGAGCAATGCCTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACg                                                                                                                                       <  1:779655/237‑1 (MQ=255)
                                                                                             tCACCCTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAACTATT‑CACCGGGAGCATTGAATCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGAGCAAGGCATCGCCATGACCCGGCCGGGGACAAAAGCGATTCGATAAAAATTACCGGGGTCCAt                                                                                                                               >  1:410958/1‑239 (MQ=255)
                                                                                             tCACCCTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCa                                                                                                                                >  1:627751/1‑238 (MQ=255)
                                                                                             tCACCCTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGAAACAAAAGAGAGTCGAAAAAAATTccaggggcccaa                                                                                                                               >  2:274174/1‑228 (MQ=255)
                                                                                               tcctttttttt‑tCCTTTTTATTTTTTCCCCTTTGCATTTTTTTTTGTTTTTCCCCCA‑AACGTTAAAAAAATTTT‑AACCGGGAACAAGTCCCTCTTTTTCCTTA‑AATAATGAAATATCAATGCGGCCCACAATTGATTTACCGGGAATATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATAc                                                                                                                             <  2:791171/235‑1 (MQ=255)
                                                                                                      tttt‑GCCAGTTAATTTTTTGCCAATTACTGATTTTGTTGTTATTCACCGA‑AACGCTAACAAAATATT‑CCCCGGGAGAAATGCCTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGAATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATaa                                                                                                                      <  1:236609/239‑1 (MQ=255)
                                                                                                            ccAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAAAGTCGATAAATATTCACGGTGTCCATAc                                                                                                                             >  1:476835/1‑227 (MQ=255)
                                                                                                                ttAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATGTT‑CACCGGGCGCATTGACTCCTTTTTCCTTA‑AATAATGAAATATCAATGCGACCGCGCATTGATTTATCGAGACATTTCAGATTGAATTCGACCGCCATTGACCAAGGCATCGCAATGACCAGGCAGGCTACAAAAGAGAGTCGATAAGACTTCACGGTGCTAATACCTGATAAATATTTTAt                                                                                                             >  2:620334/1‑239 (MQ=255)
                                                                                                                       ttggCCAATTCCTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CCCCGGGAACAATGACTCCGTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAggc                                                                                                      <  1:97763/236‑1 (MQ=255)
                                                                                                                               ttACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCgatga                                                                                              >  1:179109/1‑239 (MQ=255)
                                                                                                                                tACTGATTTTTTTGTTATTCACCGC‑AACGAAAACAAAATATT‑CCCCGGGAGCAATGACTCCTTGTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTAACGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACg                                                                                                                                       <  1:474268/197‑1 (MQ=255)
                                                                                                                                      ttttttttttATTTCCCGCAAACTCAAAAAAAATTTG‑CCCCGGGAGAAATATCTCGTTTTCCCTTA‑AATATGAAAAAATAAACGGGGCCGCGTATTATTTTAGCGAGAAGCTACAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATgccgc                                                                                        <  2:809337/239‑1 (MQ=255)
                                                                                                                                        ttgttgTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGGGTCCATACCTGATAAATATTTTATGAAAAACGGCGATGATTCCCCCaa                                                                                     >  2:489498/1‑239 (MQ=255)
                                                                                                                                                aTTCCCCGC‑AACGCTAACAAAATATT‑CCCCGGGAGCAATGCCCCCTTTTACCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAAAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATAc                                                                             <  1:489498/239‑1 (MQ=255)
                                                                                                                                                aTTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAaataaa                                                                               >  1:400445/1‑236 (MQ=255)
                                                                                                                                                      cGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGGGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACttattt                                                                       >  1:353453/1‑239 (MQ=255)
                                                                                                                                                        caacctatAAAAAAATTTT‑C‑CCGGGGGGAAAGTATCTTTTTTCCTTTAAAAAATAAAATACAAATGCGGGCCCGAATGTATTTATGGAAAAAATCCAGTTAGAATTCAACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACttattta                                                                      <  2:179109/232‑1 (MQ=255)
                                                                                                                                                        c‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATACTACttatttat                                                                     >  1:347550/1‑239 (MQ=255)
                                                                                                                                                            cgtgaaaaaaaTATT‑CACCGGGAGCAATGCCTCTTTTTTCCTTT‑AATAATGAAATATAAATGCGGTGGCGGTTTGATTTATCGAAAACTTTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAAt                                                                  <  1:274174/232‑1 (MQ=255)
                                                                                                                                                               tAACAAAATATT‑CACCGGGAGCAATGCCTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTACCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATcc                                                                <  1:53658/238‑1 (MQ=255)
                                                                                                                                                                 aCAAAATATT‑CACCGGGAGCAATGAGTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGCCAGGTTGAAAAAGAGAGCCGATAAAAATTCCAGGTTTCAATACCTGATAAATATTTTTTGCAAGGCGGCAATAATGCCCCCAAAATAAAATTTATTTTAAATCCaa                                                              >  2:147174/1‑237 (MQ=255)
                                                                                                                                                                 aCAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAAAATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGAGGCGATGATGCCGCCAAAAAAACACtttttttaataccaac                                                             >  2:831480/1‑229 (MQ=255)
                                                                                                                                                                   aaaaTATTACACCGGGAACAATGACTCATTTTTCCTTA‑AATAATGAAAATTCAAGGCGGTCGCGATTTGATTTACCGAGAACATTCATATTGAACTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATAGTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCa                                                            <  1:620334/239‑1 (MQ=255)
                                                                                                                                                                            cACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGCCAGGATACAAAAGCGAGTCGATAAATATTGACGGTGTCCTTACCTGATAACTATTTTATGAAAGCCGGCGTGGATGCGCCCAAACAAATACTTTTTTTTAAACCAGCAccacgggag                                                   >  2:773454/1‑231 (MQ=255)
                                                                                                                                                                                        tGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGAGGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTTAGCGGtt                                        >  1:168916/1‑238 (MQ=255)
                                                                                                                                                                                             ccTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTTAACGGTTACTTc                                   >  1:869137/1‑238 (MQ=255)
                                                                                                                                                                                                ttttCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAACAGCGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTTAGCGGTTACTTAAcc                                >  1:397312/1‑238 (MQ=255)
                                                                                                                                                                                                   tCCTTT‑AAAAAAGAAATATGAGGGCGGTCGCGCTTTTATTTACGAAAAATATTAAGTTTATTTGCGCCCGCCATTGCGCAAGGCATCGCCATGACAAGGCAGAATAAAAAAGAGGGTCGATAAATATTCACGGTGCCCATACCTGAAAAATATTTTAGGAAAGGCGGCGATGATGCCGCAAAAAAAATCCTTATTTATAATCCAGCACGTAGGTTGCTTTAGCGGTTACTTCCCCTGcc                            <  2:869137/239‑1 (MQ=255)
                                                                                                                                                                                                        t‑aataaTGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTAt                                                                                                             >  1:114083/1‑153 (MQ=255)
                                                                                                                                                                                                        t‑aataaTGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAACGTCGGCGATGATGCCGCAAAAGTAATACTTATTTATAATCCAGCACGTAGGTTGCGTTAGCGGTTACTTCCCCTGCCGCGAc                       >  1:502066/1‑239 (MQ=255)
                                                                                                                                                                                                        t‑aataaTGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTTAGCGGTTACTTCACCTGCCGTGAc                       >  1:716318/1‑239 (MQ=255)
                                                                                                                                                                                                        t‑aataaTGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCACAAATAATACTTATTTATAATCCAGCCCGTAGGTTCCGTTACGGGTGACTTCACCTGCCGTGAc                       >  1:254211/1‑239 (MQ=255)
                                                                                                                                                                                                                              gTCGCGCATTGATTTATCGAGAATATTCAGCTTGAATTCGACCGCCATTGCGCAAGGCATCGCCATTACCAGGCAGGATACAAAAGAGAGTCGATAACTATTCACGGTGTCCATACCTGCTAAACATTTTATGAAAGCCGGCGATGATCCCGCCAAAATCATACTTATTTTTAATCTAGACCCTAGCTTGCGTTAGCGGTTACTTACCTTGCCGTCAGACCGACTGCattttgattttg  >  2:614654/1‑239 (MQ=255)
                                                                                                                                                                                                                                                 agaACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACttattta                                                                      >  1:136436/1‑152 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                          
ATACACCGATCTATCTGTGGCAAGGACTGACGGATATCTGGTTTTAAACCAAATTTATCGACCAGTGAATCATTGATGCATGGAATGGTTTTGTCACCCTTTTTTT‑GCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGC‑AACGCTAACAAAATATT‑CACCGGGAGCAATGACTCCTTTTTCCTTT‑AATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTTAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTG  >  NC_000913/751847‑752303

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