Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,356,445 A→T 5.1% intergenic (‑34/+25) dtpC ← / ← cadA dipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,356,4450AT 5.1% 268.3 / 4.9 79intergenic (‑34/+25)dtpC/cadAdipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible
Reads supporting (aligned to +/- strand):  ref base A (39/36);  new base T (2/2);  total (41/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 = 2.51e-02

TCACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTC  >  NC_000913/4356177‑4356718
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                 
ccactgagtgggggTTGCCGAGCAGGCGGTCGGAAAGCCAGCCGCCGGGAATAGGGGGAACGTAAACCAGGGAATCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACGGAAGTACTCCCGGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCAc                                                                                                                                                                                                                                                                             <  4:36147/263‑1 (MQ=255)
           cGGTTGCCGGGCAGGCGGTCGGCAAGCCAGACGCCCAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTa                                                                                                                                                                                                                                                                   <  4:141782/274‑1 (MQ=255)
                                                gaATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCtt                                                                                                                                                                                                                             <  2:427245/275‑1 (MQ=255)
                                                gaATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCt                                                                                                                                                                                                                              <  1:157471/274‑1 (MQ=255)
                                                           aCGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCCTGTGTTCTCCTTATGAGCAAACAAGGGAAGTGGCAAGCCACTTCACTTGTACGCGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAg                                                                                                                                                                                                                   >  1:151641/1‑274 (MQ=255)
                                                           aCGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGt                                                                                                                                                                                                                       >  4:146779/1‑270 (MQ=255)
                                                           aCGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGGGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTGATTATTTTTTGCGTTCTTCTTTCACTACCTTAACGGTATAGc                                                                                                                                                                                                                  >  2:259089/1‑275 (MQ=255)
                                                               aaaCCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAg                                                                                                                                                                                                                   <  2:146780/270‑1 (MQ=255)
                                                                aaCCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAACAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTGCAATACCTTAACGGTATAGCGACCa                                                                                                                                                                                                             >  1:151175/1‑275 (MQ=255)
                                                                        aaGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATGGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTg                                                                                                                                                                                                     <  1:60086/275‑1 (MQ=255)
                                                                             tatGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTCCTTCTTTCAATAcc                                                                                                                                                                                                                               >  2:431483/1‑244 (MQ=255)
                                                                                                tGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTGCCCTTGTACGATCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCCCCGTGAATATCGGt                                                                                                                                                                             >  2:326019/1‑275 (MQ=255)
                                                                                                                      aaGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGta                                                                                                                                                                                                                      <  1:5776/212‑1 (MQ=255)
                                                                                                                       aGCTGTTGGGTGGGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAgcgc                                                                                                                                                      <  2:390704/275‑1 (MQ=255)
                                                                                                                            aTGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttt                                                                                                                                                                                                                                                    <  4:213319/176‑1 (MQ=255)
                                                                                                                            aTGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATCTTGAGCGCCGAt                                                                                                                                                  >  1:145036/1‑274 (MQ=255)
                                                                                                                                gTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACt                                                                                                                                                                                                                                                                            <  7:644709/148‑1 (MQ=255)
                                                                                                                                   agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTcaca                                                                                                                                            >  1:20864/1‑273 (MQ=255)
                                                                                                                                        gggaATGATTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCATATTTGGATCGCCACGATATAGTAGATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCa                                                                                                                                       <  4:300823/271‑1 (MQ=255)
                                                                                                                                            aatggTAAGGCACGAATGGCGTAAAAACTGAAGGACTCCCAGATTTGGGTGGCCACGAGAGATTATGTCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTGATGAGCGAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTGATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGc                                                                                                                                  <  4:437066/270‑1 (MQ=255)
                                                                                                                                            aTGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCa                                                                                                                                                                                                                                                                                   >  5:156291/1‑129 (MQ=255)
                                                                                                                                            aTGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCa                                                                                                                                                                                                                                                                                   <  8:740019/129‑1 (MQ=255)
                                                                                                                                            aTGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCa                                                                                                                                                                                                                                                                                   <  6:156291/129‑1 (MQ=255)
                                                                                                                                            aTGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCa                                                                                                                                                                                                                                                                                   >  7:740019/1‑129 (MQ=255)
                                                                                                                                               aggaaGGCACGCTTGCCTTAANNATTGANTTTCTCCCAGATTTGGTTCGCCACGGTGTAGGATATCGCGCGCTGCTGGTGGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGTGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAgg                                                                                                                               <  3:18926/271‑1 (MQ=255)
                                                                                                                                               aagaaGGCAGCCATGCTTGAAAAAGTGANNNCCTCCCAGTTTTGGTTCGGCACGATATAGTTGATCGGGCGGGGTTGTGAGGGTGTGTTAAGGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTACCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCGATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAgg                                                                                                                               <  3:145036/271‑1 (MQ=255)
                                                                                                                                                    ggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCa                                                                                                                                                                                                                                                                                   >  6:394080/1‑121 (MQ=255)
                                                                                                                                                    ggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCa                                                                                                                                                                                                                                                                                   <  5:394080/121‑1 (MQ=255)
                                                                                                                                                        cGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGcc                                                                                                                                                                                                                                                                               <  8:345948/121‑1 (MQ=255)
                                                                                                                                                        cGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGcc                                                                                                                                                                                                                                                                               >  7:345948/1‑121 (MQ=255)
                                                                                                                                                         gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGaa                                                                                                                     >  2:15171/1‑274 (MQ=255)
                                                                                                                                                              ccGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCtt                                                                                                                                                                                                                                              >  2:415669/1‑148 (MQ=255)
                                                                                                                                                               cGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTg                                                                                                                                   >  2:518719/1‑254 (MQ=255)
                                                                                                                                                               cGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAACATCTGCAGGAACTCCAGAACAGGACg                                                                                                              >  1:4440/1‑275 (MQ=255)
                                                                                                                                                                  aaaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttctt                                                                                                                                                                                                                                          <  4:415668/148‑1 (MQ=255)
                                                                                                                                                                   aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAgg                                                                                                                               <  4:518718/254‑1 (MQ=255)
                                                                                                                                                                                      aTTTGTATCGCCACGATATAGTATTGCGCGCGCGGCTTGGAGGGTGTTTGCATGTGTTCTCCTTATGAGGAAAAAAGGGAAGAGGCAAGCCACTGCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCAc                                                                                       <  1:18928/275‑1 (MQ=255)
                                                                                                                                                                                                ccACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAg                                                                                                                                                                                                                                                                                 >  6:779852/1‑79 (MQ=255)
                                                                                                                                                                                                ccACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCgg                                                                                                                                                                              >  2:198809/1‑178 (MQ=255)
                                                                                                                                                                                                ccACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTAc                                                                             >  2:411867/1‑275 (MQ=255)
                                                                                                                                                                                                ccACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCCCCTTATGAGCAAAAAAGGGAAGTGGCAAg                                                                                                                                                                                                                                                                                 <  5:779852/79‑1 (MQ=255)
                                                                                                                                                                                                 cACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAg                                                                                                                                                                                                                                                                                 >  7:257417/1‑78 (MQ=255)
                                                                                                                                                                                                 cACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAg                                                                                                                                                                                                                                                                                 <  8:257417/78‑1 (MQ=255)
                                                                                                                                                                                                  aCGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTc                                                                                                                                                                                                                                                                          <  5:1036931/84‑1 (MQ=255)
                                                                                                                                                                                                  aCGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTc                                                                                                                                                                                                                                                                          >  6:1036931/1‑84 (MQ=255)
                                                                                                                                                                                                    gATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTc                                                                                                                                                                          <  4:198808/178‑1 (MQ=255)
                                                                                                                                                                                                              atCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACgg                                                                                                                                                                                                 >  4:117420/1‑145 (MQ=255)
                                                                                                                                                                                                               tCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACg                                                                                                              <  2:466263/227‑1 (MQ=255)
                                                                                                                                                                                                                cgcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAg                                                                                                                                                                                                                                                                                 >  4:539133/1‑63 (MQ=255)
                                                                                                                                                                                                                cgcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTtcttct                                                                                                                                                                                                                                        >  4:551498/1‑104 (MQ=255)
                                                                                                                                                                                                                cgcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTcttctt                                                                                                                                                                                                                                       >  3:80918/1‑105 (MQ=255)
                                                                                                                                                                                                                cgcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTcttctt                                                                                                                                                                                                                                       >  3:93450/1‑105 (MQ=255)
                                                                                                                                                                                                                cgcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAg                                                                                                                                                                                                         >  1:81475/1‑135 (MQ=255)
                                                                                                                                                                                                                cgcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTTAAGGAACTCCAGAACCGGCCGGCTTTCTTCCGTGATCATTTCTCCCAGCATTACAACAGGAACTCCcggc                                                             >  1:133883/1‑275 (MQ=255)
                                                                                                                                                                                                                cgcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTTCcgg                                                              >  2:178922/1‑274 (MQ=255)
                                                                                                                                                                                                                cgcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCcgg                                                              >  2:412791/1‑274 (MQ=255)
                                                                                                                                                                                                                  cgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATg                                                                                                                                                                                             <  2:117421/145‑1 (MQ=255)
                                                                                                                                                                                                                  cgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTc                                                                                                                                                                          >  1:23680/1‑164 (MQ=255)
                                                                                                                                                                                                                  cgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCg                                                                                                                  >  1:155510/1‑220 (MQ=255)
                                                                                                                                                                                                                  cgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGCCGGCTTTCTTCGGTGCTCATTTCACCCGGCATTACCAGAGGAACTCCcggcgg                                                           >  1:158615/1‑275 (MQ=255)
                                                                                                                                                                                                                  cgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCcggcgg                                                           >  2:149873/1‑275 (MQ=255)
                                                                                                                                                                                                                    cgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCAc                                                                                                                                                                                                                                                                             <  2:539134/63‑1 (MQ=255)
                                                                                                                                                                                                                    cgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTGCGAGCTAATTATTTTTTGCTTTCTTCTTTCAATGCCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCCGAGGAACTCCCGGCGGGt                                                         >  2:514741/1‑275 (MQ=255)
                                                                                                                                                                                                                    cgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCaa                                                                                                                                                                                                                                   <  1:93450/105‑1 (MQ=255)
                                                                                                                                                                                                                    cgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCaa                                                                                                                                                                                                                                   <  1:80918/105‑1 (MQ=255)
                                                                                                                                                                                                                    cgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCa                                                                                                                                                                                                                                    <  2:551499/104‑1 (MQ=255)
                                                                                                                                                                                                                    cgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTg                                                                                                                                                                                                     <  3:81475/135‑1 (MQ=255)
                                                                                                                                                                                                                    cgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGt                                                         >  2:55445/1‑275 (MQ=255)
                                                                                                                                                                                                                    cgcgGCTGTGAGGGTGTCTCCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGCGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGGATAGCGGCCATCAGCCTGACGGTTTGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTCCACATAGCATCTGCAGGAACTACAGAATCGGCAGGCTTTCTTCGGTGATCATTTCACCCGGCACTACCAGAGGAACTCCAGGCGGGt                                                         >  4:139679/1‑275 (MQ=255)
                                                                                                                                                                                                                      cgGCTGTGAGGGTGTTTTCATGTGTNCTCCTTNTGANCAAAAAANGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAg                                                                                                                                                                      <  3:23680/164‑1 (MQ=255)
                                                                                                                                                                                                                         gtgtgCGGGTGTTTTCATGTGTTCTCCTGATGAGCAAAAAAGGGAATGGGAAATCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGa                                                    <  4:326018/274‑1 (MQ=255)
                                                                                                                                                                                                                          tgtgGGGGGGTTTTCAAGGGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACgg                                                     <  4:411866/273‑1 (MQ=255)
                                                                                                                                                                                                                             gAGTGTTTTTTCATGGTTTCACCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGa                                                <  4:178921/275‑1 (MQ=255)
                                                                                                                                                                                                                                      ttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCAc                                                                                                                                                                                                                                                                             >  7:26314/1‑45 (MQ=37)
                                                                                                                                                                                                                                      ttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAANCGGACGCCTTTCTTCGGTGAACATTTCACCCGCCCGTACCACAGGACTCCCCGGGAGCTGCGGACGGATCATCttt                                        >  3:117959/1‑273 (MQ=255)
                                                                                                                                                                                                                                              ttCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATa                               >  2:320735/1‑275 (MQ=255)
                                                                                                                                                                                                                                                             aaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAg                                                                                                                                                                                                         >  1:106629/1‑90 (MQ=255)
                                                                                                                                                                                                                                                             aaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTg                                                                                                                                   >  4:297463/1‑160 (MQ=255)
                                                                                                                                                                                                                                                               aaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGa                                                                                                                                                   <  6:537700/142‑1 (MQ=255)
                                                                                                                                                                                                                                                               aaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGa                                                                                                                                                   >  5:537700/1‑142 (MQ=255)
                                                                                                                                                                                                                                                                 aGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTg                                                                                                                                                                                                     <  3:106629/90‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAgg                                                                                                                               <  2:297464/160‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCAtt                                                                                           >  2:224804/1‑191 (MQ=255)
                                                                                                                                                                                                                                                                          gCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCAc                                                                                       <  4:224803/191‑1 (MQ=255)
                                                                                                                                                                                                                                                                           cAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCCTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAActtc  <  2:409186/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                 
TCACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTC  >  NC_000913/4356177‑4356718

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: