Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,093,736 T→G 23.3% V90G (GTG→GGG)  hisB → histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,093,7360TG23.3% 113.8 / 5.0 64V90G (GTG→GGG) hisBhistidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (36/10);  new base G (10/6);  total (48/16)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCTGATTAGCGAACCGCCGAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAA‑CCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCT  >  NC_000913/2093505‑2093974
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              
cccTGATTAGCGACCGGGCGAGTGATGTGCAGGTGGACCGTTTTGATAAACGCGCCTTTGAA‑CCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGAt                                                                                                                                                                                                                                          <  1:386632/238‑1 (MQ=255)
     aTTAGCGAACCGCCGAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAA‑CCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCCCAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACCGTTTGATTAAGGGCTGATTTGTc                                                                                                                                                                                                                                     >  1:899924/1‑238 (MQ=255)
     aTTAGCGAACCGCCGAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAA‑CCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCCCAACCTGATGATGCAGATCTTCACCTCCCAA‑GGCGTACAGTTTGATTAAGGGCTGATTTTTc                                                                                                                                                                                                                                     >  1:890641/1‑238 (MQ=255)
     aTTAGCGAACCGCCGAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAA‑CCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGCGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGACCTTCACCTCGCAA‑GGCGTACAGTTTGCAGAACTGCTGATTTGTcc                                                                                                                                                                                                                                    >  1:704564/1‑239 (MQ=255)
          cGAACCGCCGAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAA‑CCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCCCAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGGGCTGATTTTGCCGccc                                                                                                                                                                                                                                >  2:258114/1‑236 (MQ=255)
          cGAACCGCCGAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAA‑CCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGCTTCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGGGCTGATTTGTCCGCAc                                                                                                                                                                                                                                >  1:1038548/1‑238 (MQ=255)
            aaCCGCCGAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAA‑CCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGACTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGGGCGGATTTGTCCGCGccc                                                                                                                                                                                                                              >  1:126556/1‑237 (MQ=255)
               cgccgAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAA‑CCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCCCA‑GGCGGATTTCGATGGCCCGCCCACCCTGATGTTGCAGATCTTCACCTCGCAA‑GGCTTACAGTTTGATGAAGGGCTGGTTTGTCCGCACCTGcc                                                                                                                                                                                                                           >  1:151797/1‑238 (MQ=255)
               cgccgAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAA‑CCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGATGATTTTTCCGCACCTGcc                                                                                                                                                                                                                           >  1:927401/1‑238 (MQ=255)
               cgccgAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAA‑CCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCCCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGAGCTGATTTGTCCGCACCTGcc                                                                                                                                                                                                                           >  1:754464/1‑238 (MQ=255)
                         ttttCAGGTGGACCGTTTTGATAAACTCGCCTTTGAA‑CCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGt                                                                                                                                                                                                                >  1:796380/1‑239 (MQ=255)
                             ggggtggccttttttGAAAACCCCG‑CTTTGAACCCGGGGGTGACCCCGGAATGGCAGAAGCGACAAAAAGGGGGC‑TAAAAGCTGGGGATGATCACTAATGAGGTGGGTCTGGGAAAAAAAAGTTTCCCACA‑GGCGGATTTCGAGGGCCCGCACAACCTGATGATGAAGATCTCCACCCCGCAA‑GGCGTACAGTTTGATGAAGTGCTGTTTTGTCCGCACCGGCCCGCCGATGAGTGCGa                                                                                                                                                                                                            <  2:880110/228‑1 (MQ=255)
                                   gACCGTTTTGATAAACTCGCCTTTGAA‑CCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGGGCTGATTTGTCCGCACCTGCCCCCCGATGAGGGCGACTGCCg                                                                                                                                                                                                      >  2:810774/1‑239 (MQ=255)
                                        ttttGATAAACTCGCCTTTGAA‑CCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAg                                                                                                                                                                                                  >  1:328238/1‑238 (MQ=255)
                                                aaCTCGCCTTTGAA‑CCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCCCAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACCGTTTGCTGAAGGGCGGATTTGTCCCCCCCCGCCCCCCGAGTAGGGGGACTGCCGAAACCCGaaaaaa                                                                                                                                                                                         >  2:950137/1‑236 (MQ=255)
                                                  cTCGCCTTTGAA‑CCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCCCAA‑GGCGTAACGGTTGATGAAGTGCTGATTTGTCCGCCCCTGCCCGCCGCTGAGGGCGACTGCAGGAAACCGaaaagaaa                                                                                                                                                                                       >  2:602954/1‑234 (MQ=255)
                                                       cTTTGAA‑CCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACtgg                                                                                                                                                                                  >  1:603309/1‑239 (MQ=255)
                                                              cccGGGCGTAACCCCGGAACTGCGAAAGCTGAAAAAGGCGGGT‑AAAAAGCTGGTAATGATCATTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACG‑GGCGGATTCCGATGGCCCGCACACCCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGaa                                                                                                                                                                             <  1:357124/237‑1 (MQ=255)
                                                              acgggggGTGACCCGGGACGTGCTGAAGC‑GGAAAAAGGGGGCAGAAGAGCGGGTGATGACCAATAAGCGGGATGGTCTTGAAACAAAAATTTCCCCAAA‑GGCGAATTCGGATGGCCGGCACACCGTGATGATGCAGATCTTC‑CCCGGCAAGGGGGTACAGTTTGATGAAGGGCTGTTTTGTCCGCCCCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGaa                                                                                                                                                                             <  2:867741/235‑1 (MQ=255)
                                                                  ggCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGCTTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTa                                                                                                                                                                        >  1:377543/1‑239 (MQ=255)
                                                                  ggCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTa                                                                                                                                                                        >  1:927538/1‑239 (MQ=255)
                                                                  ggCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGACTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGCTGAAGGGCTGATTTGTCCGCACCTGACCGCCGCTGAGTGCGACGGCCGTAAGCCGAAAGGAAAACTGGTGGAAcagta                                                                                                                                                                        >  1:519446/1‑235 (MQ=255)
                                                                       gATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTg                                                                                                                                                                    >  1:292547/1‑238 (MQ=255)
                                                                                    ctgAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGCTGATGCAGATCTTCCCCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGCTCAGTGCGACTGCCGTAAGCCGAAAGCAAAACTGGGGGAACGTTATCTGGCGGAGCAAACgca                                                                                                                                                      >  2:382430/1‑237 (MQ=255)
                                                                                         gCTGCAAAAAGCGGGT‑TAAAAGCGGGGTATGATCACTAATCAGGATGGTCTGGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCGCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCCCCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGAt                                                                                                                                                  <  1:950137/238‑1 (MQ=255)
                                                                                             cAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGCGGAACGTTATCTGGCCGAGCAAGCGATGGGTcgc                                                                                                                                               >  1:583886/1‑237 (MQ=255)
                                                                                                aaaGCGGGC‑TAAAAGCGGGTGGTGAACAATAACCAGGATGGTCTTGAAACAAAAAGTTCCCCAAG‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTaa                                                                                                                                          <  1:382430/239‑1 (MQ=255)
                                                                                                   agacagttaaaaGTTGTTGATGATCACTAACCGGGGT‑GTCTTGGAAAAAAAAGTTTTCCACAGGGCGGATTGCGTGGGCCCGCAAACCGGGATGTGGAAGTTCTTCCCCCCGAAG‑GGCGTACGTTTGGATGAAGTGCTGATTTGCCCGCCCCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTaa                                                                                                                                          <  2:196591/227‑1 (MQ=255)
                                                                                                                  ggTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGGAAAACTGGTGGAACGTTATCTGGCTGAGCCAGGCAGGGAACCGGCCAAAAGTTAtggtgatgggg                                                                                                                         >  2:522222/1‑230 (MQ=255)
                                                                                                                   gTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGGGCGGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGGAAAACTGTGGGAACCGTTCCTGGCTGGGAAGGGGGGGGAAGGGCCTAACACTTTAGTGATTGGCg                                                                                                                         >  2:914027/1‑238 (MQ=255)
                                                                                                                    tgatgaTCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTAtgtg                                                                                                                                >  1:196591/1‑230 (MQ=255)
                                                                                                                        gatCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGGGGAACGTTATCTGGCTGAGCAAGCGATGGAACCGGCTAACAGTTAttgg                                                                                                                                >  2:740663/1‑223 (MQ=255)
                                                                                                                         atCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATcgcgc                                                                                                                  >  1:17217/1‑239 (MQ=255)
                                                                                                                             cTAATCAGGATGGTCTTGGCACACAACGTTTCACACA‑GGCGGATTTAGATGGCCCGCGCCACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGGCAGCCGAAAGAAAAACTGGCGGGCGGTTATCTGCCTGCGCAAACGATGGCGCGCGCTCACGATTATGAGATTGGGGATcgcgcaccc                                                                                                              >  2:120718/1‑235 (MQ=255)
                                                                                                                                aTCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTAGTGGATTGGCGATGCGGcgaccgac                                                                                                           >  1:1033833/1‑239 (MQ=255)
                                                                                                                                          tCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGGAAAACTGGGGGAACGTTATCTGGCTGAGCAAGGGAAGGATCGCGCTACAAGTATTGTGTTTGGGGATCGCGCGACCGAAATTCAACTgg                                                                                                 >  2:792285/1‑239 (MQ=255)
                                                                                                                                            ttGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGGGATTGGCGATCGCGCGACCGACATTCAACTGGCg                                                                                               >  1:518261/1‑239 (MQ=255)
                                                                                                                                            ttGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGCGGAACGTTATCTGGCTGGGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGCCATTCAACTGGCg                                                                                               >  1:608361/1‑239 (MQ=255)
                                                                                                                                            ttGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGAACAGTTTGATGAAGTGCTGATTTGTCCGCACCCGCCCGCCGCTGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGGGGAACGTTATCTGGCTGAGCAAGCGATGGAACGCGCTAAACGTTATGGGATTGGAGATCCCGCGACCGACATTAAACAgggg                                                                                               >  2:182927/1‑237 (MQ=255)
                                                                                                                                               gAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTg                                                                                                                                                                                                                                               >  1:61212/1‑92 (MQ=255)
                                                                                                                                               gAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTg                                                                                                                                                                                                                                               <  2:61212/92‑1 (MQ=255)
                                                                                                                                                   acaAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGaaaa                                                                                          >  1:867300/1‑237 (MQ=255)
                                                                                                                                                    cagaGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGGGGTTGGGGGACGCGCCACCGCAATTAAACGGGGGGaaaaaag                                                                                       >  2:913915/4‑236 (MQ=255)
                                                                                                                                                    caAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAC‑GGCGTACAGTTTGATGAAAAGCTGATTTGTCCGCCCCTGCCCGCCGATTAGTGCGCATGCCGGAAGACGAAAGTAAAACTGGTGGAACGTTAGCTGGCGGAGCAAGCGATGGAACGCGCTAACGGTTAAGTGCTTAGCGATCGCGCGGCCCACAACCAAATCGCGGGAAACAt                                                                                       >  2:130162/1‑239 (MQ=255)
                                                                                                                                                    caAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACa                                                                                                                                        >  1:1038866/1‑190 (MQ=255)
                                                                                                                                                    caAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGGGGAACGTTATCTGGCTGAGCACGCGATGGAACCGCCAAACAATTATGGGATTGGGGTTGGCGGGGCCGCAATTCAAATGGCGGGAAAAAt                                                                                       >  2:740481/1‑239 (MQ=255)
                                                                                                                                                            cccAGA‑GGCGGTTTCGGATGGCCCCAAAAACCTGTTGATGAAAATCTTCACCCCGAAA‑GGCGTACATTTGGATGAAGTGCTGATTTGCCCGCCCCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTa                                                                               <  1:258114/239‑1 (MQ=255)
                                                                                                                                                            cccACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGAAGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGCAAAACATGGGCAtt                                                                                >  1:864551/1‑238 (MQ=255)
                                                                                                                                                                ca‑GGCGGTTTTCGATGGCCCGCACAGCCTGATTGTGCAGATCCTCCCCTCGCAA‑GGCGTACAGTTTGATGAATTTCTGATTTTCCCGCCCCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTgg                                                                           <  1:330176/239‑1 (MQ=255)
                                                                                                                                                                   gggggttttGGAGGGCCCGCCAAACCTGATGAGGCAGCTCTTCACCTCGAAA‑GGCGTACATTTTGATGAAGGGCTGTTTTGTCCGCACCGGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGtt                                                                         <  1:336096/233‑1 (MQ=255)
                                                                                                                                                                                  cccGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGGGCGACCGACATTCAACCTGCGGAAAACAAGGGCATTCCGGGTTTACGccagaccccca                                                          >  2:400433/1‑229 (MQ=255)
                                                                                                                                                                                   ctccacaCCCTGATGATGCAGATCTTCCCCTGGCAG‑GGGGTCGTGTTGGATGAAGGTCTGTTTTGTCCGCCCCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACcgcg                                                          <  1:400433/235‑1 (MQ=255)
                                                                                                                                                                                   cggaacaACCGGAGGAGGCAGATCTTCACCCGGAAA‑GGCGTACAGTTTGATGAAGGGCTAATTGGCCCGCACCTGACCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACcgcg                                                          <  1:947396/234‑1 (MQ=255)
                                                                                                                                                                                           ccTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTAtg                                                                                                                                  >  1:61078/1‑158 (MQ=255)
                                                                                                                                                                                            gtgatgaGGAATATCTTCCACTCGCAG‑GGGGTTAAGTTGGAGTAGGGGGTGATTTGCCCGAACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGaa                                                <  1:522222/238‑1 (MQ=255)
                                                                                                                                                                                                     cAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGACCTGGCCAAt                                       >  1:109041/1‑239 (MQ=255)
                                                                                                                                                                                                     cAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAAt                                       >  1:1008830/1‑239 (MQ=255)
                                                                                                                                                                                                          cTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACGGGTTTACGCCACGACCCCGAAACCCTGAACGGGCCAATAATTg                                  >  2:854947/1‑239 (MQ=255)
                                                                                                                                                                                                                  cGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGACCa                           >  1:531058/1‑238 (MQ=255)
                                                                                                                                                                                                                  cGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGCCGATCGGGCGACGGACATTCAACTGGCGGAAAACATGGGCCTTACTGGTTTACGCTACGACCGCGAAACCCTGACCTGGCCAATGATggggcgtaac                          >  2:621668/1‑232 (MQ=255)
                                                                                                                                                                                                                   cgcggggtctACTGTTTGATGGTGGGTTGTTTGGCCCGCCCCTGCCCGCCGATGAGGGCGACTGCCGTAACCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCaa                          <  2:183308/230‑1 (MQ=255)
                                                                                                                                                                                                                   ccaagggCGTACGTTTGGATGAAGTTCTGATTTGCCCGCCCCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCaa                          <  1:257019/234‑1 (MQ=255)
                                                                                                                                                                                                                    cAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATAATTGGCGGACAAcc                        >  2:257019/1‑238 (MQ=255)
                                                                                                                                                                                                                                   tgatgaAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCCCGAAACCCTGAACTGCCCAATGATTGGGGAGAAAACCACCAGACGTGccc          >  2:799919/1‑239 (MQ=255)
                                                                                                                                                                                                                                           tGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGGGGGCAACTCACCAGACGGGACCGTTaaccc  >  2:865198/1‑235 (MQ=255)
                                                                                                                                                                                                                                           tGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCCACGACCGCGAAACCCTGAACTGGCCAATGATTGGGGCGCAACTACCCAGACGGAACCGTTACGc   >  2:737814/1‑238 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              
CCCTGATTAGCGAACCGCCGAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAA‑CCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGC‑TACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACA‑GGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAA‑GGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCT  >  NC_000913/2093505‑2093974

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: ‑