Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,512,767 C→T A215T (GCG→ACG)  fecC ← ferric citrate ABC transporter membrane subunit FecC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,512,7670CT100.0% 208.5 / NA 60A215T (GCG→ACG) fecCferric citrate ABC transporter membrane subunit FecC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/28);  total (32/28)

CCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCT  >  NC_000913/4512472‑4513059
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                    
ccTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCTCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTcagc                                                                                                                                                                                                                                                                                                  <  2:363002/300‑1 (MQ=255)
            tctcCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCTCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCACCCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGaaa                                                                                                                                                                                                                                                                                      >  2:732479/1‑298 (MQ=255)
            tctcCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCTCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACa                                                                                                                                                                                                                                                                                      >  2:758380/1‑300 (MQ=255)
                 ggggAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCTCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGAc                                                                                                                                                                                                                                                                                 <  2:375460/300‑1 (MQ=255)
                  gggAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCTCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACt                                                                                                                                                                                                                                                                                >  1:509506/1‑300 (MQ=255)
                      aGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCTCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACTTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAg                                                                                                                                                                                                                                                                            <  1:491627/300‑1 (MQ=255)
                         ccAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCTCCAGCAGCANGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCNGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGa                                                                                                                                                                                                                                                                         <  2:785223/300‑1 (MQ=255)
                                      gAGTACATCTGCCAGCAGCATCAGCGTGGCCTCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTAGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCaa                                                                                                                                                                                                                                                            >  2:498826/1‑300 (MQ=255)
                                                    gcagcaTCAGCGTGGCCTCCAGCATCATGCTCACTGGCAGGACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCATAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACAt                                                                                                                                                                                                                                              <  2:715746/300‑1 (MQ=255)
                                                     cagcaTCAGCGTGGCCTCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCATCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATc                                                                                                                                                                                                                                             >  1:496633/1‑300 (MQ=255)
                                                             gCGTGGCCTCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCAACGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCCCTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGCCCACCACCGGCAAGAGCTGCCAGACATCCTGCAAGc                                                                                                                                                                                                                                     >  2:466271/1‑300 (MQ=255)
                                                                  gcctccAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGGTGAACAAACGTAGCCTAGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCAGCAAGAGCTGCAAGACATCCTGCCAGCGGGCg                                                                                                                                                                                                                                >  2:350905/5‑300 (MQ=255)
                                                                  gcctccAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCg                                                                                                                                                                                                                                <  2:422789/296‑1 (MQ=255)
                                                                        agcagcATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGac                                                                                                                                                                                                                          <  1:703412/300‑1 (MQ=255)
                                                                            gcaTGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACANATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTc                                                                                                                                                                                                                      <  2:199823/300‑1 (MQ=255)
                                                                            gcaTGCTCACGGGCAGTACGTTTCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTTACGCACGCACCAACCAGAAGCAGCACTAACANATGGATGACCAAACTTAGCCTCTTCAGGTTCACTCCCAGCGTATGGTCGGTGCTTTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACGTGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTc                                                                                                                                                                                                                      <  2:201046/300‑1 (MQ=255)
                                                                                 cTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGcat                                                                                                                                                                                                                 >  2:819890/1‑298 (MQ=255)
                                                                                       ggCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCAACGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTgcaagc                                                                                                                                                                                                           >  2:448279/1‑300 (MQ=255)
                                                                                             aCGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTaa                                                                                                                                                                                                     >  1:834625/1‑300 (MQ=255)
                                                                                             aCGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCAAGCCAGTaa                                                                                                                                                                                                     >  1:834592/1‑300 (MQ=255)
                                                                                             aCGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCATGCAAGCAAGTaa                                                                                                                                                                                                     >  1:422895/1‑300 (MQ=255)
                                                                                                  gcgcTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTaaaaaat                                                                                                                                                                                                <  1:193449/300‑4 (MQ=255)
                                                                                                  gcgcTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGAt                                                                                                                                                                                                <  1:448518/300‑1 (MQ=255)
                                                                                                       gATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATCGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCTTCATCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGt                                                                                                                                                                                           >  1:773139/1‑300 (MQ=255)
                                                                                                                cTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCCACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCCGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGAt                                                                                                                                                                                  >  1:776020/1‑300 (MQ=255)
                                                                                                                 tGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACNAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATc                                                                                                                                                                                 <  2:755219/300‑1 (MQ=255)
                                                                                                                  gCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGCGGGCCACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCt                                                                                                                                                                                >  2:317253/1‑300 (MQ=255)
                                                                                                                             cgcCAGATGTGGCACCAGCCGACCGATAAACGCCACCGGACCTGCTCCGCTGACGCACGCACCAACCAAAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCCGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGATGAGGCTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCAAGACATCCTGCAAGCGGGCGTGGGACACTCAGCCTGCCAGCAAGTAAAAGATACCGTAAGCATGATCTTCGGCaagcag                                                                                                                                                                     >  2:381452/1‑300 (MQ=255)
                                                                                                                                  gATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCACCACGACAGGGACTGCAGTGACCACCAACGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAggg                                                                                                                                                                >  1:326196/1‑300 (MQ=255)
                                                                                                                                   aTGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAATCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGt                                                                                                                                                               <  2:127340/300‑1 (MQ=255)
                                                                                                                                   aTGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGt                                                                                                                                                               >  2:822676/1‑300 (MQ=255)
                                                                                                                                    tgtgGCACCAGCAGACCGATAAACGCCACCAGACCTGCTACGCTGACGCACGCACCAACCAGAGGCAGCACTAACATATTGAGGACCAAACTTAGCCTCGTCAGGTTCACTCCCAGCGTGTGGTCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCAGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTg                                                                                                                                                              <  2:330293/300‑1 (MQ=255)
                                                                                                                                          aCCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCgg                                                                                                                                                        <  1:471049/300‑1 (MQ=255)
                                                                                                                                                  aCCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCATAAGCAGCACTAACATATTGATGACCAAACCTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTGCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACACCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGcc                                                                                                                                                <  1:353019/300‑1 (MQ=255)
                                                                                                                                                  aCCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGcc                                                                                                                                                <  1:327716/300‑1 (MQ=255)
                                                                                                                                                                      cTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTAAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAGAAGATGCCGTAAGCATGATCTTCGGCAAGCAGCAGGGTGATGCGGGTAAGGACCATACAAAAGGCCGACAgcg                                                                                                                            >  2:432599/1‑300 (MQ=255)
                                                                                                                                                                                 acgcacgcACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCcgcg                                                                                                                 <  1:252600/300‑1 (MQ=255)
                                                                                                                                                                                  cgcacgcaCCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGa                                                                                                                >  2:515814/1‑300 (MQ=255)
                                                                                                                                                                                  cgcacgcaCCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCAACGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCACATACAAAAGGCCGACAGCGCGATACCcgcgc                                                                                                                >  2:761418/1‑299 (MQ=255)
                                                                                                                                                                                                       gcaCTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCAt                                                                                           >  1:645026/1‑300 (MQ=255)
                                                                                                                                                                                                            aaCATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCACCAGCAGGGTGATGCGGGTACGGCACATNCAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTa                                                                                      >  2:620587/1‑300 (MQ=255)
                                                                                                                                                                                                              cATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATg                                                                                    >  1:716385/1‑300 (MQ=255)
                                                                                                                                                                                                                 atTGCTGACCAAACGTAGCCTCATCGGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTGGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACTACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCTTAAGCATGATCTTTGGCCAGCAGCAGGGTGATGCGGTTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACg                                                                                 <  2:291026/300‑1 (MQ=255)
                                                                                                                                                                                                                                         cAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTGGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCATTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTTGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGAc                                                         <  2:251912/300‑1 (MQ=255)
                                                                                                                                                                                                                                               cACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACcagcag                                                   <  1:695429/300‑1 (MQ=255)
                                                                                                                                                                                                                                                       gCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCa                                           >  1:636365/1‑300 (MQ=255)
                                                                                                                                                                                                                                                           aTGGTCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAANACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTTGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTNTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCAcgcc                                       <  2:791094/300‑1 (MQ=255)
                                                                                                                                                                                                                                                           aTGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCATTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCAcgcc                                       <  2:244492/300‑1 (MQ=255)
                                                                                                                                                                                                                                                           aTGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCAcgcc                                       <  1:278025/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGNCAGCCAGTAAAAGATGCCGTAAGCATNATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCg                                >  2:578121/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                         cGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAACAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCACGCATGCcg                         >  2:286246/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                         cGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGAGTAAGGCACATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGCATGACGAAATACGCCTACTACGGTCATGACCAGCAGCCAGCTCACACCGCCCACGCATGCcg                         >  2:538325/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                               ggTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCATCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGatga                   >  1:778282/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                gTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTACTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGatgaa                  >  2:774161/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                    aGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTTGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGAc              <  2:758019/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      cAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACag            <  1:283332/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       aGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACaga           >  1:738521/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                        ggTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACagag          >  2:811893/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                               ttGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCTTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGcc   <  2:488635/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                tGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAATGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCt  <  2:349630/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                    
CCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCT  >  NC_000913/4512472‑4513059

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