Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,514,983 A→C V565V (GTT→GTG fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,514,9830AC95.7% 217.6 / ‑2.5 69V565V (GTT→GTGfecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base C (36/30);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 5.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCAT  >  NC_000913/4514740‑4515224
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
cacaTACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTg                                                                                                                                                                                                                                             <  1:297995/250‑1 (MQ=255)
      cGCATAGCTGGCATAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGCTCGTAACGTGCCTGCGTTTCCAGCCCGGTATAGCGCGTTTTGCCAAGTGCAGTGACGGTGTCGTTGGACTGGTTGGAGTCGTACTGATTGTTAACGTTAATCAGGAACAGCCCCATTTCCGACGTCAGCGCGCAGTCGTCGTACCGGATAACGAGTTCCCAGTTTCGCGCTTTTTCCGGGTCCACATTGCCGAttt                                                                                                                                                                                                                                      >  2:276605/1‑251 (MQ=255)
        cATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGc                                                                                                                                                                                                                                    >  1:36945/1‑251 (MQ=255)
            ggtGGCGGAGATGGAAACGTTTTCAAGCCTTGGCGTTAGCGTCCCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGGGACGGTGTCGTTGGTCTGGTTGGAGTCGTGCTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAg                                                                                                                                                                                                                                <  2:380188/249‑1 (MQ=255)
               ggcttAGATGGAAACGTAATCAATCGTGGGCGTTAGCGTACCCAGATCGGAAAGTGCCTGCGTTTCCAGCCCGGTATGGCGCGGTTTGCCACGGGCAGGGACGGTGTCGTTGGTCGGGTTGGAGTCGTTCTGATTGTTAAAGTGAATCATGACCAGCCCCATTTCCGCCTTCAGCGCGCCGTCTTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGcca                                                                                                                                                                                                                             <  2:167223/247‑2 (MQ=255)
               ggCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCTGGACCCGCGTTCACAGGGTCTCGCTTTTTCCGGTTCCACTTTGCCGCTTTGCACAGCCt                                                                                                                                                                                                                             >  2:211589/1‑251 (MQ=255)
                  gTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCTCATTGCCGCTTTGCACAGCCTTGc                                                                                                                                                                                                                          <  1:317213/251‑1 (MQ=255)
                      aTGGACACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGACGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTCACGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAAt                                                                                                                                                                                                                      >  1:346673/1‑251 (MQ=255)
                      aTGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAAt                                                                                                                                                                                                                      >  1:191175/1‑251 (MQ=255)
                        ggAAACGTTATCAAGCGTTGGCGTTAGAGTACCCAGAGCGTAACGTGCATGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTCACGGTTTCGTTGGCCTGGTTGGGGTCGTACTGATTGTTAAAGATAATCAGGAACGGCCCCCTTTCCGCCGTCAGCGCGCCTTCGGCGTAGCGGGGACCGAGTTCCCCGGTTCGCCCTTTTTCCGGTTCCACGTTGCCGCttttgacagcttttctactcc                                                                                                                                                                                                                    >  2:78510/1‑233 (MQ=255)
                         gAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGCTTGTTAAAGGTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTACCGCTTTGCACCGCCTTGCCAATCTg                                                                                                                                                                                                                   >  1:936/1‑251 (MQ=255)
                             cGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGctg                                                                                                                                                                                                               >  1:297892/1‑251 (MQ=255)
                                  gtcagtcttggttttggtcttcccAGTTCGTAACTTTCCTGCGTTTCCAGCCCGGTATGTCGTGTTTTGCCTCGTGCGGTGACGGGGGCGTGGGGCTGGTTTGAGTCGTCCTGATTGTTCAAGTTATTCAGGAATAGCCCTATTTCCGCCGTCAGCGCGTCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCTGTTCCACTTTGCCGCTTTGCACAGCCTTGTCAATCTGGctgtactg                                                                                                                                                                                                          <  2:211957/230‑1 (MQ=255)
                                   cAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGctgtactgt                                                                                                                                                                                                         <  1:130021/251‑1 (MQ=255)
                                                   cGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGc                                                                                                                                                                                         >  1:81874/1‑251 (MQ=255)
                                                   cGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGc                                                                                                                                                                                         >  1:134644/1‑251 (MQ=255)
                                                         cAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTc                                                                                                                                                                                     >  1:170365/1‑249 (MQ=255)
                                                           gATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGTTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGTTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAgtg                                                                                                                                                                                 <  1:216305/251‑1 (MQ=255)
                                                                  cggtgccttctgttccctcccGGTATGGTGCGTTTTGTGCCGTGCAGTGACGGTCTCGTTGGTCTGGTTGGAGTCGTACTTATAGTTAAAGTTAATCAGGACCAGCCCCATTCCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCTCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCAt                                                                                                                                                                          <  2:211530/233‑1 (MQ=255)
                                                                  aCGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCAt                                                                                                                                                                          <  1:72723/251‑1 (MQ=255)
                                                                    gTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCCATCTGGCTGTACTGTACGGTGCCGGACGAGCCTTACGTGTTTGccta                                                                                                                                                                         >  1:182540/1‑247 (MQ=255)
                                                                              ttCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGc                                                                                                                                                              >  1:219234/1‑251 (MQ=255)
                                                                                 cAGCCCGGTATGGCGCGTTTTGCCTCGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTgt                                                                                                                                                           <  1:96344/251‑1 (MQ=255)
                                                                                   gCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTgtca                                                                                                                                                         <  1:384889/251‑1 (MQ=255)
                                                                                    cccGGGATGGCGCGTTTTGCCACGGGCAGTGACGGGGTCTTTGGTCTGGTTGGAGTCGTACGTAGTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTgtcag                                                                                                                                                        <  2:173000/251‑1 (MQ=255)
                                                                                      cGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGCGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCTCTTTGCACAGCCTTGCCCATCTGGCTGTACTGTACGGTGCCGCACGAGCTTTCCGTGTTTGCATAAAGATTCCAGCTgtcagtc                                                                                                                                                      >  2:91725/1‑251 (MQ=255)
                                                                                       ggTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTgtcagtc                                                                                                                                                      >  1:291424/1‑250 (MQ=255)
                                                                                        gTATGGCGCGTTTTTCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAATTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGGGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTgtcagtc                                                                                                                                                      <  1:107348/249‑1 (MQ=255)
                                                                                                 gTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTc                                                                                                                                                                                     >  1:145254/1‑209 (MQ=255)
                                                                                                 gTTTTGCCACGTGCAGGGACGGTGGCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTc                                                                                                                                                                                     <  2:145254/209‑1 (MQ=255)
                                                                                                    ttgacacggtaaGTGACGGTGTCGGTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAAACAGTAACAGCCCCATGTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCa                                                                                                                                         <  2:115266/239‑1 (MQ=255)
                                                                                                      gccacgggcagtgtagggggcGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACATCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATGGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGt                                                                                                                                      <  2:145275/231‑1 (MQ=255)
                                                                                                           gTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAAc                                                                                                                                 >  1:380553/1‑251 (MQ=255)
                                                                                                                  gACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGaa                                                                                                                          <  1:91725/251‑1 (MQ=255)
                                                                                                                             tGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGtt                                                                                                                <  1:319876/250‑1 (MQ=255)
                                                                                                                             tGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGtt                                                                                                                >  1:253475/1‑250 (MQ=255)
                                                                                                                               gTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTtat                                                                                                              >  1:281369/1‑250 (MQ=255)
                                                                                                                                       ggAGTCGTACTTATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGGCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCAct                                                                                                      <  1:253088/250‑1 (MQ=255)
                                                                                                                                       ggAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGcc                                                                                                                                                                                                 <  1:103615/159‑1 (MQ=255)
                                                                                                                                       ggAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGcc                                                                                                                                                                                                 >  2:103615/1‑159 (MQ=255)
                                                                                                                                          gTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCActtctt                                                                                                  >  1:384187/1‑251 (MQ=255)
                                                                                                                                           tCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTATTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACttcttc                                                                                                 >  1:146471/1‑251 (MQ=255)
                                                                                                                                            ctaactaaTTGATAAAGTTAATCAGGAGCAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGCAACGAGTTCCCCGGTTAGCGCTTTTTCCGGGTCCACATTGCCGCTTTGCACCGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGACCGAGCCTTCCGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTGCAACGCCGGAGGCGTTGCGTGATAGCTCCAttcttct                                                                                                >  2:74235/8‑250 (MQ=255)
                                                                                                                                              tACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTcgtg                                                                                              >  1:195884/1‑251 (MQ=255)
                                                                                                                                                 tGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCGGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGCTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTAAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTcgtgcgg                                                                                           >  1:149883/1‑250 (MQ=255)
                                                                                                                                                  gATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTcgtgcgt                                                                                           <  1:67492/250‑1 (MQ=255)
                                                                                                                                                  gATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTcgtgcg                                                                                            <  1:96858/249‑1 (MQ=255)
                                                                                                                                                         aagtgcaaccatgtaCAGCCCCTTTTCCGTCCACAGCGCGCCGTCGTCGTAGCGGGTCCCGCGTTCCCAGGTTCGCGCTTTGTCCGGTGCCGCATTGCCGCATTGTACAGCCGTGCCAATCTGGCTGTTCTGTACGGTTCCGAACTAGCCTTAAGTGTTTGCATAAATATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGTCGGAAGCGGTTGGTTATAGCTCACGTCTTCGTGCGTGCTTGTGa                                                                                   <  2:121018/237‑1 (MQ=255)
                                                                                                                                                         aaaGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCtgtg                                                                                    <  1:305325/250‑1 (MQ=255)
                                                                                                                                                         aaaGTTAACCCGGAACAGCCCCATTGCCGCCGTTGTGGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTGCTATCTTTTTCCGGTGTCACATTGCCGCTTTGCACAGCCTTGCCAAGCTGGCTGTACTGGACGGTGCAGAACGAGCCTTCAGTGTTTGCATAAAGATTCAAGCTGTCAGTCAGGTGATAGAGCTCGTTCAACGCCGGAAGCGGTGCGGTATAGCGCACTTCTTCGTGCGTGCCtgtg                                                                                    <  2:375152/250‑1 (MQ=255)
                                                                                                                                                                 tCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCgttg                                                                           >  1:272795/1‑251 (MQ=255)
                                                                                                                                                                     gaccAGCCCCATTGCCGCCGTCTGCGCGCCGTCGTTGGAGCGGGTTCCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTTCGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCt                                                                       <  2:359244/248‑1 (MQ=255)
                                                                                                                                                                        cAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGt                                                                    >  1:73472/1‑251 (MQ=255)
                                                                                                                                                                           ccccATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATg                                                                 >  1:30234/1‑251 (MQ=255)
                                                                                                                                                                            cccttttCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTCTTTCCGGTTCCACATTGCCGCTTTGCACAGCATTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATCAAGCTTCCAGCTGTCAGTCAGCGGACAGCGCCCGTTCAACCCCGGAAGCGGTGCGTTCTAGATCACTTCTTCGTGCGGGTCTGTGCTGGCGTTGCTCTGGTAtt                                                                 >  2:40404/5‑249 (MQ=255)
                                                                                                                                                                                  tCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATCGCTCCCTTCTTCGTGCGTGCCTGTGCTGACGTTGTTCTGGTATGACTCgtt                                                          >  2:292051/1‑249 (MQ=255)
                                                                                                                                                                                       cgTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGtt                                                     <  1:247354/251‑1 (MQ=255)
                                                                                                                                                                                                  cgtcgtcgtAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCAt                                           >  1:74890/1‑250 (MQ=255)
                                                                                                                                                                                                     cgtcgtAGCGGGTACCGAGTTCCCAGTTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCTGTGCGt                                                                                                                 <  2:58596/177‑1 (MQ=255)
                                                                                                                                                                                                     cgtcgtAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGt                                                                                                                 >  1:58596/1‑177 (MQ=255)
                                                                                                                                                                                                           ggCGGGTACCGTGTTCCCAGGGTCGCGCTTTTTCCGGTTGCACATTGCCGCTTTGCACAGCCTCGCCCTTCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGATTGCACAAAGCTTACAGCTGTCAGTCAGGTGATAGAGCACGTTCCACGACGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCTGGCGTg                                 <  2:113279/250‑1 (MQ=255)
                                                                                                                                                                                                           aGCGGGTACCGAGTTCCCAGGTTCGAGCAGTTTTCGGATCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGACGAACGGTCCTTCAGGGTTTGCATAAAGATTCCAGGTGTCAGTCAGGTGATAGCGCACGTTCAACGCGGGAAGCGGATCGTTAGAGATCACTTCTTCGTACGTGCTTGTGAAGGCGGTGTTCTTGTATGACTCGAGATGtgtgtaaaccaaagagggcgtg                                 >  2:103281/1‑230 (MQ=255)
                                                                                                                                                                                                               ggTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATgg                             <  1:5177/251‑1 (MQ=255)
                                                                                                                                                                                                               ggTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGACGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGGTGATCTGGTATGACTCGCTATGTTCGAAACGCATACCCGGCGTGATgg                             >  1:76751/1‑251 (MQ=255)
                                                                                                                                                                                                                gTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGCGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTAAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGAGGGAGTTGTTATGGTAGGACTCGATATGTTCGAAACGCATACCCGGCGTGCTGGt                            >  2:292938/1‑251 (MQ=255)
                                                                                                                                                                                                                  aCCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTcc                          <  1:55485/251‑1 (MQ=255)
                                                                                                                                                                                                                              ggTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGt              >  1:50599/1‑251 (MQ=255)
                                                                                                                                                                                                                                 tCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCg            <  1:216827/250‑1 (MQ=255)
                                                                                                                                                                                                                                          ttCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCAt  >  1:351664/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
CACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCAT  >  NC_000913/4514740‑4515224

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