Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 5,439 A→C 25.7% H69P (CAC→CCC)  yaaX → DUF2502 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,4390AC25.7% 117.5 / 4.1 71H69P (CAC→CCC) yaaXDUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (29/20);  new base C (7/15);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 4.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGA‑TGCAATCTATCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGT  >  NC_000913/5201‑5653
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                         
aCCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGA‑TGCAATCTATCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCa                                                                                                                                                                                                                           <  2:298635/239‑1 (MQ=255)
     tAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGA‑TGCAATCTATCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTac                                                                                                                                                                                                                      >  1:899400/1‑239 (MQ=255)
      aaCGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGA‑TGCAATCTATCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTac                                                                                                                                                                                                                      >  1:768666/1‑238 (MQ=255)
           gCAATGATAAAAGGAGTAACCTGTGAAAAAGA‑TGCAATCTATCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCCCTGGCACCTACAC‑GGa                                                                                                                                                                                                                >  1:218043/1‑238 (MQ=255)
            cAATGATAAAAGGAGTAACCTGTGAAAAAGA‑TGCAATCTATCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAAATACATTATGAATGGGCA‑GGCAATCGGTGGACACTACAC‑GGAc                                                                                                                                                                                                               >  1:42268/1‑238 (MQ=255)
                         aGTAACCTGTGAAAAAGA‑TGCAATCTATCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGACTGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CAccgcc                                                                                                                                                                                                 >  1:912166/1‑238 (MQ=255)
                           tAACCTGTGAAAAAGA‑TGCAATCTATCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACGGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCGCGGCTGGTGGAACCAAAATTATGAATGGCGG‑GGCAATGCCTGGCCCCACAAC‑GGAACGCCGC‑CCCcgccgc                                                                                                                                                                                               >  2:554169/1‑238 (MQ=255)
                           tAACCTGTGAAAAAGA‑TGCAATCTATCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGCATGGCGA‑GGCAAGCGCTGGCCCCTCCAC‑GGACCCCCCC‑CCCccccgc                                                                                                                                                                                               >  2:459826/1‑238 (MQ=255)
                           tAACCTGTGAAAAAGA‑TGCAATCTATCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCGGGCCCCTACAC‑GGACCGCCGC‑CAccgcccc                                                                                                                                                                                               >  1:680497/1‑236 (MQ=255)
                           tAACCTGTGAAAAAGA‑TGCAATCTATCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGGGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACcgccgc                                                                                                                                                                                               >  1:1316165/1‑238 (MQ=255)
                           tAACCTGTGAAAAAGA‑TGCAATCTATCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATAGTGGCTATTAAGTGGATGGAGGTCACTGACCGG‑ACACCGCCTGGGGGGAAAAAAATTATGATTGGCCG‑GGCAAACGGCGGCACCTACAC‑GGACCGCCGC‑CACgcccgc                                                                                                                                                                                               >  2:983682/1‑238 (MQ=255)
                             accctggAAAAAGATTGCAATTTATCGGCCCCCCATTTCCCCTGGTTC‑GGGGCGCCCCCAGGGCAGCCAGGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATGGGGGATCGTGATAATCGTGGCTATTACGGGGATGGAGGTCACGGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCgcgc                                                                                                                                                                                             <  2:899400/233‑1 (MQ=255)
                              cggagatgtggaaaagagacagaTCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCgcgc                                                                                                                                                                                             <  1:539110/216‑1 (MQ=255)
                              cctttaaaaaaGA‑TGAAATCTTCCGT‑CCCTCCCTTTCCCTGGGTCGTGGCGCCTCCCATGGCAGACCAGGCTGGGGAAATTACGTTAGCCCAGTCAATAAAATTACAAATAGGCGTGCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGAAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGcca                                                                                                                                                                                           <  2:1209427/233‑1 (MQ=255)
                              ccTGTGAAAAAGA‑TGCAATCTATCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGAAAGGCGATCGTGGTAATCGTGGCTATTACTGGGATGGAGGTCACTGGCCCG‑ACCACGGCTGGCGG‑AACAAAACTATTGATGGTGAGGGGAAAACCTGGCACCCACAC‑GAACcccccc‑ccccccgccggc                                                                                                                                                                                            >  2:956760/1‑230 (MQ=255)
                                 ttaaaaaaGA‑TGCATTCTGTGCTGCCCCCACTTTCCCTGTGTC‑TGGCCCCTCCCAGGGCAGCAAAGGCTGCGGAAATTACGTTAGTCCCGGCAGGAAAATTACAGATAGGCGGTCGTGATAATCGTTGCTATTACTGGGATGGAGGTCACGGGGGCG‑ACCACGGCGGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGccacca                                                                                                                                                                                        <  2:695410/236‑1 (MQ=255)
                                  tgAAAAAGA‑TGCAATCTATCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACACAATTATGAATGGCGA‑GGCAATCACTGGCACCTCCAC‑GGACCCCCCC‑CACCCCCCCCcccccaa                                                                                                                                                                                       >  2:1000318/1‑238 (MQ=255)
                                         aa‑TGCAATTCTTCGTCCCCGCACTTCCCCTGGTTT‑GGGCCGCCCCCATGGCAGCCCGGGCTGGGAAAATTACGTTAGTCCCGTCAGAAAAATTACAGAGAGGCGATCGTGAAAACCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACACCATTATGAATGGCGA‑GGCAATCGCTGGCCCCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAg                                                                                                                                                                                <  2:768666/238‑1 (MQ=255)
                                            tGCAATCTATCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGc                                                                                                                                                                               >  1:271285/1‑238 (MQ=255)
                                             gCAATCTATCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATACGAAAGct                                                                                                                                                                              >  1:941051/1‑238 (MQ=255)
                                               aagcTATCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGctcc                                                                                                                                                                            <  1:806100/235‑1 (MQ=255)
                                                ctttttCGTTCTCCCTCTTTCCTTGGTTT‑TGGTCGCCCCAAGGGCACACCAGGCTGGGGAAATTCCGTTAGTCCGGTCAGTAAAATTACAGAAGGGCGATCGTGATAATCGTGGCTATTACTGGGAGGGAGGTCACTGGCGCG‑CCCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCCCCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGctcctc                                                                                                                                                                          <  2:1316165/234‑1 (MQ=255)
                                                         aCTCGCCTTTTCCCTGGTTC‑GGGGGCCCCCAATGGCAGCAAAGGCGGCGGAAATTCGGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTAGGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGatcat                                                                                                                                                                  <  1:1190131/238‑1 (MQ=255)
                                                          cTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCCA‑GGCAATCGCTGGACCCAACAC‑GGACCCCCGCACACGACCGCCCCCCCaaaaaaaacaccctcaatacaaaa                                                                                                                                                                 >  2:73833/1‑223 (MQ=255)
                                                          cTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGGGGAAACAACATTATGAAAGGCGT‑GGCAATCACGGGCCACTACAC‑GGACCCCCGC‑CACCACCGGACCCCCATAAGAAAgaacacgataataatca                                                                                                                                                                >  2:95736/1‑223 (MQ=255)
                                                           tCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAACTTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GCCAATCGCCTGCCCCTACAC‑GGACCGCCCC‑CACCCCCCCCGCCCCATAAAACAACCCCTCATCaacatacc                                                                                                                                                               >  2:698294/1‑236 (MQ=255)
                                                             gCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCAGAATCATCAcg                                                                                                                                                              >  1:754940/1‑238 (MQ=255)
                                                               aCTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CCCCGCCGCGCCCCCATAAGAAAGCTCCTCATGATCATCAcggcg                                                                                                                                                           >  1:341469/1‑239 (MQ=255)
                                                               aCTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGAGGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCTA‑AGGAATCGCTGGCACCTACAC‑GGCCCACCGC‑CCCCCCCCGGCCCCCCTAAAAAAACACCCCATCACCAACAcggcg                                                                                                                                                           >  2:1344023/1‑239 (MQ=255)
                                                                         gTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTAGTAATGGCGA‑GGCACACGCTGGCACCCACAC‑GGACCGCCCC‑CACCCCACCCCCACCCAAAAAAAGCCCCCCAtctccacaaccgccggcccgagc                                                                                                                                                  >  2:195104/1‑216 (MQ=255)
                                                                                gTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTCAAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGGCATGGTCCCGGcca                                                                                                                                           >  1:568921/1‑237 (MQ=255)
                                                                                  cGCCCCCATGGCAGACCAGGCTGCGAAAATTACGTTAGCCCCGGCAGTAAAATTACAGAAAGGCAATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGCCCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCaaa                                                                                                                                          <  1:643941/238‑1 (MQ=255)
                                                                                    cTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGGCCAGGGCCCAGCAAACa                                                                                                                                        >  1:1149783/1‑238 (MQ=255)
                                                                                    cTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑AACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGGCCCGGCaaaaa                                                                                                                                        >  1:930589/1‑236 (MQ=255)
                                                                                             cagcaCAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTCACC‑GGACCGCCGC‑CACCCCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGGCACGGTCCAGGAAAACAACAcccccca                                                                                                                              >  2:1046722/1‑234 (MQ=255)
                                                                                               gcaCAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑AGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCCCCGCGCCACCATAAGAAAGCCCCTCATGATCATCACGGCGGGCACGGTCCCGGCcaaaagacacactcaaa                                                                                                                            >  1:451967/1‑225 (MQ=255)
                                                                                                 acaGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCCCCATAAGAAAGCTCCTCATGATCATCACGGCGGGCACTGTCCACGCAAACATCACCGCCAAAAGa                                                                                                                          >  1:324268/1‑239 (MQ=255)
                                                                                                     gCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGGCAGGGGTCAGGCCAACACCACCCCCAAATCACaaa                                                                                                                      >  1:1106669/1‑239 (MQ=255)
                                                                                                                             tCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCCTTGTCCCGGCAAAAATCACCCCTAAATGACAAAATCCGGGGTACAATCCGGCATTc                                                                                               >  1:372738/1‑238 (MQ=255)
                                                                                                                                 tAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTGCCGGCAAACATCACCGCTAAATGACCAAAGCCGGGGAACAAACCCGCCTTCAggc                                                                                           >  1:1123087/1‑236 (MQ=255)
                                                                                                                                       tacggtagggCGATCGTTATAATCGTGTCTTTTACGGGGGTGGGGGCCACTGGCGGG‑ACCACGGTGGGGGGAACAAACATTATGAATGGCGG‑GGCAATCGCTGGCCCCTACAC‑GGTCCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑Tg                                                                                   <  2:271285/231‑1 (MQ=255)
                                                                                                                                                          tAATCGGGGCTTTTCCGGGGATGGAGGTCCCTGGCGCG‑ACCCCGGCTGGTGGAAACAACCTTATGAAGTGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTAt                                                                <  1:1248190/239‑1 (MQ=255)
                                                                                                                                                          atagcggtgtattgacggggATGGAGGTCCGGCGCGCC‑GCCAGGGGTGGCAGAAACAACATTAGGAAGGGCGA‑GGCAATCGCTGGCCCCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTAt                                                                <  2:912166/222‑1 (MQ=255)
                                                                                                                                                          aaagcGTGGCTTTTCCGGGGATGGAGGTCACGGGCGCG‑ACCGCGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCCCCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTAt                                                                <  1:554169/235‑1 (MQ=255)
                                                                                                                                                                ttgataaaaagagagggggggccACGGGCCCG‑CCCCCGGCTGGTGGAAAAAACATTAAGAATGGCGA‑GGCCATCGGTGGCACCTACAC‑GGACCGCCCC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGcc                                                          <  2:452316/217‑1 (MQ=255)
                                                                                                                                                                 ggCTTTTCCTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGcc                                                          <  1:459826/238‑1 (MQ=255)
                                                                                                                                                                 ggCTTTTCCGGGGATGGAGGTC‑CTGGGCCGCACCACGGCTGGTGGAAACAACATTATGAAGGGCGA‑GGCAATCGCTGGCCCCCACCC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGcc                                                          <  1:1231706/238‑1 (MQ=255)
                                                                                                                                                                 ggCTTTTACTGGGATGGAGGTCACTGGCGCC‑GCCCCGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACCC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGcc                                                          <  1:73833/238‑1 (MQ=255)
                                                                                                                                                                    tttccggggggggggtgcccggggCGCG‑‑CCACGGGGGGGGAAAACAAGTTTATGAATGGCGA‑GGCAATCGCTGGCCCCTACGC‑GGCCCGCCGC‑CCCCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGt                                                     <  2:1149783/217‑1 (MQ=255)
                                                                                                                                                                        acggggATGGGGTGCACTGGCGCG‑ACCACGGCTGGGGGAAACAACATTAGGAAGGGCAG‑GGAAATCGCTGGCACCTACAC‑GGCCCGCCGC‑CACCGCCGGGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGTTaa                                                  <  1:793615/236‑1 (MQ=255)
                                                                                                                                                                        aCTGGGTGGGAGGTCAGTGGCCAG‑ACCACGGCTGGGGAAAACAACATTATGAAGGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGTTaa                                                  <  1:1000318/239‑1 (MQ=255)
                                                                                                                                                                          tggggtggggGTCCCGGGCGCG‑CCCGCGGCGGGGGGAACAAACATTTTGGATGGCAG‑GGCAATCGCTGGCCCCT‑CACGGGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAAtt                                                <  2:1236930/239‑1 (MQ=255)
                                                                                                                                                                          tGGGATGGAGGCTACTGGCGCC‑ACCAGGGCTGGTGGAAACAACATTATGAAGGGCGA‑GGAAAGCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAAtt                                                <  1:716152/239‑1 (MQ=255)
                                                                                                                                                                          ggggttggggtgcccTGGGGCG‑ACCCCGGCTGGGGGAACCCACATTTAGAAGGGCGG‑GGCAATCGCTGGCCCCTACAC‑GGCCCGCCGC‑CCCCGCCGGGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAAtt                                                <  2:320473/225‑1 (MQ=255)
                                                                                                                                                                              aTGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGCAAg                                                                                                                                                                                >  2:1170346/1‑108 (MQ=255)
                                                                                                                                                                              aTGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAg                                                                                                                                                                                <  1:1170346/108‑1 (MQ=255)
                                                                                                                                                                                ggAGGTCCCTGGGGCG‑CCCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCa                                           <  1:1046722/238‑1 (MQ=255)
                                                                                                                                                                                 gAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGGAAACATCACCGCTAAATGACAAATGCCGGGGAACAAACCGGCATTCAGCGCCTGA‑TTCGACGCCGGGGCGGCTTTTCAGGCCCAAGGTAATTTTGCaaa                                         >  1:1164901/1‑238 (MQ=255)
                                                                                                                                                                                        gtGGGGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGG‑GGCAATCGCTGGCCCCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGa                                  <  1:378576/238‑1 (MQ=255)
                                                                                                                                                                                         ctggcacacccgcggtgggggaaaaCAACATTATGAAGGGCGA‑GGCAATCGCTGGCCACTCCAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGa                                  <  1:195104/217‑1 (MQ=255)
                                                                                                                                                                                           gcgcG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCCAGCAAACATCCCCGCTAAATTAAAAAAGCCCGGGAACAAACCGGCATTCAACGCCTGA‑TGCGACCCTGGCCGGCCTTATCAGGCCTACGCTAATTCTGcccggctcttttac                               >  1:533756/1‑226 (MQ=255)
                                                                                                                                                                                           gcgcG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTCAATCACAAATGCCGGGGAACAATCCGGCATTCAGCGCCTGA‑TGCGGCGCTGGCGCGGCTTATCAGGCCTACGGTAAtt                                                >  1:452316/1‑222 (MQ=255)
                                                                                                                                                                                                ggccACGGTTTGGGGAAACAACATTAGGAATGGCAA‑GGCAATCGCTGGCCCCTACCC‑GGCCCGCCGC‑CACCGCCGGGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCa                           <  1:954073/237‑1 (MQ=255)
                                                                                                                                                                                                gccctgggtgggtggAAAAAACATTAGAAATGGGGG‑GGCAATCGCGGGCCCCTCACC‑GGCCCGCCCC‑CCCCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCa                           <  2:372738/229‑1 (MQ=255)
                                                                                                                                                                                                 ccgccGGTTGGTCGAAACAAAATGATGAAGGGCGA‑GGCAATCGCTGGCACCTCCAC‑GGACCGCCGA‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCa                           <  1:1188737/234‑1 (MQ=255)
                                                                                                                                                                                                  cccgggCCGGTGGAAACAAAATTAGGAATGGCGA‑GGCAATCGGTGGCCCCCACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATg                         <  1:8789/235‑1 (MQ=255)
                                                                                                                                                                                                            tggAAACACCATTAAGAATGGCGA‑GGCAATCGCTGGCCCCTACCC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCCCGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAg                <  1:476628/238‑1 (MQ=255)
                                                                                                                                                                                                             ggAAACAACATTAAGAATGGCTA‑TGCAATCGCTGGCACCCACCC‑GGACCCCCGC‑CACCGCCGCGACACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg               <  1:962449/238‑1 (MQ=255)
                                                                                                                                                                                                                     caTTATGAATGGCGA‑GGCAACCGCTGGCACCTACCC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGcc                                                                                         <  1:97259/157‑1 (MQ=255)
                                                                                                                                                                                                                     caTTAAGAATGGCGA‑GGCAATCGCTGGCACCTCCAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCCCCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAAtt                                                <  1:190891/197‑1 (MQ=255)
                                                                                                                                                                                                                         aTGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGGTATGGGCCAGGGCAAAATCACCCCTCAAAAACAAAACCCGGGGAACAAACAAGCATTCAGCCACCCA‑TACGACACCGGCGCGAGCTTTCAGGCCTCCTCGAATTCTGCACTATATTTAATACGCATAACTTTGTTGGccagaaaaaaaga  >  2:252355/1‑227 (MQ=255)
                                                                                                                                                                                                                          tGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGGCCAGGCAAACATCCCCGCTCAAAAACAAATGACGGGGAACAATCCCGCATTCAGACCCGCACTGCG‑CGCTGGGGCGCCTTTTCAGGACCCCGTTGATTTTGCAATAAATTTAATCTGCATGCTTTTTTTGGCAGGAGAAggggg  >  2:925635/1‑235 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                         
ACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGA‑TGCAATCTATCGTACTCGCACTTTCCCTGGTTC‑TGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCG‑ACCACGGCTGGTGGAAACAACATTATGAATGGCGA‑GGCAATCGCTGGCACCTACAC‑GGACCGCCGC‑CACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGT  >  NC_000913/5201‑5653

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