Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,514,559 T→C S707G (AGT→GGT)  fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,514,5590TC97.3% 108.7 / ‑6.1 37S707G (AGT→GGT) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (19/17);  minor base G (0/1);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTAC  >  NC_000913/4514317‑4514795
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                            
gccAGCAAAAAGGCCGGGAAGCGCGGCGTTAACACAGCAGAACTGTAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGTTGAAGTAGTCCTGGTCGAAGATGTTTGTCACACCGAAGGCCAGGTTTAGAGCTGTCATCTGGGGGCCAAAGTCGGACGCCACGCGTGCGGCCCAGGGGATGAAGCCGGGGATGGGGCCGGTACCGCCGTCgg                                                                                                                                                                                                                                      <  2:99543/250‑1 (MQ=255)
  cAGCTAAAAGACCGGCATGCCGGGCGTTCACACATCAGAACTTCAACGAACCCTGCATAAACAGCGTGCGCTGCTGGCATGCATAGCTGCATTTGTTGTTGTCGTCATAAGAGAGGCTGATGTAGTCCTGGTCGTAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCCCGCGTGAGCCCCAGAGCATGAAGCCCGGACTGCGGCCGGTACCGCCGTCGgcc                                                                                                                                                                                                                                    >  2:96915/1‑250 (MQ=255)
       aaaaGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGctc                                                                                                                                                                                                                                 <  2:212616/249‑1 (MQ=255)
            cccGGCACGCCGGGCGTTAACCAATCAGAACTTCAACGACCCCTGCATAAACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGAGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGAACCGCGTGCGCCCCAGAGCAGGAAGCCGGGAATGCGGCCGGTCCCGCCGTCGGCGATCTCTTTcc                                                                                                                                                                                                                          >  2:206010/1‑250 (MQ=255)
            cccGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGCATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTc                                                                                                                                                                                                                           >  1:8974/1‑250 (MQ=255)
               ggCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATACGAGCGGATGAAGTAGTCCTGGCCGAAGACGATTTTCACACCGACTGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATCCGCCACGCGGGCGCCCCAGAGCATGAAGCCGGGATTGCGGCCGGTCCCGCCGGCGGCGCTCTCTTTCCCCg                                                                                                                                                                                                                       >  2:244057/1‑251 (MQ=255)
                 cAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGGCGTCATAAGAGCGGAGGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTa                                                                                                                                                                                                                     <  2:33256/251‑1 (MQ=255)
                      cGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGCTGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGCATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGc                                                                                                                                                                                                                >  1:98757/1‑251 (MQ=255)
                      cGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGc                                                                                                                                                                                                                >  1:263352/1‑251 (MQ=255)
                             tAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATcc                                                                                                                                                                                                         >  1:120184/1‑251 (MQ=255)
                                 acaTCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCAGCTGGCCTGCATAGATGCCTTTGTTGTTGTAGTCCTAAGAGCGGATGCAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCa                                                                                                                                                                                                      >  2:27744/1‑250 (MQ=255)
                                 acaTCAGAACTTCAACGACCCCTGCAAATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTGTGTGGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGAGCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGTGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCaa                                                                                                                                                                                                     <  2:98757/251‑1 (MQ=255)
                                      aGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAActgg                                                                                                                                                                                                >  1:235769/1‑251 (MQ=255)
                                               acgtccccTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGGCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGctggactgg                                                                                                                                                                                       <  2:235769/247‑1 (MQ=255)
                                                    cccTGCATATACAGCGTGCGGGTCTGGCTTGCATAGATGCCTTTGTCGTTTTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATTTGCGGGCCAAAGTCATACGCCACGCGTGCGCCTCATAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTATTACCGTATTAGCGTTATCCGTAAACTGGCTGGAGTGGAAATc                                                                                                                                                                                  <  2:138230/251‑1 (MQ=255)
                                                               cAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTc                                                                                                                                                                           >  1:260894/1‑247 (MQ=255)
                                                                             ggCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGtt                                                                                                                                                          <  2:263352/250‑1 (MQ=255)
                                                                                     aagaggCCTGTGTTGTTGTCGTCAGAAGAGCGGATGAATTAGGCCTGGTCGAAGATGTTTTTCACCCCGAATGCCAGGTTCAGATCTGCAATCTGCGGCCCAAAGTCATACGCCACGCGTGCGCCCCAGATCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCTCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGATCGTCCAGTTTCCTGGCtt                                                                                                                                                 <  2:8974/246‑1 (MQ=255)
                                                                                         tGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGt                                                                                                                                               >  1:159730/1‑249 (MQ=255)
                                                                                         tGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGt                                                                                                                                               >  2:208948/1‑249 (MQ=255)
                                                                                                           cATAAGAGCGGATGAAGTAGTCCTGGTTGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGGGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCTGTACGGCCGTCGGCGCGCTGGTTCACCGTATTGGTGTTATCCGCAAACTGGCTGGACTGGAAATCGCGGTGCAGATGGAGCGTGCAGTGTCCTGGCTTGTAGTCCACGCCCAGCGTGcc                                                                                                                            <  1:28440/250‑1 (MQ=255)
                                                                                                               agagCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTAt                                                                                                                       >  1:23770/1‑251 (MQ=255)
                                                                                                                    ggATGCAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAAAGCCAGGTTCAGATCTGCCATCTGCGGGGCAAAGACATACGCAACGCGTGAGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGGACCGCCGTCGGCGATCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATAGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCAAGCCCAGCGTGCCTTTCTGTTTc                                                                                                                  >  2:257250/1‑251 (MQ=255)
                                                                                                                              gtCCGGGTCGAAGATGGGTTTCACACCGAATGCCAGGTTCCGACCTGTCATCTGCGGGCCAAAGTCATACGCCACGCGGGCGCCCCAGAGCATGAAGCCGGTAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTTTTCAGATTGAACGGCCAGGTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTGTATGTTTCGGTGAGAATg                                                                                                        <  2:211214/251‑1 (MQ=255)
                                                                                                                                 cTGGTCGAAGATGGTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGAGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGACGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCACACTGGACGGAATGGAAATCGCTGTTCATATTGAACGTCAAAGTTCCTGTATTGTAGTCCCCACCCACATTGCCTTTATTTTTCgggggaaagggt                                                                                                      >  2:278234/1‑242 (MQ=255)
                                                                                                                                              tttttCACACCGAATGCCAGGGTCAGATCTGCCATCTGCGGGCCAAAGTCATAAGCCACGCGTGCGCCCCAGCGCATGAAGCCGGGAATGCGGCCGGTACCGCAGTCGGCGCTCTCTTTCACCGTATGGGCGTTATCCGCAAACTGGCTGGCCTGGAAAACGCTGTTCAGATTGTAAGTCCAGTATCCTGGCTTATTGTACACGCCCAGCGTGCCTTTATGTTTTGGGGAGCATGGTACCAAGTTGCCgg                                                                                         >  2:172201/1‑249 (MQ=255)
                                                                                                                                                     cacCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCGCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCg                                                                                 <  2:120184/251‑1 (MQ=255)
                                                                                                                                                             gCCAGGTTCAGATCTGCCGTCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAATTCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTtct                                                                          <  1:270737/250‑1 (MQ=255)
                                                                                                                                                             gCCAGGTTCAGATCTGCCACCGGCGGGCCAAAGACATACGCCACGCGGGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGTCGCTCTCTTTCACCGTAGTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCGTGTAGTCCACGCCCAGCTTGTCTTTAGGTTTCGGGGAGAATGGTACCAGATGGCCGTAGGTGTCGCCTTtct                                                                          <  2:28898/250‑1 (MQ=255)
                                                                                                                                                                aGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTcgcg                                                                      <  1:208948/251‑1 (MQ=255)
                                                                                                                                                                              cATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTTCCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTAt                                                                                                                       <  1:106101/188‑1 (MQ=255)
                                                                                                                                                                              cATCTACGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTGTTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTCGTCCACGCCCAGCGTGCCTTTAt                                                                                                                       >  2:106101/1‑188 (MQ=255)
                                                                                                                                                                                           aaGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCt                                           <  1:27744/251‑1 (MQ=255)
                                                                                                                                                                                                          cgTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTCCCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCCGCGTGCCTTTATGTTTCGGGGAGAATGGGACCAGATTGCCGTATGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATAAGCATAGCTGGAGTAGATGGCAAc                            >  2:46615/1‑251 (MQ=255)
                                                                                                                                                                                                           gTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACg                           <  1:178723/251‑1 (MQ=255)
                                                                                                                                                                                                           gTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACg                           >  1:25836/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aaTGCGGCCGGTACCGCCGTCGGCGCTCTCTTTCACCGTAGTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTAc  <  1:40677/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                            
TCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTAC  >  NC_000913/4514317‑4514795

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