Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,513,841 A→G R156R (CGT→CGC fecB ← ferric citrate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,513,8410AG98.4% 199.5 / NA 64R156R (CGT→CGCfecBferric citrate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (32/30);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

GTGGGCAACCAGCAGCCAGGCAGGATTGACCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTT  >  NC_000913/4513601‑4514074
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          
gggggCAACCGGCAGCCGGGCAGGATTGACCGCCAGCGGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCGCCTGATTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCa                                                                                                                                                                                                                                 <  2:74980/249‑1 (MQ=255)
     cAACCAGCAGCCAGGCAGGATTGACCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCctctct                                                                                                                                                                                                                            >  1:230023/1‑251 (MQ=255)
                 gggcagtATTGACCGCCAGCTGTTGCTCCGGGCCGATGGACGGCATGGACGCCCCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCGCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTAccc                                                                                                                                                                                                                <  2:396384/250‑1 (MQ=255)
                   gcaggATTGACCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACccac                                                                                                                                                                                                              >  1:167762/1‑251 (MQ=255)
                      gtatgacccGCCACCAGTTGCCCCAGGCCGATGGACGGCTTGGTCTCTCCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTTATGTGCCAAAGGCCTCGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCAt                                                                                                                                                                                                           <  1:43314/244‑1 (MQ=255)
                        aTTGACCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGGTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCGCTCGCTTTTTACCCACCAttt                                                                                                                                                                                                         >  1:348620/1‑251 (MQ=255)
                           gACCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCg                                                                                                                                                                                                       <  1:330500/250‑1 (MQ=255)
                           gACCGCCAGCAGTTGCTCCAGGCCGAGGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCg                                                                                                                                                                                                       <  1:316028/250‑1 (MQ=255)
                             ccGCCAGCAGTTGCTCCAGGCCGAGGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTACGTTGTTCCAGGCGTGCCTGCATCCCTCGCTTTTTACCCACCATTTCGCCg                                                                                                                                                                                                    <  1:417231/251‑1 (MQ=255)
                                     agTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGc                                                                                                                                                                                            >  1:158941/1‑251 (MQ=255)
                                      gtggCTCCAGGCCGATGGCCGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGcc                                                                                                                                                                                           <  1:330377/248‑1 (MQ=255)
                                      gTTGCTCCAGGCCGATGGACGGCATGGACGCCCCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGcc                                                                                                                                                                                           <  1:142157/251‑1 (MQ=255)
                                          cTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGc                                                                                                                                                                                         >  1:269749/1‑249 (MQ=255)
                                            ccAGGCCGGTGTACGGCAGTGACCCCCCCGCCATCGCAGCGGGAACTTTCAGCCCCATAGAGGCCAGCACGCTGCCGGTCCAGTTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTTCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGAt                                                                                                                                                                                     <  2:437254/251‑1 (MQ=255)
                                             cAGGCCGATGGCCGGCATGGACCCCCCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGAtt                                                                                                                                                                                    <  1:369672/251‑1 (MQ=255)
                                             cAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGAtt                                                                                                                                                                                    >  2:472380/1‑251 (MQ=255)
                                              aGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCCGATTg                                                                                                                                                                                   >  2:196427/1‑251 (MQ=255)
                                                gCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGACTTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGc                                                                                                                                                                                  >  2:435766/1‑250 (MQ=255)
                                                                gACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCCAATTTTCCGCGTAggg                                                                                                                                                                 >  2:191976/1‑250 (MQ=255)
                                                                 ggacACCCGCTAGCGCCGCGGTAACGTGCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGtt                                                                                                                                                                <  2:401565/248‑1 (MQ=255)
                                                                 ccGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGtt                                                                                                                                                                <  1:196427/250‑1 (MQ=255)
                                                                     aCCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGt                                                                                                                                                            >  1:4252/1‑251 (MQ=255)
                                                                         gCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACTCGGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTGTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCg                                                                                                                                                        <  2:268901/251‑1 (MQ=255)
                                                                            aTCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCggg                                                                                                                                                      >  2:374372/1‑250 (MQ=255)
                                                                                cAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTa                                                                                                                                                 <  2:100928/251‑1 (MQ=255)
                                                                                  gCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAag                                                                                                                                               <  1:268607/251‑1 (MQ=255)
                                                                                      gAACGTTCAACCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGCGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTAACCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGTCTTAagcagc                                                                                                                                           >  2:472776/1‑251 (MQ=255)
                                                                                        aCGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGCTGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATCGCCGCAGTTTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGGCTTAagaagc                                                                                                                                           >  2:284954/1‑249 (MQ=255)
                                                                                            tCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTAc                                                                                                                                      <  1:191976/250‑1 (MQ=255)
                                                                                               gCCCCAGAGATGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGgc                                                                                                                                  <  1:267151/251‑1 (MQ=255)
                                                                                               cgacgcttttccgttctaCGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGgc                                                                                                                                  <  2:257831/234‑1 (MQ=255)
                                                                                                ccccAGAGCGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCAGCATCTCACGCTTTTTCCCCACCATTTCGCCGATGATAACCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGAGGGACTTAAGCAGCAGAAACGgcg                                                                                                                                 >  2:349463/1‑251 (MQ=255)
                                                                                                  ccAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGcgcg                                                                                                                               >  1:474445/1‑251 (MQ=255)
                                                                                                      agagGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGCTAGCCGCAGATTGCAAATTTTCAGCGTCGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGAttt                                                                                                                           >  2:135350/1‑251 (MQ=255)
                                                                                                        agGCCAGACCGCTGCAGGTCCAGGTCTCCTGAGACTGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCAACGCGTGTCCCTTTGGGCAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCACTGATAGCCGCCGATTGCAATTTTCCAACGTAGGTTTCGTTGCGGGAATTAAGAACAAGTACCGGCGGGAttggc                                                                                                                         >  2:280704/1‑248 (MQ=255)
                                                                                                          gccagccagCTGCCTGTCCAGGTCTCCTGAGTTTGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTtgctg                                                                                                                       <  2:287308/243‑1 (MQ=255)
                                                                                                           ccAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCCTCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGGTTtgctgc                                                                                                                      >  2:111608/1‑251 (MQ=255)
                                                                                                             aGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCa                                                                                                                     >  1:414379/1‑250 (MQ=255)
                                                                                                                aCGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGgc                                                                                                                                  >  1:257831/1‑234 (MQ=255)
                                                                                                                     gccgtcccatgTCTCCTGATTATGCAGGTTGAAGTGCGGGTCGCTGGATGTGCCGAAGGCCACGCGTGTCCCTTTTGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGGGCCTGCATCTCGCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAATTTTTCAGCGTAGGTTTCGTTGCGGTACTTAAGCAGCAGTACCGGCGCCATTTTCTGCAATGCGAt                                                                                                              <  2:231668/239‑1 (MQ=255)
                                                                                                                     gCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCCTCCTCTCTTTATGTTGTTCCAGGCGTACCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGCTAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGt                                                                                                            >  2:317155/1‑251 (MQ=255)
                                                                                                                     gCCGGTCCAGGTCCCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGAt                                                                                                              <  2:147581/249‑1 (MQ=255)
                                                                                                                       cGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACGTAAGCAGCAGTTCCGGCGCGATTTGCTGCAAGGCGATGTaa                                                                                                          <  1:369226/251‑1 (MQ=255)
                                                                                                                                tctcCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGc                                                                                                  >  1:421829/1‑250 (MQ=255)
                                                                                                                                         tATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAATGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGActg                                                                                        >  1:397298/1‑251 (MQ=255)
                                                                                                                                           ggCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGGGTCCCTCTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGGTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCCCCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATGTGCTGCAAGGCGATGTAAACCCCCGCATGGCGActgct                                                                                      <  2:379083/250‑1 (MQ=255)
                                                                                                                                            gCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGActgctg                                                                                     >  1:107089/1‑251 (MQ=255)
                                                                                                                                              aGGTTGAATGGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCCCTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTc                                                                                   <  2:76907/251‑1 (MQ=255)
                                                                                                                                               ggTTGAATTGCTGTTCGCGTGATGTGCCAAGGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCg                                                                                  <  2:158941/251‑1 (MQ=255)
                                                                                                                                                    aaTTGCTGTTCGCGTAATGTGCAAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCATCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCCCCATTTCGCCGATGATAGCCGCCGATTGCAACTTTTCAGCGTAGGTTTCGTTGCGGGACTTAAACAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCACGCACGGCGCCTGCTGTCggaca                                                                              >  2:329202/1‑247 (MQ=255)
                                                                                                                                                           gTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAgg                                                                      <  1:472380/251‑1 (MQ=255)
                                                                                                                                                                 gTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAAGCGTGCCAGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTCGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCTATGTAAACCCCCGCATGGAGACTGCTGTCGGCAATGTTCAGGTCTGGt                                                                >  2:174141/1‑251 (MQ=255)
                                                                                                                                                                  tGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGtt                                                               >  1:404468/1‑251 (MQ=255)
                                                                                                                                                                  tGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGACCTCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATCGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCGGCAAGGCGATGTAAAACCCCGCATGGCGACTGCTGTCGGAAATGATCAGGTCTGGtt                                                               >  1:429626/1‑251 (MQ=255)
                                                                                                                                                                       tgCCAAAGGACACGCTTGACCCGTTGTGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCCCCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCTTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCag                                                           <  2:295712/250‑1 (MQ=255)
                                                                                                                                                                          cAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCg                                                       >  1:407285/1‑251 (MQ=255)
                                                                                                                                                                             aGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCa                                                    <  1:76823/251‑1 (MQ=255)
                                                                                                                                                                                    gcgTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTCCCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCtt                                             <  2:116485/251‑1 (MQ=255)
                                                                                                                                                                                     cgTGTCCCTTTGGGACGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTc                                            >  1:403123/1‑251 (MQ=255)
                                                                                                                                                                                        gtCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGCGCGGCAATGGCTTcc                                           >  1:216272/1‑249 (MQ=255)
                                                                                                                                                                                                gggAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCCGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCt                                 >  1:336879/1‑251 (MQ=255)
                                                                                                                                                                                                               ccACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCg                   >  1:239138/1‑250 (MQ=255)
                                                                                                                                                                                                                          cccTCTCTTTATGTTGTGTCAGCGTGGCCTGCAGCTCTCGCTTTTGACCCACCATTGCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTCCCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCGGGCGGCAATGATCAGGTCTGGTTGCAGAGCGGCAATGGCTGCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCAc       <  2:200085/250‑1 (MQ=255)
                                                                                                                                                                                                                                  ttATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGttt  <  1:171731/248‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          
GTGGGCAACCAGCAGCCAGGCAGGATTGACCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTT  >  NC_000913/4513601‑4514074

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