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→GGGsrkACpx 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

TGGAAGCGGTTGGGCGTTATTTAGGGCGTATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTA  >  NC_000913/4042778‑4043240
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     
tGGAAGGGGTTGGGCTTTTTTTAGGGCCTAGGCCCCAGTCGGGGCGTAAACAGCTTTTTATCCATCGCCCGACCATCGTTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATCCCTTCCGGTTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGga                                                                                                                                                                                                                                       <  2:681179/236‑1 (MQ=255)
 ggAAGCGGTTGGGCGTTATTTAGGGCGTATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCCACCTGGCGGGGCGAATTGACCGATCTGTGGCGACATGgagcc                                                                                                                                                                                                                                    >  2:533778/1‑236 (MQ=255)
         ttGGGCGTTATTTAGGGCGTAGTCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACAGTCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGGGGACGCTACACTGATCCCTTCCGGGTGGAAAGCGGAATTCTTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGc                                                                                                                                                                                                                           <  2:698977/238‑1 (MQ=255)
                        ggCGTATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTAAACGATTTCGGGGTACATGGGGAAAGCCCACGCGGGGAAAATTTCTgg                                                                                                                                                                                                          >  2:896039/1‑239 (MQ=255)
                          cGTATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTCCAGACCACTGGCGGGGAGAATTCACCGCTTTTGCGCCCACATGGGTACTCCCACCCCGGGAATATCTCTGgcc                                                                                                                                                                                                        >  2:472750/1‑238 (MQ=255)
                             ttgggccAGACGGGGCACAAACAGCTTTTTATCCAGCGCCCGACCATCGGTTTGAACGAATATCTCATTGGGCCACGCAAGCTTTTTGAGGACGCTACACGGATCCCTTCCGGTTTGAAAGCGGCATTCCGGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGa                                                                                                                                                                                                     <  2:322661/234‑1 (MQ=255)
                             aTGCCCCAGACGGGGCGAAAACAGCTTTTTATCCATCGCCCGACCTTCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGa                                                                                                                                                                                                     <  2:885422/239‑1 (MQ=255)
                             aTGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGGATTTACCGTTCTCCGGGCTCACGGGGACTCACCAGCCGGGGATATTTTTTGGcgcgc                                                                                                                                                                                                     >  2:11976/1‑238 (MQ=255)
                             aTGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCCTTTTGCGGCTCCATGGAGAATGCCACACCGGGGATATTCTCTGGcgcgc                                                                                                                                                                                                     >  2:932623/1‑238 (MQ=255)
                             aTGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCCTTTTGCGGCTACAATGAGGCTGGCACGCCggggaaaatttcttggccca                                                                                                                                                                                                      >  2:319140/1‑222 (MQ=255)
                             aTGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGAAAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTACCCTTCTTCGGCCACGGTGAGACTGGCCCGCCCGGGAAATTttctttgcgcg                                                                                                                                                                                                     >  2:222477/1‑232 (MQ=255)
                              tGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCCGTCTTCGGCTCCCAGAGGACTGCCACGCCGGGGAAATTTCCTTGcgcggt                                                                                                                                                                                                    >  2:423764/1‑237 (MQ=255)
                                   cAGGCGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGGCCCACGCAAGCTGTTTGAGGACGCTACACTGATCCCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTcc                                                                                                                                                                                               <  2:204958/239‑1 (MQ=255)
                                    ggagggggCGCAAACAGGTTTTTAACCAGCGCCCGACCACCGGTTTGAACGAAATTCTCATTGAGCCACGCAAGCGGGTTGGGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAGTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCa                                                                                                                                                                                              <  2:521089/235‑1 (MQ=255)
                                      aCGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTCCCGTTCTGCGGCCACCTGGGGACTGCCACCCCGGGGATATTTTCTCGCGCGCTGGTGCAAt                                                                                                                                                                                            >  2:203049/1‑239 (MQ=255)
                                      aCGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGGTCCCTGGAGAATGCCACCCCGGGAATATTTCTTTGGCCGATGGTccc                                                                                                                                                                                              >  2:347430/1‑236 (MQ=255)
                                      aCGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGAGGAATGCCACGCCGGGAAAATTTTCTGTGGCGATGGGCCaaa                                                                                                                                                                                            >  2:937689/1‑238 (MQ=255)
                                       cGGGGCGCAAACAGTTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCa                                                                                                                                                                                              <  2:92236/236‑1 (MQ=255)
                                          ggCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGGTACATGGGGACTGCCCCGCCGGGGGAATTTTCTgggggcggggtccaaacttt                                                                                                                                                                                        >  2:937422/1‑221 (MQ=255)
                                            cgaaaaCCGCTTTTTATCCATCCCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGAAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAAtgtttg                                                                                                                                                                                       <  2:1075338/235‑1 (MQ=255)
                                                           acccACCGCCCGATCATGGGTTTGAACGAATTTCTCATTGGCCCACGCAAGCTGTTTGAGGACGCTCACCTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGc                                                                                                                                                                       <  2:329414/237‑1 (MQ=255)
                                                            gccATCGCCCGGACATGGGTTGGAACGAATATCTCATTGGGCCACGCAAGCTGTTTGAGGCCGCTACACTGATACCTTCCGGGTTGAAAGCGGAATTCCTGAAAGCGACAGATGAGCTGATTGCCGGCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCa                                                                                                                                                                      <  2:391890/238‑1 (MQ=255)
                                                             ccATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGGTTTCACACTTTTGCGGGCCCATTGGGACTGCACCCCCGGGGAGAATTTCTTGCGCGCTGGGCCAAAGGATGTGTATATGgatggggggc                                                                                                                                                                     >  1:617568/1‑233 (MQ=255)
                                                                 cctcgtccatgtgtttttaccGAATTTCTTATGGGGCCACCCATGCTTTTTGAGACCGCTCAACTGGCCCCTTCCGGGTGGAAGGGGGCATTCCGAAAAGCAACAGAGGCGGTTATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTaa                                                                                                                                                                 <  1:98769/219‑1 (MQ=255)
                                                                        ccATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACCCTACCCTGATACCTTCCGGGTTGAAAGCGGCATTCCGAAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCa                                                                                                                                                          <  2:904240/239‑1 (MQ=255)
                                                                              gTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGc                                                                                                                                                                                                                                <  2:1115563/164‑1 (MQ=255)
                                                                                         tatCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCCCCCCGGGGATATTCTCTGGCGCGATGGTCCAATGTTTTTTTATCTGGGTGATGCACGTAATGGTCCCGCCGCTCCGGATTtttg                                                                                                                                         >  2:10385/1‑239 (MQ=255)
                                                                                            ctcATTGAGCCACGCAAGGTGTTGGAGGTCGCTTCACTGATCCCTGCCGGTTTGAAAGCGGCATTCCTGAAAGCGACAGATGATCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCCGTAATGGTCCAGCCGTTCAGGATTTGTGGATg                                                                                                                                     <  2:295568/239‑1 (MQ=255)
                                                                                             ctattgaccccgccacgCTTTTTGAGGACCTTACACGTATCCCTCCCGGGTTGAAGGCGGCTTTCCTGAAAGCGACAGATGAGGTGTTTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCCCCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATg                                                                                                                                     <  1:937689/223‑1 (MQ=255)
                                                                                                      ccagcacgttGTTTGAGGACGCAACGAGTATTCTGTCGGGGTGAAAAGCGGCATCCCTGAAAGCGAAGGGTGAGGTGATTGCCGCCGTTACAGCCCACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAAt                                                                                                                            <  1:423764/230‑1 (MQ=255)
                                                                                                                    gaggcaccTACGCTGATCCCTCCGGGGTTGAAAGGGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGa                                                                                                              <  2:213640/232‑1 (MQ=255)
                                                                                                                        aCGCTGCACTGATTCCTTCCGGGTTGAAAGCGGCATCCGTGAAAGCGACAGATGGGCTGATTGCCGCCGTTACAGCACCCTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGagca                                                                                                           <  2:374131/238‑1 (MQ=255)
                                                                                                                          gttaaccttatcccTTCGGGGTTGAAAGCGGCATTCCTGAAAGCAAAAGATGAGCTGTTTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGagcagc                                                                                                         <  2:862001/226‑1 (MQ=255)
                                                                                                                              cacTGATCCCTTCGGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAAt                                                                                                                                                                                            <  2:1183006/151‑1 (MQ=255)
                                                                                                                                    ttccTTCCGGGGTTTAACGCGGCTTTCTAGAAAGCGAAGGTGGAGTTGATCGCCGCCGTTACAGCCCACGGGCGGGATGATTTCACCGTTCTGCGGCTACATGGAGATTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACt                                                                                               <  1:1201197/237‑1 (MQ=255)
                                                                                                                                      ccTTCGGGGTTGAAAGCGGTATTCCGGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                            <  2:380811/239‑1 (MQ=255)
                                                                                                                                      ccTTCGGGGGTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGAGTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCCGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                            <  2:75887/239‑1 (MQ=255)
                                                                                                                                      ccTTCCGGGTTGAAGGCGGCATTCCTGAAAGCGACAATTGGGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                            <  2:617568/239‑1 (MQ=255)
                                                                                                                                      ccTTCCGGGTTGAAAGCGGCATTCCTGAAAACGGAGGTGGAGGTAATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTTCCACTCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                            <  2:922686/239‑1 (MQ=255)
                                                                                                                                      ccTTCCGGGTTAAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                            <  2:1138842/239‑1 (MQ=255)
                                                                                                                                       cctgcggtgTTGAAAGCGTCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACGCTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                            <  2:415429/231‑1 (MQ=255)
                                                                                                                                              ggTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTTCGGCTACATTGAGACTGCCACGCCGGGGATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGGCCAGCCGGTCAGGATTTGTGGGTGTTGCTCAATGGCGGTAAAACCCGGCAGCGGGTGCAAcctgtaactttt                                                                                     >  2:491441/1‑228 (MQ=255)
                                                                                                                                                  acaaGCGGACTTGCTGAAAGCGACAGAGGGGGTGATTGCCGCCTTTACAGCCCACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGGCTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTg                                                                                 <  2:916980/237‑1 (MQ=255)
                                                                                                                                                   aaaGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTCTCTGGCGCGATGGTCCAATGTTTTTTGATCTGGGTGATGCACGGAATGGTCCAGCCGTTCAGGATTTTTGGATGTTGCTCAATGGGGATAAAACCGGGCAGCGGGTGCacttgttcttttacccat                                                                                 >  2:262021/1‑221 (MQ=255)
                                                                                                                                                    aggCGGCATTCCTGAACGCGACAGACTGGTGATTTGCGTGAGTTACAGCACACTGGCGGGAAGAGTTCCCCGTTCTGGGGGTACATGGAGTCTGCCACGCCGGGAATTTTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGTTGATGCACGTAATGGTACAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGa                                                                                <  2:456430/236‑1 (MQ=255)
                                                                                                                                                    aaGGGGCATTCCTGAAAGGGAAAGATGAGCTGATTGGCCACGTTACAGCACACTGGCAGGAAGATTTCCCCGGTCTGCGGCTACATGGATACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCGGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGa                                                                                <  2:781535/238‑1 (MQ=255)
                                                                                                                                                          cATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGGCCAATGTTTGTTGGTCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGGTTTGTGGGTGTTGGTCAATTGCGGTAAAGCCGAGCAGCGGGTGCAACTGGGAACTATTATTTAAGGTTAt                                                                         >  2:843735/1‑239 (MQ=255)
                                                                                                                                                          cATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTACCGTTCTGCGGGCACATGGAGAATGCCACGCCGGGGATATTTTCTGGCGCGGGGGTCCCATGGTTGTTGATCTGGATGTGGCACGGAATGGTCCAGCCGTTTCGGATTTTTGGTGTTTTGTCAAATGCGATAAGACGGAGCGGCGGGTGCAACGGGgaaaaaataatgtagcgaat                                                                         >  1:1235730/1‑222 (MQ=255)
                                                                                                                                                                       gACAGATGAGCTAATTGCCGGCGGTACAGCACCCTGGCGGGAAGATTTCACCGTTCTTGGGCTACAGGGATACTGACACGCCGGGAATATTCTTTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGc                                                             <  2:1292083/238‑1 (MQ=255)
                                                                                                                                                                         acatgcgtcgttgttGCCGCCTTTACCGCACACGGGCGGGAAGTGTTCACCGTTCTGGGGCTACAGGGGGACTGCCACGTCGGGATTATTCTCTGGCGTGTTTGGCCAATGTTTGTTGATCTGTATGATGCACGTAATGGTCCAGCCGTTCAGTATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGc                                                             <  1:533778/224‑1 (MQ=255)
                                                                                                                                                                               ggCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGa                                                     <  2:1235730/237‑1 (MQ=255)
                                                                                                                                                                                    aTTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTTTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTTGCGCGATGGGCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGGTTTTTGGGTGTTGCTCCATGGCGATAAAACCGGGCAGCGGGGGCAACTTGGAACATTTATTTAAACTTTTTAAGAAttttagaggttgtgaaccc                                                >  2:1198153/1‑222 (MQ=255)
                                                                                                                                                                                     tttGCCGCCGTTACAGCACAGTGGCGGGAAGATTGTCCGGTTCTCGGGCTGCAGGGAGACTGCCACGTCGGGAATAGTCTCTGGCGGGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGc                                               <  1:347430/238‑1 (MQ=255)
                                                                                                                                                                                     gtttccggtgTTAAAGAACCCTGGGGTAAAGAGTTCGACGTTGTGCGTGTCCACGGAGACCTTCCCGCAGGGAATACTTTATGGCGCTATGTTCCAATGTTTGTTGAGCGGGATGATGCACGTAATGGCCCATCCGGTCAGGATTTTTGGATGTTGCTTAATGGCGATAAAGACGAGCAGTGGATGCAACTGGAAACCATTCTTGAAGTTTATGAAGAATTTAGCGAGTTCGACACCGc                                               <  1:1071178/230‑1 (MQ=255)
                                                                                                                                                                                             cgTTAGAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGga                                                                                                                                                                             <  2:1074621/105‑1 (MQ=255)
                                                                                                                                                                                               gtaccacacaCGTGCGGGGAGATTGCTCCGTGCTGCGGCTATGTGGGGATCGCCCAGACGTGAATTGTCTCTGTCGCTAGGGGCAAATGTTTGGTGGTCTGGATGGTGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACGGGAAACTATTATTGAAGCATATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGa                                     <  1:1212630/233‑1 (MQ=255)
                                                                                                                                                                                                tacatgccacTGGGGGGAAATTTTCCCCTTTTGGCGGCTACATGGGGTCTGCCCCGCCGGGAATATTCTCGGGCGCGATGGTCCAATTGTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGa                                     <  1:11488/231‑1 (MQ=255)
                                                                                                                                                                                                                    aTTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGGGcc                >  2:423918/1‑239 (MQ=255)
                                                                                                                                                                                                                     tttCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCa               <  2:493587/239‑1 (MQ=255)
                                                                                                                                                                                                                           cGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGtt         <  2:166554/239‑1 (MQ=255)
                                                                                                                                                                                                                           cGTTCGGCGTCTACATGGAGCCTGCCCCGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGtt         <  2:1272441/239‑1 (MQ=255)
                                                                                                                                                                                                                           cGTTCGGCGGCTACATGGATACTGCCACGCCGGGAATATTCCCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGtt         <  2:454294/239‑1 (MQ=255)
                                                                                                                                                                                                                            gTTCTGCGGCTACATGGAGACTGCCACGGCGGGAATATTCTCTGGCCCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGttt        <  2:964339/239‑1 (MQ=255)
                                                                                                                                                                                                                                    ggCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTtta  <  2:828132/238‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     
TGGAAGCGGTTGGGCGTTATTTAGGGCGTATGCACCAGACGGGGCGCAAACAGCTTTTTATCCATCGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTA  >  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: