Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,048,134 A→C N186K (AAT→AAG gfcC ← putative O‑antigen capsule production periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,048,1340AC98.6% 247.3 / ‑2.9 70coding (558/747 nt)gfcCputative O‑antigen capsule production periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (32/37);  minor base . (0/1);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-01

GGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGA  >  NC_000913/1047904‑1048363
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                     
ggCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTtcttct                                                                                                                                                                                                                                   >  1:12587/1‑235 (MQ=255)
ggCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTtcttct                                                                                                                                                                                                                                   >  1:285554/1‑235 (MQ=255)
          aGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCt                                                                                                                                                                                                                        >  1:627286/1‑236 (MQ=255)
            gCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGc                                                                                                                                                                                                                      >  1:48908/1‑236 (MQ=255)
                cTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTg                                                                                                                                                                                                                       <  2:627286/231‑1 (MQ=255)
                      aGATAAGAATTCTTCTTCATCAGTCAGGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCggggg                                                                                                                                                                                                            <  1:501831/236‑1 (MQ=255)
                       gATAAGAATTCTTCTTCATCAGTCAGGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGt                                                                                                                                                                                                           <  1:600990/236‑1 (MQ=255)
                          aaGAATTCTTCTTCATCAGTCAGGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGAt                                                                                                                                                                                                        >  1:426862/1‑236 (MQ=255)
                          aaGAATTCTTCTTCATCAGTCAGGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGAt                                                                                                                                                                                                        >  2:720791/1‑235 (MQ=255)
                              aTTCTTCTTCATCAGTCAGGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACG‑‑‑TTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAg                                                                                                                                                                                                 <  2:589508/236‑1 (MQ=255)
                                   tcttcATCAGTCAGGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAACTCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGAt                                                                                                                                                                                               <  2:416222/236‑1 (MQ=255)
                                    cttcATCAGTCAGGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATa                                                                                                                                                                                              >  1:91833/1‑236 (MQ=255)
                                    cttcATCAGTCAGGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATa                                                                                                                                                                                              <  2:12587/236‑1 (MQ=255)
                                    cttcATCAGTCAGGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTATTGTCCGCGCCTGCAAGGCGGGGGGGACCTTGcgcgta                                                                                                                                                                                              >  1:382715/1‑231 (MQ=255)
                                           agtcagGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTGTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTa                                                                                                                                                                                       <  2:659200/236‑1 (MQ=255)
                                                       gTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCggggg                                                                                                                                                                                                            <  1:674104/203‑1 (MQ=255)
                                                       gTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCggggg                                                                                                                                                                                                            >  2:674104/1‑203 (MQ=255)
                                                       gTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACCTTCTTCTTGTCCGCGCCTGCAAGTCGGGGGCGAGATTGGAGGTAATTGGTAACGCTAAGgcc                                                                                                                                                                            >  2:389477/1‑235 (MQ=255)
                                                        ttGCTTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGgcctg                                                                                                                                                                          <  2:487750/236‑1 (MQ=255)
                                                           cGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATa                                                                                                                                                                                              >  1:58628/1‑213 (MQ=255)
                                                           cGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATa                                                                                                                                                                                              <  2:58628/213‑1 (MQ=255)
                                                           cGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGt                                                                                                                                                                                        >  2:315860/1‑219 (MQ=255)
                                                           cGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGt                                                                                                                                                                                        <  1:315860/219‑1 (MQ=255)
                                                                 aaCAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAgg                                                                                                                                                                 >  1:736545/1‑236 (MQ=255)
                                                                                gTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCg                                                                                                                                                      <  1:758752/232‑1 (MQ=255)
                                                                                 tCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGcccc                                                                                                                                                 <  1:732337/236‑1 (MQ=255)
                                                                                      tAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGGGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGaa                                                                                                                                             >  2:227590/1‑235 (MQ=255)
                                                                                             gcattgtttcCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCAcc                                                                                                                                     <  2:35139/230‑1 (MQ=255)
                                                                                              cATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCAccc                                                                                                                                    >  1:35139/1‑236 (MQ=255)
                                                                                              cATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCAcc                                                                                                                                     >  1:208829/1‑235 (MQ=255)
                                                                                              cATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCAcc                                                                                                                                     <  2:208829/235‑1 (MQ=255)
                                                                                                  ttttaCGGGCAGTACATGCACAGAAAAGCCCAACCAAGGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCcagca                                                                                                                                <  2:279809/236‑1 (MQ=255)
                                                                                                  ttgtccGGGCATTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCcag                                                                                                                                  <  2:265477/229‑1 (MQ=255)
                                                                                                     ttCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCa                                                                                                                                       <  2:488661/226‑1 (MQ=255)
                                                                                                     ttCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCa                                                                                                                                       >  1:488661/1‑226 (MQ=255)
                                                                                                        gggCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTg                                                                                                                                                                                                          <  1:350629/156‑1 (MQ=255)
                                                                                                        gggCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTg                                                                                                                                                                                                          >  2:350629/1‑156 (MQ=255)
                                                                                                        gggCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGGACccaccaaggta                                                                                                                           >  2:54248/1‑230 (MQ=255)
                                                                                                           cAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTAATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAg                                                                                                                       <  1:421637/236‑1 (MQ=255)
                                                                                                                    cgcAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGGGATCTTGAAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGAc                                                                                                              >  1:222353/1‑236 (MQ=255)
                                                                                                                     gcAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACg                                                                                                             <  2:91833/236‑1 (MQ=255)
                                                                                                                     gaaGAAACGCCCAACCAAAGCTGGCTGCCAGTCGGTGGTTCCACCTGTCTTTTATTCCACAGGGCATCAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACg                                                                                                             <  2:382715/234‑1 (MQ=255)
                                                                                                                      cAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGGAATAGTCACCg                                                                                                                 >  1:394348/1‑231 (MQ=255)
                                                                                                                       aGAAAAGCCCAACCAAAGCGGGCTGCCAGTCGGTGGTTCCACATCTCTTTTATTCCACAGGGCAAAAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGa                                                                                                               <  2:394348/232‑1 (MQ=255)
                                                                                                                             gCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTa                                                                                                        <  1:227590/233‑1 (MQ=255)
                                                                                                                             gCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACg                                                                                                      <  2:48908/235‑1 (MQ=255)
                                                                                                                                  accaaAGCTCGCTGCCGGGCGGTGGTTCCACTTGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGtat                                                                                                  <  2:222353/234‑1 (MQ=255)
                                                                                                                                   ccaaAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGGCCCCGAAACCGCACCCAGCAGAGTAATAGTCACCGGACCCTGTAAGGtctaaa                                                                                               >  2:718438/1‑234 (MQ=255)
                                                                                                                                        gctggctgCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGGAATAGTCACCGGACGCTGTACGGTATAGAACgg                                                                                           >  1:748414/1‑234 (MQ=255)
                                                                                                                                                         gtttCCACATGTCTTCTATCCCACAGGGCAGCAGGGGCGACGACAGTCTCACCCGCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACg                                                                          <  2:664697/233‑1 (MQ=255)
                                                                                                                                                         ggtTCCGCAGGTCTTTTATGCCACAGGGCAACAGGGGCGAAGACAGTCTCACCCTCAGGCGTAATCAACAGAACGTTCTTCTTGGCCGCGCCTGCAAGGCGGGGGTGATCTGGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACg                                                                          <  2:157341/235‑1 (MQ=255)
                                                                                                                                                          gtTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCt                                                                                                           >  1:79632/1‑201 (MQ=255)
                                                                                                                                                          gtTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCt                                                                                                           <  2:79632/201‑1 (MQ=255)
                                                                                                                                                          gtTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACAGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGCATATTCCCCGACCaaacgg                                                                        >  2:435051/1‑232 (MQ=255)
                                                                                                                                                                caTGTCTTTTAGTCCACAGGTCACCAGGGGCGACGCCAGTCTCCCCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTg                                                                  <  2:748414/236‑1 (MQ=255)
                                                                                                                                                                    tCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCt                                                                                                           >  1:610435/1‑191 (MQ=255)
                                                                                                                                                                    tCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATaga                                                                                               <  1:81168/203‑1 (MQ=255)
                                                                                                                                                                    tCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATaga                                                                                               >  2:81168/1‑203 (MQ=255)
                                                                                                                                                                    gcTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCt                                                                                                           <  2:610435/190‑1 (MQ=255)
                                                                                                                                                                     cTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCGCCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGtt                                                             <  2:285554/236‑1 (MQ=255)
                                                                                                                                                                                caGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTc                                                  >  1:758097/1‑236 (MQ=255)
                                                                                                                                                                                caGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGT‑‑CCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGt                                                >  2:338647/1‑236 (MQ=255)
                                                                                                                                                                                   ggCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTa                                               >  1:207146/1‑236 (MQ=255)
                                                                                                                                                                                             gCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAg                                      >  1:112235/1‑235 (MQ=255)
                                                                                                                                                                                             gCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGGATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTTTCATCCACTCGACCAAAACCac                                      >  2:59152/1‑234 (MQ=255)
                                                                                                                                                                                               gacgacAGTCTCACCCTGAGGCGTAAGCACCATAACGTTCTTCTTGGCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAggg                                    <  2:207146/235‑1 (MQ=255)
                                                                                                                                                                                                acgacAGTCTCACCCTCAGGCGTAATCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGGAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCAACCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAg                                      >  1:153062/1‑232 (MQ=255)
                                                                                                                                                                                                                       aaTCACCATAACGTTCTTCTTGTCCGTGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGATACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCCCCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGAc              <  2:426862/233‑1 (MQ=255)
                                                                                                                                                                                                                       aaTCACCATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGAc              <  2:112235/233‑1 (MQ=255)
                                                                                                                                                                                                                             cATAACGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGatat     >  1:214868/1‑236 (MQ=255)
                                                                                                                                                                                                                                gaCGTTCTTCTTGGCCGCGCCGGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATattga  <  1:59152/235‑1 (MQ=255)
                                                                                                                                                                                                                                aaCGTTCTTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATattga  <  1:435051/236‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                     
GGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGA  >  NC_000913/1047904‑1048363

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: