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→GTGfecAferric 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

GTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCC  >  NC_000913/4514737‑4515226
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
gTTCACCTACGCGGAGCTGGCGTCGATGGAAACGTTATCAAGCGCTGCCGTTAGCGTCCCCAGATCGTAACGGGCCTCCGTTTCCGGCCCGGTACGGCGCGATTTGCAACGTGCAGTGACGGTGTCGTTGGGCTGGTTGGAGTCGTAGTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCTGCGCGCCGTCGTCGTAGCGGGTACCGCGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACAt                                                                                                                                                                                                                                                 >  2:118528/1‑251 (MQ=255)
   cacaTACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTg                                                                                                                                                                                                                                               <  1:91151/250‑1 (MQ=255)
         cGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCttt                                                                                                                                                                                                                                        <  1:129849/251‑1 (MQ=255)
            aTAGCTGGAGTAGATGGAAACGTTATCACGCGTTGGAGGTAGCGTACCCAGCTCGTACCGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCAGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGTGTACACATTGCCGCTTTGc                                                                                                                                                                                                                                      >  2:131517/1‑250 (MQ=255)
             tAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGGGACGGAGACGTTGGTCTGGTTGGGGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTAGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGATTTGcac                                                                                                                                                                                                                                    >  1:267809/1‑251 (MQ=255)
              aGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGcac                                                                                                                                                                                                                                    >  1:262940/1‑250 (MQ=255)
              aGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGcac                                                                                                                                                                                                                                    <  2:79722/250‑1 (MQ=255)
              aGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGcac                                                                                                                                                                                                                                    <  1:105334/250‑1 (MQ=255)
              aGCTGGCGAAGATGGGAACGTTATCAAGCGCTGGCGTTAGCGTACCCAGATCGTACCGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGACACGTGCAGTGACGGTGACGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCTCGCTTTTTCCGGTTCCACCTTGCCGCTTTGcac                                                                                                                                                                                                                                    >  2:234920/1‑250 (MQ=255)
                     gTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGTCCTGTAGGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTTGTCTGGTTGGAGTCGTACTGATTGGTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATGGCCGCTTTGCACAGCCTTGc                                                                                                                                                                                                                            <  2:120434/251‑1 (MQ=255)
                       aGATGGAAACGTTATCCAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGAGTTTACAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCAGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCGACATTGCAGCTTTGCACAGCCTTggc                                                                                                                                                                                                                           >  1:236728/1‑248 (MQ=255)
                       aGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTTCCCATATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGGTAATCAGGAACAGCCCCATGTCCGCCGTCAGCGCGCCGTCGGCGTAGCGGGGACCTAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACAGTTCCGCTTTGCGCAGCCTTGCCa                                                                                                                                                                                                                          <  2:134781/251‑1 (MQ=255)
                             aaaCGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTgg                                                                                                                                                                                                                    >  2:219233/1‑251 (MQ=255)
                                 gTTATCAAGCGTTGGCGTTAGCCTACCCAGATCGTAACGTGCCTGCGATTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGCACAGCCCCATTTCCGACGTCAGCGCGACGTCGTCGTCGCGGGGACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACCTTGCAGCTTTGCCCAGCCTTGcccacatggctgt                                                                                                                                                                                                                >  2:174661/1‑240 (MQ=255)
                                   tATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGctgtac                                                                                                                                                                                                              >  1:129971/1‑251 (MQ=255)
                                      cAAGCGTTGGCGTTAGCGTCCCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGctgtactgt                                                                                                                                                                                                           <  2:3799/251‑1 (MQ=255)
                                       aaGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTGCCACATTGCCCCTTTGCACAGCCTTGCCAATCTGGctgtactgta                                                                                                                                                                                                          <  2:257179/251‑1 (MQ=255)
                                            gtGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTg                                                                                                                                                                                                     >  2:47231/2‑251 (MQ=255)
                                              ggCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGAGTCGTTGGTCTGGTTGGAGCCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTCCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGAACTGTACGGTGcc                                                                                                                                                                                                   >  1:23476/1‑251 (MQ=255)
                                                cGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGa                                                                                                                                                                                                 <  1:219233/251‑1 (MQ=255)
                                                       gTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGcc                                                                                                                                                                                          >  1:202404/1‑251 (MQ=255)
                                                                    aaCGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGc                                                                                                                                                                              >  1:44796/1‑250 (MQ=255)
                                                                              cGTTTCCAGCCCGGTATGGCGCGTTTTTCCACGTGCAGTGACGGTGTCGTTGGTCTGGGTGGAGTCGTACTGATTGTTAAAGTTAAACAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTcc                                                                                                                                                                   <  1:55265/251‑1 (MQ=255)
                                                                              cGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTcc                                                                                                                                                                   <  2:93269/251‑1 (MQ=255)
                                                                               gTTTCCCGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTAGCTGGTCTGGTTGGAGTCGTACTGATTGTTACAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTTCTGTCCGGTGACGAACTTGCCTTACGTGTTTGCATAAAGATTccc                                                                                                                                                                  >  2:32174/1‑250 (MQ=255)
                                                                                        ccGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTgtcagt                                                                                                                                                         >  1:192904/1‑251 (MQ=255)
                                                                                                gcgcGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGAt                                                                                                                                                 <  2:101155/251‑1 (MQ=255)
                                                                                                          ccACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGtt                                                                                                                                       <  1:249807/251‑1 (MQ=255)
                                                                                                             cGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCa                                                                                                                                     >  2:24245/1‑250 (MQ=255)
                                                                                                              gTGCAGGGACGGGGTCGTTGGTCTGGTTGGAGGCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTGCCTGTTCCACATTGCCGCTTTGCACAGCCTTTCCAATCTTGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTGTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAAc                                                                                                                                   <  2:36441/251‑1 (MQ=255)
                                                                                                               tGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACg                                                                                                                                  >  2:88715/1‑251 (MQ=255)
                                                                                                                 cAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACg                                                                                                                                         <  2:71071/242‑1 (MQ=255)
                                                                                                                 cAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACg                                                                                                                                         >  1:71071/1‑242 (MQ=255)
                                                                                                                       cGGTGGCGTTGGTATGGTTGGAGCCGTACTGATTGTTAAAGTTAATCAGGAACAGACCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCACAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCAGAACGAGCCTTCAGTGTTTGCATAAGGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCACCGCCGGAAgg                                                                                                                          >  2:108625/1‑250 (MQ=255)
                                                                                                                        ggTGTCGTTGGTCTGGTTGGAGTCGTCCTTATTGTTAGAGTGAATCAGGAACAGCCCCAGTGCCGCCGTCAGCGCGCCGTCGTCGTAGCGTGTATCGAGTTCCCAGGTTCGCTCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGGCAATCGGGCTGTACTGTACGGTGCCGAACGAGCCTGCAGTGTTTGCATAAAGATTCCAGTTGGCAGACAAGTGATAGAGGACGTTCAACGCCGGAAGCg                                                                                                                         <  2:22927/251‑1 (MQ=255)
                                                                                                                         gtgtCGGTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATATAGCACGTTCAACGCCGGAAGCgg                                                                                                                        <  2:37224/251‑1 (MQ=255)
                                                                                                                             cGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGc                                                                                                                     <  2:212967/250‑1 (MQ=255)
                                                                                                                              gTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCCGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGt                                                                                                                   >  1:41366/1‑251 (MQ=255)
                                                                                                                                         tGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGAGAGAGCACGTTCAACGCCGGAAGCTGTGCGTGATAGCTCAc                                                                                                         <  2:89296/250‑1 (MQ=255)
                                                                                                                                             gTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGa                                                                                                                                                  <  2:87555/205‑1 (MQ=255)
                                                                                                                                               cGTACTGATTGTTAAAGTTAAGCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTcg                                                                                                  <  2:201584/251‑1 (MQ=255)
                                                                                                                                                 tactgcttgctcaaGTTAATCAGGAACAGCCCCATTTCCGCCGCCAGCGCGCCGTCGTCGTAGCGGGTCCCGCGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGCGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCAAGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTcgtg                                                                                                >  2:133383/13‑251 (MQ=255)
                                                                                                                                                     gATTGGTAAAGTTAATCAGGAACAGCCCCATTGCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGAGTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTcgtgcgt                                                                                             <  2:87472/250‑1 (MQ=255)
                                                                                                                                                                 tAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCg                                                                                >  1:86347/1‑251 (MQ=255)
                                                                                                                                                                       ggAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCGTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCCTCGTGCGTGCCTGTGATGGCGTTGTTc                                                                          <  2:16135/251‑1 (MQ=255)
                                                                                                                                                                          aCAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGTCGTTGTTCTgg                                                                       >  2:7954/1‑251 (MQ=255)
                                                                                                                                                                             gCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTAt                                                                    <  1:152304/251‑1 (MQ=255)
                                                                                                                                                                             gCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTAt                                                                    <  1:264382/251‑1 (MQ=255)
                                                                                                                                                                             gCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTAt                                                                    <  1:25818/251‑1 (MQ=255)
                                                                                                                                                                                ccATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGAc                                                                 >  1:170535/1‑251 (MQ=255)
                                                                                                                                                                                      ccgccgTCAGCGCGCCGTCGTCGTAGCGGGTAGCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGGTTGCATAAAGATTCCAGCTGTCAGTCATGTGATAGAGCACGTGCAACGCCGGAAGCGGTGCGTTATAGCTCAGTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGat                                                            <  1:240205/250‑1 (MQ=255)
                                                                                                                                                                                                  cgcCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGGGCGTGCCCGTGATGGCGTTGTTCTGGTTTGACTCGATATGGTCGAAACGc                                               <  2:119784/251‑1 (MQ=255)
                                                                                                                                                                                                      gtcgtcgtAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATAc                                           <  2:95452/251‑1 (MQ=255)
                                                                                                                                                                                                      gtcgtcgtAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCGTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATAc                                           <  1:24245/251‑1 (MQ=255)
                                                                                                                                                                                                           cgtAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACATCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATCAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCgg                                       <  1:143605/250‑1 (MQ=255)
                                                                                                                                                                                                                  ggTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATgg                               >  2:103276/1‑251 (MQ=255)
                                                                                                                                                                                                                      ccGAGTTCCCAGGTTCGCGCTTGTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGTTGCCGAACGAGTCTTCAGTGTTTGCATAAAGATTCCAGCTGTCATTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCTCTTCTTCGTGCTTGCCTGTGATGGCGTTGTCCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGTTCCa                           <  1:131517/251‑1 (MQ=255)
                                                                                                                                                                                                                        gAGTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGAGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGTTATGTTCGAAACGCATACCCGGCGTGATGGTCCAGt                         >  2:60447/1‑251 (MQ=255)
                                                                                                                                                                                                                          gTTCCCAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTg                       <  1:23102/251‑1 (MQ=255)
                                                                                                                                                                                                                               cAGGTTCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGAt                  <  1:214332/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tCGCGCTTTTTCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCg              <  1:168349/250‑1 (MQ=255)
                                                                                                                                                                                                                                       cgcTTTTTCCGGTTCAACATTGCCGACTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCCTCAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGGGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGCTATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGtttt          >  2:40038/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ttttCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGCTTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGt                                    >  2:77929/1‑221 (MQ=255)
                                                                                                                                                                                                                                           ttttCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGt                                    <  1:77929/221‑1 (MQ=255)
                                                                                                                                                                                                                                             ttCCGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGGGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTTATGGCGGTGTTCTGGGATGACTCGATATGTTCGAAACGCATACCCGTCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCa     <  1:241551/250‑1 (MQ=255)
                                                                                                                                                                                                                                               ccGGTTCCACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATcc  >  2:166762/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
GTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCC  >  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: