Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,043,013 A→G 15.6% G200G (GGA→GGG)  srkA → Cpx stress response Thr/Ser protein kinase, MazF antagonist protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,043,0130AG15.6% 133.5 / 2.8 62G200G (GGA→GGG) srkACpx stress response Thr/Ser protein kinase, MazF antagonist protein
Reads supporting (aligned to +/- strand):  ref base A (11/40);  new base G (8/3);  total (19/43)
Fisher's exact test for biased strand distribution p-value = 2.01e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGAAGCGGTTGGGCGTTATTTAGGGCGTATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTA  >  NC_000913/4042778‑4043240
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     
tGGAAGGGGTTGGGCTTTTTTTAGGGCCTAGGCCCCAGTCGGGGCGTAAACAGCTTTTTATCCATCGCCCGACCATCGTTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATCCCTTCC‑GGTTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGga                                                                                                                                                                                                                                       <  2:681179/236‑1 (MQ=255)
 ggAAGCGGTTGGGCGTTATTTAGGGCGTATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCCACCTGGCGGGGCGAATTGACCG‑ATCTGTGGCGACATGgagcc                                                                                                                                                                                                                                    >  2:533778/1‑236 (MQ=255)
         ttGGGCGTTATTTAGGGCGTAGTCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACAGTCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGGGGACGCTACACTGATCCCTTCC‑GGGTGGAAAGCGGAATTCTTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGc                                                                                                                                                                                                                           <  2:698977/238‑1 (MQ=255)
                        ggCGTATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTAAACG‑ATTTCGGGGTACATGGGGAAAGCCCACGCGGGGAAAA‑TT‑TCTgg                                                                                                                                                                                                          >  2:896039/1‑239 (MQ=255)
                          cGTATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTCCAGACCACTGGCGGGGAGAATTCACCGCTTTTGCGCCCACATGGGTACTC‑CCACCCCGGGAATA‑‑TCTCTGgcc                                                                                                                                                                                                        >  2:472750/1‑238 (MQ=255)
                             ttgggccAGACGGGGCACAAACAGCTTTTTATCCAGCGCCCGACCATCGGTTTGAACGAATATCTCATTGGGCCACGCAAGCTTTTTGAGGACGCTACACGGATCCCTTCC‑GGTTTGAAAGCGGCATTCCGGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGa                                                                                                                                                                                                     <  2:322661/234‑1 (MQ=255)
                             aTGCCCCAGACGGGGCGAAAACAGCTTTTTATCCATCGCCCGACCTTCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGa                                                                                                                                                                                                     <  2:885422/239‑1 (MQ=255)
                             aTGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGGATTTACCG‑TTCTCCGGGCTCACGGGGACTC‑ACCAGCCGGGGATA‑TTTTTTGGcgcgc                                                                                                                                                                                                     >  2:11976/1‑238 (MQ=255)
                             aTGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCC‑TTTTGCGGCTCCATGGAGAATG‑CCACACCGGGGATA‑TTCTCTGGcgcgc                                                                                                                                                                                                     >  2:932623/1‑238 (MQ=255)
                             aTGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCC‑TTTTGCGGCTACAATGAGGCTG‑GCACGCCggggaaaatttcttggccca                                                                                                                                                                                                      >  2:319140/1‑222 (MQ=255)
                             aTGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGAAAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTACCC‑TTCTTCGGCCACGGTGAGACTG‑GCCCGCCCGGGAAA‑TTttctttgcgcg                                                                                                                                                                                                     >  2:222477/1‑232 (MQ=255)
                              tGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCC‑GTCTTCGGCTCCCAGAGGACTG‑CCACGCCGGGGAAA‑TTTCCTTGcgcggt                                                                                                                                                                                                    >  2:423764/1‑237 (MQ=255)
                                   cAGGCGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGGCCCACGCAAGCTGTTTGAGGACGCTACACTGATCCCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTcc                                                                                                                                                                                               <  2:204958/239‑1 (MQ=255)
                                    ggagggggCGCAAACAGGTTTTTAACCAGCGCCCGACCACCGGTTTGAACGAAATTCTCATTGAGCCACGCAAGCGGGTTGGGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGAGTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCa                                                                                                                                                                                              <  2:521089/235‑1 (MQ=255)
                                      aCGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTCCCG‑TTCTGCGGCCACCTGGGGACTG‑CCACCCCGGGGATA‑TTTTCTCGCGCGCTGGTGCAAt                                                                                                                                                                                            >  2:203049/1‑239 (MQ=255)
                                      aCGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCG‑TTCTGCGGGTCCCTGGAGAATG‑CCACCCCGGGAATATTTCTTTGGC‑CGATGGTccc                                                                                                                                                                                              >  2:347430/1‑236 (MQ=255)
                                      aCGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGAGGAATG‑CCACGCCGGGAAAA‑TTTTCTGTGGCGATGGGCCaaa                                                                                                                                                                                            >  2:937689/1‑238 (MQ=255)
                                       cGGGGCGCAAACAGTTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCa                                                                                                                                                                                              <  2:92236/236‑1 (MQ=255)
                                          ggCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGGTACATGGGGACTG‑CCCCGCCGGGGGAA‑TTTTCTgggggcggggtccaaacttt                                                                                                                                                                                        >  2:937422/1‑221 (MQ=255)
                                            cgaaaaCCGCTTTTTATCCATCCCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGAAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAAtgtttg                                                                                                                                                                                       <  2:1075338/235‑1 (MQ=255)
                                                           acccACCGCCCGATCATGGGTTTGAACGAATTTCTCATTGGCCCACGCAAGCTGTTTGAGGACGCTCACCTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGc                                                                                                                                                                       <  2:329414/237‑1 (MQ=255)
                                                            gccATCGCCCGGACATGGGTTGGAACGAATATCTCATTGGGCCACGCAAGCTGTTTGAGGCCGCTACACTGATACCTTCC‑GGGTTGAAAGCGGAATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGGCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCa                                                                                                                                                                      <  2:391890/238‑1 (MQ=255)
                                                             ccATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGGTTTCACAC‑TTTTGCGGGCCCATTGGGACTG‑CACCCCCGGGGAGA‑ATTTCTTGCGCGCTGGGCCAAAGGATGTGTATATGgatggggggc                                                                                                                                                                     >  1:617568/1‑233 (MQ=255)
                                                                 cctcgtccatgtgtttttaccGAATTTCTTATGGGGCCACCCATGCTTTTTGAGACCGCTCAACTGGCCCCTTCC‑GGGTGGAAGGGGGCATTCCGAAAAGCAACAGAGGCGGTTA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTaa                                                                                                                                                                 <  1:98769/219‑1 (MQ=255)
                                                                        ccATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACCCTACCCTGATACCTTCC‑GGGTTGAAAGCGGCATTCCGAAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCa                                                                                                                                                          <  2:904240/239‑1 (MQ=255)
                                                                              gTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑c                                                                                                                                                                                                                                <  2:1115563/164‑1 (MQ=255)
                                                                                         tatCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCCCCCCGGGGATA‑TTCTCTGGCGCGATGGTCCAATGTTTTTTTATCTGGGTGATGCACGTAATGGTCCCGCCGCTCCGGATTtttg                                                                                                                                         >  2:10385/1‑239 (MQ=255)
                                                                                            ctcATTGAGCCACGCAAGGTGTTGGAGGTCGCTTCACTGATCCCTGCC‑GGTTTGAAAGCGGCATTCCTGAAAGCGACAGATGATCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGC‑CGTAATGGTCCAGCCGTTCAGGATTTGTGGATg                                                                                                                                     <  2:295568/239‑1 (MQ=255)
                                                                                             ctattgaccccgccacgCTTTTTGAGGACCTTACACGTATCCCTCCC‑GGGTTGAAGGCGGCTTTCCTGAAAGCGACAGATGAGGTGT‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCCCCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATg                                                                                                                                     <  1:937689/223‑1 (MQ=255)
                                                                                                      ccagcacgttGTTTGAGGACGCAACGAGTATTCTGTCG‑GGGTGAAAAGCGGCATCCCTGAAAGCGAAGGGTGAGGTGA‑TTGCCGCCGTTACAGCCCACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAAt                                                                                                                            <  1:423764/230‑1 (MQ=255)
                                                                                                                    gaggcaccTACGCTGATCCCTCCG‑GGGTTGAAAGGGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGa                                                                                                              <  2:213640/232‑1 (MQ=255)
                                                                                                                        aCGCTGCACTGATTCCTTCC‑GGGTTGAAAGCGGCATCCGTGAAAGCGACAGATGGGCTGA‑TTGCCGCCGTTACAGCACCCTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGagca                                                                                                           <  2:374131/238‑1 (MQ=255)
                                                                                                                          gttaaccttatcccTTCG‑GGGTTGAAAGCGGCATTCCTGAAAGCAAAAGATGAGCTGT‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGagcagc                                                                                                         <  2:862001/226‑1 (MQ=255)
                                                                                                                              cacTGATCCCTTCG‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAAt                                                                                                                                                                                            <  2:1183006/151‑1 (MQ=255)
                                                                                                                                    ttccTTCCGGGGTTTAACGCGGCTTTCTAGAAAGCGAAGGTGGAGTTGA‑TCGCCGCCGTTACAGCCCACGGGCGGGATGATTTCACCG‑TTCTGCGGCTACATGGAGATTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACt                                                                                               <  1:1201197/237‑1 (MQ=255)
                                                                                                                                      ccTTCG‑GGGTTGAAAGCGGTATTCCGGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                            <  2:380811/239‑1 (MQ=255)
                                                                                                                                      ccTTCG‑GGGGTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑GTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCCGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                            <  2:75887/239‑1 (MQ=255)
                                                                                                                                      ccTTCC‑GGGTTGAAGGCGGCATTCCTGAAAGCGACAATTGGGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                            <  2:617568/239‑1 (MQ=255)
                                                                                                                                      ccTTCC‑GGGTTGAAAGCGGCATTCCTGAAAACGGAGGTGGAGGTAA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTT‑CCACTCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                            <  2:922686/239‑1 (MQ=255)
                                                                                                                                      ccTTCC‑GGGTTAAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                            <  2:1138842/239‑1 (MQ=255)
                                                                                                                                       cctgcggtgTTGAAAGCGTCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACGCTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                            <  2:415429/231‑1 (MQ=255)
                                                                                                                                              ggTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTTCGGCTACATTGAGACTG‑CCACGCCGGGGATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGGCCAGCCGGTCAGGATTTGTGGGTGTTGCTCAATGGCGGTAAAACCCGGCAGCGGGTGCAAcctgtaactttt                                                                                     >  2:491441/1‑228 (MQ=255)
                                                                                                                                                  acaaGCGGACTTGCTGAAAGCGACAGAGGGGGTGA‑TTGCCGCCTTTACAGCCCACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGGCTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTg                                                                                 <  2:916980/237‑1 (MQ=255)
                                                                                                                                                   aaaGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGGATA‑TTCTCTGGCGCGATGGTCCAATGTTTTTTGATCTGGGTGATGCACGGAATGGTCCAGCCGTTCAGGATTTTTGGATGTTGCTCAATGGGGATAAAACCGGGCAGCGGGTGCacttgttcttttacccat                                                                                 >  2:262021/1‑221 (MQ=255)
                                                                                                                                                    aggCGGCATTCCTGAACGCGACAGA‑CTGGTGATTTGCGTGAGTTACAGCACACTGGCGGGAAGAGTTCCCCG‑TTCTGGGGGTACATGGAGTCTG‑CCACGCCGGGAATT‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGTTGATGCACGTAATGGTACAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGa                                                                                <  2:456430/236‑1 (MQ=255)
                                                                                                                                                    aaGGGGCATTCCTGAAAGGGAAAGATGAGCTGA‑TTGGCCACGTTACAGCACACTGGCAGGAAGATTTCCCCG‑GTCTGCGGCTACATGGATACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCGGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGa                                                                                <  2:781535/238‑1 (MQ=255)
                                                                                                                                                          cATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGGCCAATGTTTGTTGGTCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGGTTTGTGGGTGTTGGTCAATTGCGGTAAAGCCGAGCAGCGGGTGCAACTGGGAACTATTATTTAAGGTTAt                                                                         >  2:843735/1‑239 (MQ=255)
                                                                                                                                                          cATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTACCG‑TTCTGCGGGCACATGGAGAATG‑CCACGCCGGGGATA‑TTTTCTGGCGCGGGGGTCCCATGGTTGTTGATCTGGATGTGGCACGGAATGGTCCAGCCGTTTCGGATTTTTGGTGTTTTGTCAAATGCGATAAGACGGAGCGGCGGGTGCAACGGGgaaaaaataatgtagcgaat                                                                         >  1:1235730/1‑222 (MQ=255)
                                                                                                                                                                       gACAGATGAGCTAA‑TTGCCGGCGGTACAGCACCCTGGCGGGAAGATTTCACCG‑TTCTTGGGCTACAGGGATACTG‑ACACGCCGGGAATA‑TTCTTTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGc                                                             <  2:1292083/238‑1 (MQ=255)
                                                                                                                                                                         acatgcgtcgttgttGCCGCCTTTACCGCACACGGGCGGGAAGTGTTCACCG‑TTCTGGGGCTACAGGGGGACTG‑CCACGTCGGGATTA‑TTCTCTGGCGTGTTTGGCCAATGTTTGTTGATCTGTATGATGCACGTAATGGTCCAGCCGTTCAGTATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGc                                                             <  1:533778/224‑1 (MQ=255)
                                                                                                                                                                               ggCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGa                                                     <  2:1235730/237‑1 (MQ=255)
                                                                                                                                                                                    a‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTTTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTTGCGCGATGGGCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGGTTTTTGGGTGTTGCTCCATGGCGATAAAACCGGGCAGCGGGGGCAACTTGGAACATTTATTTAAACTTTTTAAGAAttttagaggttgtgaaccc                                                >  2:1198153/1‑222 (MQ=255)
                                                                                                                                                                                     tttGCCGCCGTTACAGCACAGTGGCGGGAAGATTGTCCGG‑TTCTCGGGCTGCAGGGAGACTG‑CCACGTCGGGAATA‑GTCTCTGGCGGGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGc                                               <  1:347430/238‑1 (MQ=255)
                                                                                                                                                                                     gtttccggtgTTAAAGAACCCTGGGGTAAAGAGTTCGACG‑TTGTGCGTGTCCACGGAGACCT‑TCCCGCAGGGAATA‑CTTTATGGCGCTATGTTCCAATGTTTGTTGAGCGGGATGATGCACGTAATGGCCCATCCGGTCAGGATTTTTGGATGTTGCTTAATGGCGATAAAGACGAGCAGTGGATGCAACTGGAAACCATTCTTGAAGTTTATGAAGAATTTAGCGAGTTCGACACCGc                                               <  1:1071178/230‑1 (MQ=255)
                                                                                                                                                                                             cgTTAGAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGga                                                                                                                                                                             <  2:1074621/105‑1 (MQ=255)
                                                                                                                                                                                               gtaccacacaCGTGCGGGGAGATTGCTCCG‑TGCTGCGGCTATGTGGGGATCG‑CCCAGACGTGAATT‑GTCTCTGTCGCTAGGGGCAAATGTTTGGTGGTCTGGATGGTGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACGGGAAACTATTATTGAAGCATATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGa                                     <  1:1212630/233‑1 (MQ=255)
                                                                                                                                                                                                tacatgccacTGGGGGGAAATTTTCCCCT‑TTTGGCGGCTACATGGGGTCTG‑CCCCGCCGGGAATA‑TTCTCGGGCGCGATGGTCCAATTGTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGa                                     <  1:11488/231‑1 (MQ=255)
                                                                                                                                                                                                                    aTTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGGGcc                >  2:423918/1‑239 (MQ=255)
                                                                                                                                                                                                                     tttCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCa               <  2:493587/239‑1 (MQ=255)
                                                                                                                                                                                                                           cG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGtt         <  2:166554/239‑1 (MQ=255)
                                                                                                                                                                                                                           cG‑TTCGGCGTCTACATGGAGCCTG‑CCCCGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGtt         <  2:1272441/239‑1 (MQ=255)
                                                                                                                                                                                                                           cG‑TTCGGCGGCTACATGGATACTG‑CCACGCCGGGAATA‑TTCCCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGtt         <  2:454294/239‑1 (MQ=255)
                                                                                                                                                                                                                            g‑TTCTGCGGCTACATGGAGACTG‑CCACGGCGGGAATA‑TTCTCTGGCCCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGttt        <  2:964339/239‑1 (MQ=255)
                                                                                                                                                                                                                                    ggCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTtta  <  2:828132/238‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     
TGGAAGCGGTTGGGCGTTATTTAGGGCGTATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCC‑GGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGA‑TTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCG‑TTCTGCGGCTACATGGAGACTG‑CCACGCCGGGAATA‑TTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTA  >  NC_000913/4042778‑4043240

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