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.4% 203.9 / NA 66V565V (GTT→GTG) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base C (25/37);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 7.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCC  >  NC_000913/4514737‑4515226
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
gTTCACCTACGCGGAGCTGGCGTCGATGGAAACGTTATCAAGCGCTGCCGTTAGCGTCCCCAGATCGTAACGGGCCTCCGTTTCCGGCCCGGTACGGCGCGATTTGCAACGTGCAGTGACGGTGTCGTTGGGCTGGTTGGAGTCGTAGTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCTGCGCGCCGTCGTCGTAGCGGGTACCGCGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACAt                                                                                                                                                                                                                                                 >  2:118528/1‑251 (MQ=255)
   cacaTACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTg                                                                                                                                                                                                                                               <  1:91151/250‑1 (MQ=255)
         cGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCttt                                                                                                                                                                                                                                        <  1:129849/251‑1 (MQ=255)
            aTAGCTGGAGTAGATGGAAACGTTATCACGCGTTGGAGGTAGCGTACCCAGCTCGTACCGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCAGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGTGTACA‑CATTGCCGCTTTGc                                                                                                                                                                                                                                      >  2:131517/1‑250 (MQ=255)
             tAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGGGACGGAGACGTTGGTCTGGTTGGGGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTAGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGATTTGcac                                                                                                                                                                                                                                    >  1:267809/1‑251 (MQ=255)
              aGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGcac                                                                                                                                                                                                                                    >  1:262940/1‑250 (MQ=255)
              aGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGcac                                                                                                                                                                                                                                    <  2:79722/250‑1 (MQ=255)
              aGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGcac                                                                                                                                                                                                                                    <  1:105334/250‑1 (MQ=255)
              aGCTGGCGAAGATGGGAACGTTATCAAGCGCTGGCGTTAGCGTACCCAGATCGTACCGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGACACGTGCAGTGACGGTGACGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCTCGCTTTTTCCG‑GTTCCACCTTGCCGCTTTGcac                                                                                                                                                                                                                                    >  2:234920/1‑250 (MQ=255)
                     gTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGTCCTGTAGGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTTGTCTGGTTGGAGTCGTACTGATTGGTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATGGCCGCTTTGCACAGCCTTGc                                                                                                                                                                                                                            <  2:120434/251‑1 (MQ=255)
                       aGATGGAAACGTTATCCAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGAGTTTACAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCAGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCGACATTGCAGCTTTGCACAGCCTTggc                                                                                                                                                                                                                           >  1:236728/1‑248 (MQ=255)
                       aGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTTCCCATATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGGTAATCAGGAACAGCCCCATGTCCGCCGTCAGCGCGCCGTCGGCGTAGCGGGGACCTAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACAGTTCCGCTTTGCGCAGCCTTGCCa                                                                                                                                                                                                                          <  2:134781/251‑1 (MQ=255)
                             aaaCGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTgg                                                                                                                                                                                                                    >  2:219233/1‑251 (MQ=255)
                                 gTTATCAAGCGTTGGCGTTAGCCTACCCAGATCGTAACGTGCCTGCGATTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGCACAGCCCCATTTCCGACGTCAGCGCGACGTCGTCGTCGCGGGGACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACCTTGCAGCTTTGCCCAGCCTTGcccacatggctgt                                                                                                                                                                                                                >  2:174661/1‑240 (MQ=255)
                                   tATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGctgtac                                                                                                                                                                                                              >  1:129971/1‑251 (MQ=255)
                                      cAAGCGTTGGCGTTAGCGTCCCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGctgtactgt                                                                                                                                                                                                           <  2:3799/251‑1 (MQ=255)
                                       aaGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTGCCACATTGCCCCTTTGCACAGCCTTGCCAATCTGGctgtactgta                                                                                                                                                                                                          <  2:257179/251‑1 (MQ=255)
                                            gtGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTg                                                                                                                                                                                                     >  2:47231/2‑251 (MQ=255)
                                              ggCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGAGTCGTTGGTCTGGTTGGAGCCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTCCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGAACTGTACGGTGcc                                                                                                                                                                                                   >  1:23476/1‑251 (MQ=255)
                                                cGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGa                                                                                                                                                                                                 <  1:219233/251‑1 (MQ=255)
                                                       gTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGcc                                                                                                                                                                                          >  1:202404/1‑251 (MQ=255)
                                                                    aaCGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGc                                                                                                                                                                              >  1:44796/1‑250 (MQ=255)
                                                                              cGTTTCCAGCCCGGTATGGCGCGTTTTTCCACGTGCAGTGACGGTGTCGTTGGTCTGGGTGGAGTCGTACTGATTGTTAAAGTTAAACAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTcc                                                                                                                                                                   <  1:55265/251‑1 (MQ=255)
                                                                              cGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTcc                                                                                                                                                                   <  2:93269/251‑1 (MQ=255)
                                                                               gTTTCCCGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTAGCTGGTCTGGTTGGAGTCGTACTGATTGTTACAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTTCTGTCCGGTGACGAACTTGCCTTACGTGTTTGCATAAAGATTccc                                                                                                                                                                  >  2:32174/1‑250 (MQ=255)
                                                                                        ccGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTgtcagt                                                                                                                                                         >  1:192904/1‑251 (MQ=255)
                                                                                                gcgcGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGAt                                                                                                                                                 <  2:101155/251‑1 (MQ=255)
                                                                                                          ccACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGtt                                                                                                                                       <  1:249807/251‑1 (MQ=255)
                                                                                                             cGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCa                                                                                                                                     >  2:24245/1‑250 (MQ=255)
                                                                                                              gTGCAGGGACGGGGTCGTTGGTCTGGTTGGAGGCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTGCCT‑GTTCCACATTGCCGCTTTGCACAGCCTTTCCAATCTTGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTGTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAAc                                                                                                                                   <  2:36441/251‑1 (MQ=255)
                                                                                                               tGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACg                                                                                                                                  >  2:88715/1‑251 (MQ=255)
                                                                                                                 cAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACg                                                                                                                                         <  2:71071/242‑1 (MQ=255)
                                                                                                                 cAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACg                                                                                                                                         >  1:71071/1‑242 (MQ=255)
                                                                                                                       cGGTGGCGTTGGTATGGTTGGAGCCGTACTGATTGTTAAAGTTAATCAGGAACAGACCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCACAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCAGAACGAGCCTTCAGTGTTTGCATAAGGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCACCGCCGGAAgg                                                                                                                          >  2:108625/1‑250 (MQ=255)
                                                                                                                        ggTGTCGTTGGTCTGGTTGGAGTCGTCCTTATTGTTAGAGTGAATCAGGAACAGCCCCAGTGCCGCCGTCAGCGCGCCGTCGTCGTAGCGTGTATCGAGTTCCCAGGTTCGCTCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGGCAATCGGGCTGTACTGTACGGTGCCGAACGAGCCTGCAGTGTTTGCATAAAGATTCCAGTTGGCAGACAAGTGATAGAGGACGTTCAACGCCGGAAGCg                                                                                                                         <  2:22927/251‑1 (MQ=255)
                                                                                                                         gtgtCGGTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATATAGCACGTTCAACGCCGGAAGCgg                                                                                                                        <  2:37224/251‑1 (MQ=255)
                                                                                                                             cGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGc                                                                                                                     <  2:212967/250‑1 (MQ=255)
                                                                                                                              gTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCCGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGt                                                                                                                   >  1:41366/1‑251 (MQ=255)
                                                                                                                                         tGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGAGAGAGCACGTTCAACGCCGGAAGCTGTGCGTGATAGCTCAc                                                                                                         <  2:89296/250‑1 (MQ=255)
                                                                                                                                             gTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGa                                                                                                                                                  <  2:87555/205‑1 (MQ=255)
                                                                                                                                               cGTACTGATTGTTAAAGTTAAGCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTcg                                                                                                  <  2:201584/251‑1 (MQ=255)
                                                                                                                                                 tactgcttgctcaaGTTAATCAGGAACAGCCCCATTTCCGCCGCCAGCGCGCCGTCGTCGTAGCGGGTCCCGCGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGCGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCAAGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTcgtg                                                                                                >  2:133383/13‑251 (MQ=255)
                                                                                                                                                     gATTGGTAAAGTTAATCAGGAACAGCCCCATTGCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGAGTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTcgtgcgt                                                                                             <  2:87472/250‑1 (MQ=255)
                                                                                                                                                                 tAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCg                                                                                >  1:86347/1‑251 (MQ=255)
                                                                                                                                                                       ggAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCGTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCCTCGTGCGTGCCTGTGATGGCGTTGTTc                                                                          <  2:16135/251‑1 (MQ=255)
                                                                                                                                                                          aCAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGTCGTTGTTCTgg                                                                       >  2:7954/1‑251 (MQ=255)
                                                                                                                                                                             gCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTAt                                                                    <  1:152304/251‑1 (MQ=255)
                                                                                                                                                                             gCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTAt                                                                    <  1:264382/251‑1 (MQ=255)
                                                                                                                                                                             gCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTAt                                                                    <  1:25818/251‑1 (MQ=255)
                                                                                                                                                                                ccATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGAc                                                                 >  1:170535/1‑251 (MQ=255)
                                                                                                                                                                                      ccgccgTCAGCGCGCCGTCGTCGTAGCGGGTAGCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGGTTGCATAAAGATTCCAGCTGTCAGTCATGTGATAGAGCACGTGCAACGCCGGAAGCGGTGCGTTATAGCTCAGTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGat                                                            <  1:240205/250‑1 (MQ=255)
                                                                                                                                                                                                  cgcCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGGGCGTGCCCGTGATGGCGTTGTTCTGGTTTGACTCGATATGGTCGAAACGc                                               <  2:119784/251‑1 (MQ=255)
                                                                                                                                                                                                      gtcgtcgtAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATAc                                           <  2:95452/251‑1 (MQ=255)
                                                                                                                                                                                                      gtcgtcgtAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCGTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATAc                                           <  1:24245/251‑1 (MQ=255)
                                                                                                                                                                                                           cgtAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACATCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATCAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCgg                                       <  1:143605/250‑1 (MQ=255)
                                                                                                                                                                                                                  ggTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATgg                               >  2:103276/1‑251 (MQ=255)
                                                                                                                                                                                                                      ccGAGTTCCCAGGTTCGCGCTTGTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGTTGCCGAACGAGTCTTCAGTGTTTGCATAAAGATTCCAGCTGTCATTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCTCTTCTTCGTGCTTGCCTGTGATGGCGTTGTCCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGTTCCa                           <  1:131517/251‑1 (MQ=255)
                                                                                                                                                                                                                        gAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGAGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGTTATGTTCGAAACGCATACCCGGCGTGATGGTCCAGt                         >  2:60447/1‑251 (MQ=255)
                                                                                                                                                                                                                          gTTCCCAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTg                       <  1:23102/251‑1 (MQ=255)
                                                                                                                                                                                                                               cAGGTTCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGAt                  <  1:214332/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tCGCGCTTTTTCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCg              <  1:168349/250‑1 (MQ=255)
                                                                                                                                                                                                                                       cgcTTTTTCCG‑GTTCAACATTGCCGACTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCCTCAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGGGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGCTATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGtttt          >  2:40038/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ttttCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGCTTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGt                                    >  2:77929/1‑221 (MQ=255)
                                                                                                                                                                                                                                           ttttCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGt                                    <  1:77929/221‑1 (MQ=255)
                                                                                                                                                                                                                                             ttCCG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGGGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTTATGGCGGTGTTCTGGGATGACTCGATATGTTCGAAACGCATACCCGTCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCa     <  1:241551/250‑1 (MQ=255)
                                                                                                                                                                                                                                               ccG‑GTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATcc  >  2:166762/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
GTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCG‑GTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCC  >  NC_000913/4514737‑4515226

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