Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,516,171 G→C 100% N169K (AAC→AAG)  fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,516,1710GC100.0% 223.1 / ‑6.7 60N169K (AAC→AAG) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (29/30);  minor base T (0/1);  total (29/31)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCT  >  NC_000913/4515925‑4516400
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                     
ggCACGGGTAACAAAGTTCCCCACGCCGCCCACGCTCTGCGGTCCGCAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAg                                                                                                                                                                                                              <  1:434914/272‑1 (MQ=255)
    cgtgtcaCAAAGTTCACCACGCCGCCCACGCTCTGCGGGCCGTAACGCACCGCACCACGACCGCGTACCACGTCAATGGCATTCATGTTGCCGAGCGAAACGGGAGCCAGTGGAAGCGTCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCAGCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAggg                                                                                                                                                                                                            <  1:486138/264‑1 (MQ=255)
    cGGGTAACAAAGTTCCCCACGCCGCCCACGCTCTGCGGTCCGTACCGCACCGCACCCCCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGATTTCCGGCGCGCTGACGCCAggg                                                                                                                                                                                                            <  1:436121/270‑1 (MQ=255)
                    ccacGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAggg                                                                                                                                                                                                            <  1:690429/254‑1 (MQ=255)
                                                         accaccaccGCGTAACACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGGTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGGCGTGGCTGCCGGGGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCCCACGCATTGTGGTTGCGCCGGTTTTGGGGAAATCCTcc                                                                                                                                                      >  1:1096053/1‑270 (MQ=255)
                                                                   cgTACCACGTCAATGGCATCCCTGTTGCCGAGCGAAACGGGCGCCCGTGAAAGCTGCGGCTGACCGTAAGGGGCGCAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGCCTGGCGAGGCGCGGGCTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGgcg                                                                                                                                                                                                                          >  1:1068576/1‑193 (MQ=255)
                                                                      accgcggaaaTGGCGTCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCGGCGGCTGACCGGAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGAtcacgtc                                                                                                                                        <  2:759824/264‑1 (MQ=255)
                                                                      aCCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTTACGCCAGGGGTGCGGTTAAGTACCTCAAGCATGGTGGTTGCGCCGGGTTTGGCGAAATCCTCACGGCGGGtcacgt                                                                                                                                         >  1:759825/1‑271 (MQ=255)
                                                                      aCCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGGTGCGGGTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGAtcacggc                                                                                                                                        >  1:154333/1‑270 (MQ=255)
                                                                      aCCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTAACTCACGCATGGTGGGTGCGCCGGGTTTGGCGAAATCCTTCAGGCGGacgcc                                                                                                                                           >  1:997424/1‑265 (MQ=255)
                                                                       ccACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGGGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGgcgc                                                                                                                                                                                                                         <  2:706447/190‑1 (MQ=255)
                                                                             cAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCTTCAGCGCGGATCAGGTCAcgcgc                                                                                                                                  >  1:288846/1‑271 (MQ=255)
                                                                                        aTGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCCTGGTGGTTGCGCCGGGTTTGGCGAAATCCTCACGGCGCATCACGCCCCGCGCGCCCGCAGGtgc                                                                                                                      >  1:54501/1‑270 (MQ=255)
                                                                                                       aCGGGGGCCGGGGAACGCGGCGGCGTACCGGAAGGGGGGAAGGGGCCGGGGAGGCCGTCCATCAGGACGGGCGGGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGtttt                                                                                                        <  2:112843/271‑1 (MQ=255)
                                                                                                                   gAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAgg                                                                                                                                                                                                             >  1:184532/1‑158 (MQ=255)
                                                                                                                    aaaGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAggg                                                                                                                                                                                                            <  2:184532/158‑1 (MQ=255)
                                                                                                                                  ccGTAAGGGGCGGGTGGTGCGGGGTGGCCCTCCATCAGGACGGTCGGGCGGCGGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCCTCGCCAGGTCGTGGCTGCCGGTGCGCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGTTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGa                                                                             <  1:433866/271‑1 (MQ=255)
                                                                                                                                    gTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAgg                                                                                                                                                                                                             >  1:667900/1‑141 (MQ=255)
                                                                                                                                     tAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAggg                                                                                                                                                                                                            <  2:667899/141‑1 (MQ=255)
                                                                                                                                       gtgggCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCGGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGTGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACg                                                                        <  1:104013/269‑1 (MQ=255)
                                                                                                                                              aaGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCt                                                                                                                                                        >  1:274853/1‑184 (MQ=255)
                                                                                                                                              aaGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGc                                                                >  1:151040/1‑272 (MQ=255)
                                                                                                                                               aGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTc                                                                                                                                                       <  2:274852/184‑1 (MQ=255)
                                                                                                                                                       gggATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAgg                                                                                                                                                                                                             >  2:267808/1‑122 (MQ=255)
                                                                                                                                                        ggATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTc                                                                                                                                                                                                                               >  1:588640/1‑103 (MQ=255)
                                                                                                                                                        ggATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAggg                                                                                                                                                                                                            <  1:267809/122‑1 (MQ=255)
                                                                                                                                                         gATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTcc                                                                                                                                                                                                                              <  2:588639/103‑1 (MQ=255)
                                                                                                                                                              cGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGa                                                                                                                                                                                                                    >  1:422147/1‑108 (MQ=255)
                                                                                                                                                              cGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCt                                                                                                                                                                                            >  2:660810/1‑132 (MQ=255)
                                                                                                                                                               gTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGAc                                                                                                                                                                                                                   <  2:422146/108‑1 (MQ=255)
                                                                                                                                                               gTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTc                                                                                                                                                                                           <  1:660811/132‑1 (MQ=255)
                                                                                                                                                                tCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAgg                                                                                                                                                                                                             >  2:131448/1‑113 (MQ=255)
                                                                                                                                                                 ccATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAggg                                                                                                                                                                                                            <  1:131448/113‑1 (MQ=255)
                                                                                                                                                                     cAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAgg                                                                                                                                                                                                             >  2:258190/1‑108 (MQ=255)
                                                                                                                                                                      aGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAggg                                                                                                                                                                                                            <  1:258190/108‑1 (MQ=255)
                                                                                                                                                                             tCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGAGGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGGTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTGCAAATTAGTCGTTTTCACGGGCCTCACCCAGCCAGTCGCCGCCCACGGTCAGGGCATCTTCTTTTTGTTGGGGCGCGGGCGCCAga                                 >  1:643548/1‑271 (MQ=255)
                                                                                                                                                                               gAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAgg                                                                                                                                                                                                             >  1:67173/1‑98 (MQ=255)
                                                                                                                                                                                aGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGgcg                                                                                                                                                                                                                          >  1:215120/1‑84 (MQ=255)
                                                                                                                                                                                aGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAggg                                                                                                                                                                                                            <  2:67173/98‑1 (MQ=255)
                                                                                                                                                                                aGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAgg                                                                                                                                                                                                             >  1:717051/1‑97 (MQ=255)
                                                                                                                                                                                aGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATg                                                                                                                                                                                                         >  1:245589/1‑101 (MQ=255)
                                                                                                                                                                                aGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGt                                                                                                                                                                                                >  1:773308/1‑110 (MQ=255)
                                                                                                                                                                                aGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACACCAGGGATGCGGGTAAGTACCTCACGCGTGGTGGTTGCGCCGGGTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCCCATCACCCAGCCAGTCGc                                                                                >  1:1055195/1‑222 (MQ=255)
                                                                                                                                                                                 gCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGgcgc                                                                                                                                                                                                                         <  2:215120/84‑1 (MQ=255)
                                                                                                                                                                                 gCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAggg                                                                                                                                                                                                            <  2:717050/97‑1 (MQ=255)
                                                                                                                                                                                 gCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGc                                                                                                                                                                                                        <  2:245589/101‑1 (MQ=255)
                                                                                                                                                                                 gCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTa                                                                                                                                                                                               <  2:773307/110‑1 (MQ=255)
                                                                                                                                                                                 gCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTcc                             <  1:596263/272‑1 (MQ=255)
                                                                                                                                                                                          gAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCCGCCAGTCGCCGACCACGGTCAggg                                                                 >  2:826210/1‑227 (MQ=255)
                                                                                                                                                                                           aGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGc                                                                <  1:826211/227‑1 (MQ=255)
                                                                                                                                                                                                     tCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCCCGGGCAGGGCATCTTCCTTTGGTTCGGGCGCGGGCTCCGTCCCCCAGCCGTTAGTTCCCacgggt         >  2:796896/1‑267 (MQ=255)
                                                                                                                                                                                                             gggtggCCTAAGCTCTTCGTCGGTTCGTGGGTGCCGGGGCCTTGGTGTTCCGGCGCGCTGACGCCGGGGAGGCGGTTAAGTACCTCACGCATGGTGGTGGCGCCGGTTTTGGCGAAATCCTCACTGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCt  <  2:997423/266‑1 (MQ=255)
                                                                                                                                                                                                                aTGCCAAAGTTCATCGCCAGGTAGTGGCTGCCGGTGGCCTTGGTTTCCGgcg                                                                                                                                                                                                                          >  1:372007/1‑52 (MQ=255)
                                                                                                                                                                                                                 tGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGgcgc                                                                                                                                                                                                                         <  2:372006/52‑1 (MQ=255)
                                                                                                                                                                                                                 tGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCAc                                                                                                     >  1:380777/1‑168 (MQ=255)
                                                                                                                                                                                                                 tGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTgcgggcgcggg                                        <  1:350018/229‑1 (MQ=255)
                                                                                                                                                                                                                 tGCCAAAGTGCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTcc                <  1:613350/253‑1 (MQ=255)
                                                                                                                                                                                                                  gCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCAcg                                                                                                    <  2:380776/168‑1 (MQ=255)
                                                                                                                                                                                                                                           cTGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAg                                                                   >  2:189263/1‑176 (MQ=255)
                                                                                                                                                                                                                                            tGCCGGTGCCCTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAgg                                                                  <  1:189263/176‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                     
GGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCT  >  NC_000913/4515925‑4516400

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