Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,043,013 A→G 22.4% 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,0130AG22.4% 153.4 / 8.9 76G200G (GGA→GGGsrkACpx stress response Thr/Ser protein kinase, MazF antagonist protein
Reads supporting (aligned to +/- strand):  ref base A (12/44);  new base G (15/3);  total (28/48)
Fisher's exact test for biased strand distribution p-value = 4.08e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATGGAAGCGGTTGGGCGTTATTTAGGGCGTATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGC  >  NC_000913/4042775‑4043249
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
gggttgAAGCGGTTGGGCTTTTTTGGGGGCGTAGCACCCAGCAGGGGCGCAAAACGTTTTTTTCCCAGCGCCCTCCCATCGGTTTAAAGGAATATCCCATTGGGCCACGCAAGGTGTTTGGGGACGCTACACTGATCCCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACGGATGAGCTGTTTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGga                                                                                                                                                                                                                                                   <  2:849145/234‑1 (MQ=255)
    ggAAGCGGTTGGGCGTTATTTAGGGCGTATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGAGGGCAGGATTCTCCGTCCTGCGGGTCCATGgtgaa                                                                                                                                                                                                                                               >  1:622803/1‑237 (MQ=255)
        gCGGTTGGGCGTTATTTAGGGCGTATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAGTTCACCGTTTTGCGGGTACATGgagaagtgc                                                                                                                                                                                                                                           >  1:207702/1‑233 (MQ=255)
            ttGGGCGTTATTTAGGGTGTATCCACCAGACGGGGCCCAAAACGTTTTTTACCCATCCCCCGACCATCGGTTTGAAAGAAGATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTCCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGc                                                                                                                                                                                                                                      <  1:366209/238‑1 (MQ=255)
            ttGGGCGTTATTTAGGGCGTATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCCCACTGGCGGGAAAGATTTCCCGTTCTGCGGCCTCAAGGGGACTTcccagcc                                                                                                                                                                                                                                      >  1:227526/1‑234 (MQ=255)
            ttGGGCGTTATTTAGGGCGTATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTTTTCGGCTCCAAGGGGACTGCCAcccc                                                                                                                                                                                                                                     >  1:843320/1‑239 (MQ=255)
            ttGGGCGTTATTGAGGGCGTTTGCACCAGACGGGGCGCAAACAGCTTTTTCTCCATCGCCCGACCATCGGTTTGAACGAATATCTTATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGc                                                                                                                                                                                                                                      <  1:381800/238‑1 (MQ=255)
             ggggCGTTATTAGGGCGGTATGCAACAGTCGGGGCTCAAAAAGCTTTTTACCCATCGCCAGACAATGGGTTTGAACGAATTGCTCATTGAGCCAGGCAAGCTGTTTGAGGACGCTACACTGATCCCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACATATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCg                                                                                                                                                                                                                                    <  1:1045981/238‑1 (MQ=255)
                gCGTTATTTAGGGCGTATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGAGTTTACCGTTTTGCGGCTCCACAGGGAGTGCCACGCCcggg                                                                                                                                                                                                                                 >  1:256538/1‑238 (MQ=255)
                           ggCGTATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCCCACTGGCGGGCAGAATTCAACGTGCTGCGGCAACCAGGAGGAACCCACGCAGGGAATAAtcttcgg                                                                                                                                                                                                                      >  2:15743/1‑235 (MQ=255)
                                tATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACCCTGGCGGGAGAAATTTACCGTACTGGGGCGACCATAGAACAGCCACCCCGGGaaaaaatactgggggg                                                                                                                                                                                                                  >  2:382113/1‑225 (MQ=255)
                                tATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTAACCGTTCTCGGGGCACATGGGGAGTGTCACGCCAGGAGAATTTTCTGCcgcg                                                                                                                                                                                                                  >  1:220150/1‑239 (MQ=255)
                                 ttggcccacgacgggtgccAAACTGTTTTTTACCCCCCCCCCGACCGTGGGTTTGAACCGATATCGCTTGGCGCCACGCAACGTGTTTGAGGACGCTCCGCCGATTCCTGCCGGGTTGAAAGCGTCATTCCGAAAAGGGAAAGCTGGGCTGTGTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGcgcg                                                                                                                                                                                                                  <  2:227526/221‑1 (MQ=255)
                                 tgtcACCAGGCGGGGGGCCAACCGTTTTTATACCTGCGCCCGACCGTGGGTTGAACAGAAATCCCTATTGGGCAACGAAGGCTTTTGGAGGACGCTACACTGTTCCCTACCGGTTTGAAAGCGGAATTCCGAAAAGGGAAAGATGGGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGcgcg                                                                                                                                                                                                                  <  2:757850/235‑1 (MQ=255)
                                 aTGCACCAGATGGGTCGCAAACTGCTTTTTATCCAGCGCCCGACCTTCGGTTTGAACGAATATCTCATTGAGCCCCGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGa                                                                                                                                                                                                                 <  1:496735/239‑1 (MQ=255)
                                 aTGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTTTGCGGCCACACGGGGACTGCAACACCCGGGGAATTCTCTgggtcga                                                                                                                                                                                                                 >  1:70152/1‑234 (MQ=255)
                                 aTGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCAACCTTCTGCGGCTACATTGAGAATGCCCCCCCCGGGAGATTCTCTGGcgcgc                                                                                                                                                                                                                 >  1:1085401/1‑238 (MQ=255)
                                  tGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGGTCTGCGGCTACATGGGGACTGCCCCGCCGGGAATATTTTCTGGcgcgcg                                                                                                                                                                                                                >  1:1010910/1‑237 (MQ=255)
                                        caggcgggtccaaaaCAGCTTTTTTCCCATCGCCCGCCCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATCCCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGTGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTcc                                                                                                                                                                                                           <  2:405626/227‑1 (MQ=255)
                                        cAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAAGATTTAACGGTCTGCGGGTACCAGGGGACTGCCCCCCCCGGGATATTTTCTgggcgccggggcc                                                                                                                                                                                                           >  1:853655/1‑228 (MQ=255)
                                         cgACGGGGCCAAAATCGTTTTTTCTCCATCGCCCGCACATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATCCCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCa                                                                                                                                                                                                          <  2:256538/238‑1 (MQ=255)
                                              gggCCCAACCGGTTTTTTTCCCCTCCCCCGGCCTTGGGTTTGAACGAATTCCCCTTTGGCCCACGCAAGCTGTTTGGGGACCCTATACTGATCCCTTCCGGGTTGAAAGCGTCATTCCGGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAAt                                                                                                                                                                                                       <  2:207702/236‑1 (MQ=255)
                                               ggCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAGAAATTCACCCTTTTTCGGCCCACTGGGGACTGCCCCCCCGCGGATATtctcctgtcgggagggccaacgatt                                                                                                                                                                                                    >  2:605695/1‑218 (MQ=255)
                                                gcgcAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAAAATTTACCCTTCTGCGGCTACATGGGGACTGCCACGCCCGGGGAATTTTCCGGCGCGCTGGTTCAAtggttg                                                                                                                                                                                                  >  1:589702/1‑239 (MQ=255)
                                                   caaaacgCTTTTTACCCATCCCCCGGCCATCGGTTTGAACGAATATCCCATTGGGCCCCGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAAtgtttgtt                                                                                                                                                                                                <  1:1025971/232‑1 (MQ=255)
                                                         gCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAAATTTCACCGCTCTGCGGCTACATGGGGAATGCCACCCCCGGGCAATTTCCTGGGCCGACGGTCCACTGTtttttgtctcgt                                                                                                                                                                                         >  1:319523/1‑233 (MQ=255)
                                                                 ttcccccgccctgccttggttttGAATGATTTTCTCTTTGAGCCACGCAGGTGTGTGAGGGCCGCTACACTGTCACCTTCGGGTTTGAAGGCGGCATTCCGGAAAGGGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTAACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGc                                                                                                                                                                                  <  2:495096/219‑1 (MQ=255)
                                                                 acccAGCGCCCGACATGCGGTTTGACCGAATATCTCATTGAGCCACGAAAGCTGTTTGAGGACGCTACACGGATCCCTTCGGGTTTGAAAGCGGCATTCCTGAAAGGGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGc                                                                                                                                                                                  <  1:728063/237‑1 (MQ=255)
                                                                       cGCCCGACCTTCGTTTTGAACGAATATCTCTTGGAGCCACCAAGGTTTTTTGAGGCCGCTACCCGGATGCCTTCCGGGTGGAAAGCGGCATTCCTGAAGGCGACAGATGAGCTGATGGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTaa                                                                                                                                                                            <  2:319523/239‑1 (MQ=255)
                                                                          ccGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTGCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGa                                                                                                                                                                                                                                 <  1:463800/184‑1 (MQ=255)
                                                                               cctCGTTTTGAACGAATATCTCATTGGCCCACGGAAGTTTGTGGAGGACGCGACATTGATCCCTCCTGGGTTGAAGGCGGCATGCCGAAAAGCGACAGTGTAGCTGATTGCCGCCGTTACAGCACCCTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAg                                                                                                                                                                    <  2:220150/237‑1 (MQ=255)
                                                                                     tttGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTCCACTGATCCCTTCGGGGTTGAAGGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTc                                                                                                                                                              <  1:446374/239‑1 (MQ=255)
                                                                                      ttGAACGAAATCCCCATGGGGCCACGAAAGCTGTTGGAGGACGCTCCACTGATCCCTTCCGGGTTGAAAGGGGCATTCCTGAAAGCGACAGAGGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTc                                                                                                                                                              <  1:965985/238‑1 (MQ=255)
                                                                                        gAACGAATATCCAATGGGGCCACGCAAGCTGTTTGGGGACGCTACACTGATCCCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAgg                                                                                                                                                          <  1:382113/239‑1 (MQ=255)
                                                                                              ttttctcTTTGACCCCCGCAAGCTGTTGTAGGACCCTACCTTGATTCCTTCGGGTTTGAAGGCGGCTTTCCGAAAGGCGACAGTTGAGTTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTtgt                                                                                                                                                    <  2:853655/236‑1 (MQ=255)
                                                                                                    taaTTGAGCCCGGCAAGGTGTTTGAGGACCCGAGACTGATGCCTTCCGGGTTGAAAGCGGCATTCCGGAAGGCGACAGATGAGCTGATGGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATg                                                                                                                                               <  1:409899/237‑1 (MQ=255)
                                                                                                        tGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTACTGAAAGCGACAGATGAACTGATTGCCGCCGTTACAGCACCCGGGCGGGAAACATTTAACGCTATGCGGGTACAGAGAGAACTAGCACACGGGGAAAATACCTGCTGGCGTGGGCAAATGTTTGTTGGAATGGATGAAGCACGTAATGGGGCCGCCGGCCCGGAGTTGTgggggtggg                                                                                                                                           >  2:185596/1‑232 (MQ=255)
                                                                                                          aGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTACCGTTCTGCGGCTAAATTGGGACTGCCCCGCCGGGGATATTTTCTGGGGCGCTGGGCCCATGTTTGTTTATCTGGATGATGGACGGAATGGTCCAGCCGTTCAGGATTTTTGGGTGTTGTTc                                                                                                                                        >  1:167459/1‑239 (MQ=255)
                                                                                                                         ttGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTTTCTGGCGCGATTGTTCAATGTTTGTTTGTCTGGGTGATGCACGGAATGGTCCAGCCGTTCAGGGTTTTTTGATGTTGTTCAATGGGGTTAAAAcc                                                                                                                         >  1:829308/1‑239 (MQ=255)
                                                                                                                               gACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGGGACTGCCACGCCGGGGATATTCTCTGGCGCGAGGGGCCAATGTTTGTTGATCTGGATGATGCACGGAATGGTCCAGCCGTTCAGGATTTTTGGATGTTGTTCAATGGGGGTAAAAcccggca                                                                                                                    >  1:129724/1‑234 (MQ=255)
                                                                                                                               gACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTCTCTGGCGCGCTGGTCCCATGTTTGTTTATCTGGATGATGCACGTAATGGTCCAGCCGTGCCGGGTTTTTGGGTTTTGTTCAAAGGCGCTAAAGCCCagca                                                                                                                    >  1:217957/1‑239 (MQ=255)
                                                                                                                                aCGCTACACTGATCCCTGCCGGGTTGAATGCGTCATTCCTGAAGGCGACAGAGGTGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGagca                                                                                                                    <  1:948900/238‑1 (MQ=255)
                                                                                                                                 cGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAAAGTTTACCCTTCTTCGGCTCCCTCGAGACTGCCCCGCCCGGGAATATCTCTTTCCCGCTTGTCCACAGTTTGTTGGTTTTGGTGATGGTCGCAATGGTGCAGCCCTTCAGGATTTTTGGGTGTTTCTCAAAGGCGATAAAAACGagccgc                                                                                                                  >  2:211421/1‑239 (MQ=255)
                                                                                                                                   cTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTCTCCGTTCTGCGGGCCCATGGGGAATGCCCCCCCGGGGAAATTTTCTGTCGCGAGGGTCCCATGGTTGTTGTTCTTGGGGAGGCACGGAATGGTCCCGCCCGCCAGGGTTTTGTTGGTGTGTCTCAAGGCGGTAAAGcagcgaaaggg                                                                                                                >  2:306249/1‑229 (MQ=255)
                                                                                                                                    ttcgcggATGCCTTCCGGTTTAGATGCGGCATTCCGGAAGGCGAATGATGAGCTGGTTGCCGCCGTTACAGCCACCGGGCGGGAAGATTTCCCCGTTCGGGGGCTACATGGAGCCTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGa                                                                                                               <  2:217957/233‑1 (MQ=255)
                                                                                                                                    tACACTGATCCCTTCCGGGTTGAAAGCGGAATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTaa                                                                                                                                                                            <  1:1038332/180‑1 (MQ=255)
                                                                                                                                      cacgtaTACCTTCCGGGTTAAAAGCGGCATTCCTGAAAGCGACAGATGACCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATg                                                                                                             <  1:367201/234‑1 (MQ=255)
                                                                                                                                      cacTGCTCCCCTCCGGGTTAAAAGCGGCATCCCGTAAAGCAAGAGTGGGGCTGCTTGCCGGCGTTACAGCACCCTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATg                                                                                                             <  2:843320/239‑1 (MQ=255)
                                                                                                                                      cacTGATCCCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGa                                                                                                                                                                                                                 <  1:378560/142‑1 (MQ=255)
                                                                                                                                         tGATACCTTCCGGGTTGAAAGCGGCATGCCGGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCa                                                                                                                                                                                 <  1:1068091/170‑1 (MQ=255)
                                                                                                                                            tACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGCGACTGCCACGCCGGGGATATTCTCTGGCGCGATGGTGCAAtg                                                                                                                                                                                                      >  1:265452/1‑146 (MQ=255)
                                                                                                                                             cccTTGCGGTTGGAAAGGGGCATCCAGAAAGCGGAAAGATGAGCTTAATGCCGACGTTACAACACACTGGCGGGAAGATTTCTCCGTTCTGCGTCTACATGGCGACTGCCAAGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTgg                                                                                                      <  2:1091746/238‑1 (MQ=255)
                                                                                                                                              ctttCGGGGTTGAACGCGGATTTCCTGAAAGGGACAGTGGGGGTGTTGGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                                     <  2:129724/237‑1 (MQ=255)
                                                                                                                                              ctgtCCGGTTTAAAAGCGGAATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                                     <  1:829951/236‑1 (MQ=255)
                                                                                                                                              ccTTCGGGTTTGAAGGCGAGTTTCCTAAAGGCGACAGATGAGCTGATTGCCGCCGTTCCAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                                     <  1:686/239‑1 (MQ=255)
                                                                                                                                              ccTTCCGTGTTGAAAGCGGCATTCCTGAAACCGACAGATGAGCTGATTGCCGCCGTTACAGCACTCTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTGAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                                     <  1:93704/239‑1 (MQ=255)
                                                                                                                                              ccTTCCGGGTTGAAAGCGGAATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                                     <  1:306249/239‑1 (MQ=255)
                                                                                                                                               ctgcggtgTTGAACGAGGCATCCATTAAGGCGACAGATGTGTTGTTTCGCGCGGTTGCAGCCCCCTGGCGGGAAGATTTCCCGTTTGTGTGGCTAAATGGGGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                                     <  2:70152/231‑1 (MQ=255)
                                                                                                                                                 tCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTATGCGGCCGCAAGGAGACTTCCACGCCGGGAATAGTCTCTTGCGCGATGGTCCCACGTTTTTTGATCTTGGTGATGCTCGGACTGGTCCAGCCGTTCAGGGTTTGTGGGGGTTGCGCAATTGCGATAAAACCGAGCACCggggggagaggggaaa                                                                                                  >  2:607320/1‑225 (MQ=255)
                                                                                                                                                   cGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGGTTTCACCCGTCCGCTGCTCAATGGGGAATGCCACGCCGGGGAAATTTTCTGGCCGGATTGTCCAATGTTTGTTTATCTGGATGATGCACGTAATTGGCCCGCCGGTCTGGGTTTTTGGGTGTTTGCTCATGGCGATAAAGACGCGCAGGGGGAGCAACAGGGAAAta                                                                                                >  2:6305/1‑239 (MQ=255)
                                                                                                                                                          aaaGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCa                                                                                                                                                                                 <  1:154353/153‑1 (MQ=255)
                                                                                                                                                          aaaGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACCCCGGGAATATTTTCTGGCGCGATGGTCCAATGTTTGTTGATCTTGATGATGCACGTAATGGTCCAGCCGTTCAGGGTTTTTGGATGTTGCTCAACGGGGGTAAAACCGAGCAGCGGATGCAACTGGAAACTattattta                                                                                         >  1:2573/1‑237 (MQ=255)
                                                                                                                                                          aaaGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAAATTTCCCCCCTTCGCGGGTCCAAGGGAAATGCCACCCCGGGAATATTTCCTGGGCGGACGGGCCAACTGTTTTTTGCTCGGCGGGTGCACGAAATGGGCCAGCCGGTCGGGGATTGGGGGGTTTGTTCAAAGGGGCTAAAACCCAGCCGCGGAGGCAAATGGAAAAAATTATTGa                                                                                         >  2:831286/1‑239 (MQ=255)
                                                                                                                                                          aaaGCCGCTTTCTGGAAAGAGACAGGTGAGCTGTTGCCCGCGGTTACAGCACCCTGGCGGGAAGATTTCCCCGTTCGGCGGCTACATGGAGGCTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGa                                                                                         <  2:1085401/239‑1 (MQ=255)
                                                                                                                                                                   ttCCGGAAGGGGACAGGTGGGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCCCCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATga                                                                                <  1:211421/239‑1 (MQ=255)
                                                                                                                                                                                            aTTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGc                                                        <  1:15743/239‑1 (MQ=255)
                                                                                                                                                                                                      ttACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGa                                              <  1:185596/239‑1 (MQ=255)
                                                                                                                                                                                                                 tGGCGGGAAGATTTCACCGTTCTGCGTCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGc                                                                                      <  1:921675/188‑1 (MQ=255)
                                                                                                                                                                                                                            aTTTCCCCGTTCGGCGGCTCCATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGa                                              <  1:42284/218‑1 (MQ=255)
                                                                                                                                                                                                                                   cGTTTTGGGGCTACATGGGGGTGGCCACGCCGGGAATTTTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGtt                  <  2:829308/239‑1 (MQ=255)
                                                                                                                                                                                                                                   cGTTCTGCGGCTACATGTAGCCTGCCACGCCGGGAATATTCTGTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCCCGTAAGTGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGtt                  <  1:948163/239‑1 (MQ=255)
                                                                                                                                                                                                                                   cGTTCTGCGGCTACATGGGGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGtt                  <  1:466430/239‑1 (MQ=255)
                                                                                                                                                                                                                                   cGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGtt                  <  1:605695/239‑1 (MQ=255)
                                                                                                                                                                                                                                      tCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGACGAATTTAGCAAGTTCGACACACCTCAAATCGGACTGTTTTAAcccattttcttttaacaaacc                >  2:42284/1‑219 (MQ=255)
                                                                                                                                                                                                                                       cTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCCGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGt              >  1:563373/1‑239 (MQ=255)
                                                                                                                                                                                                                                          cGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGATTATGAAGAATTTAGCGAGTTCGAACCCCCTGAAATCGGAATGATTGAAACTTTACCCCCGCTGCGTTTGGtttt           >  2:1152/1‑238 (MQ=255)
                                                                                                                                                                                                                                                    ggAGACTGCCACGCCGGGAATTTTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGc  <  1:445251/239‑1 (MQ=255)
                                                                                                                                                                                                                                                    ggAGACTGCCACGCCGGGAATATTCTCGGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGc  <  1:85173/239‑1 (MQ=255)
                                                                                                                                                                                                                                                    ggAGACGGCCACGCCGGGAATATTCTCTGGCTCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGc  <  1:28870/239‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
AGATGGAAGCGGTTGGGCGTTATTTAGGGCGTATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGC  >  NC_000913/4042775‑4043249

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: