Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,516,110 C→T A190T (GCC→ACC)  fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,516,1100CT95.0% 160.6 / 0.0 64A190T (GCC→ACC) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (23/34);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01

AGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCG‑TTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGG  >  NC_000913/4515868‑4516355
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        
agagGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGcgcg                                                                                                                                                                                                                                                  >  1:292502/1‑251 (MQ=255)
    gTTGGGCTGAGGTGACCTCCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCGGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGTACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGtt                                                                                                                                                                                                                                              <  2:460175/251‑1 (MQ=255)
         gctggtctgaCCTTCCCCGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGc                                                                                                                                                                                                                                         <  1:181561/245‑1 (MQ=255)
         gctgagctgaCCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGc                                                                                                                                                                                                                                         >  1:317778/1‑251 (MQ=255)
             accgtcctccccACCCCCGCCTTGAGACCAACGGCCTGCGGAATGGCACGGGGAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGGCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCg                                                                                                                                                                                                                                     <  2:2797/241‑1 (MQ=255)
                 gtccTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGGTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGAt                                                                                                                                                                                                                                  <  1:271936/248‑1 (MQ=255)
                    cTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCTGGTTCAGGCCCCGGATGcc                                                                                                                                                                                                                               >  1:118326/1‑250 (MQ=255)
                     tgccccgcccgccTCGATACCAAAGTCCGGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGcc                                                                                                                                                                                                                               <  2:443022/244‑1 (MQ=255)
                       ccACGCCCGCCTCGAGACCAAAGTCCTGCTGAATGGCTCGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAg                                                                                                                                                                                                                           <  1:31879/251‑1 (MQ=255)
                                ccTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGcc                                                                                                                                                                                                                  >  1:428450/1‑251 (MQ=255)
                                        ccAAAGTCCTGCGGAATGGCGCGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCAATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGTCGGGGCTGCCGTCCATAAGGACGGACGAGCGGCTGGGGAGGCGCGGGTTCAGGCCACGGATGCaaaaatcactcgccagggtatg                                                                                                                                                                                                          >  2:417952/1‑233 (MQ=255)
                                             gTCCTGAGGAATGGCACGGGTAACAACGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGGGGGGCGCGGGTTCCGGCCCCGGGTGCCAAAGTTCATCGCCAG‑GGCGTGGcggc                                                                                                                                                                                                    >  2:233891/1‑248 (MQ=255)
                                                 ggCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGt                                                                                                                                                                                                <  1:499275/250‑1 (MQ=255)
                                                         ggCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGTGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGGGCCgttgt                                                                                                                                                                                        >  2:411755/1‑251 (MQ=255)
                                                         ggCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCgttgt                                                                                                                                                                                        >  1:421828/1‑251 (MQ=255)
                                                                gTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCgg                                                                                                                                                                                 >  1:516/1‑251 (MQ=255)
                                                                  aaCAAAGTGCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCCCCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCTGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGgc                                                                                                                                                                                <  2:118326/250‑1 (MQ=255)
                                                                          tCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGCTGCCGTCCATCAGGACGGTCGAGCGGCTGGTGCGGCGCGGGTTCAGGCCCCGGATGCCACAGTTCATCGCCCG‑GTCGGGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTggcg                                                                                                                                                                       >  2:26132/1‑248 (MQ=255)
                                                                          tCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGAc                                                                                                                                                                        <  1:448693/250‑1 (MQ=255)
                                                                            accacGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGcc                                                                                                                                                                     >  1:436333/1‑251 (MQ=255)
                                                                             ccacGCCGCCCACGCTCGGCTGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGTCCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCa                                                                                                                                                                    <  1:68654/251‑1 (MQ=255)
                                                                              cacGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCGTCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGACTGCCAAAGTTCATCGCCCG‑GTCGCGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGcc                                                                                                                                                                     >  1:433785/1‑249 (MQ=255)
                                                                              cacGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAg                                                                                                                                                                   <  1:472415/251‑1 (MQ=255)
                                                                                  ccgccCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTTCCGTTGTTTTCCGGGGCGCTGACGCCAGGGAt                                                                                                                                                               <  1:154516/251‑1 (MQ=255)
                                                                                  ccacccccgcccggCGGGCCGTACCGCGCCGCACCCGCACCGCGTTCTACGTCAAGGGCCTCCATGTTGCCGCGCGAAACGGGAGCCAGGGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGTACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGAGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGAt                                                                                                                                                               <  2:286115/238‑1 (MQ=255)
                                                                                   cgccCACGCTCGGCGGTCCGTACCGCCCCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATg                                                                                                                                                              <  2:153610/251‑1 (MQ=255)
                                                                                        ccGCTCTGCGGTCCGTAACGCCCCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑TTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGtt                                                                                                                                                         <  1:255650/250‑1 (MQ=255)
                                                                                           ctctGCGGTCCGTAACGCACCGCACCAACACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACAGTAAGGGGCGAAGGGGACGGGGATGCAGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGGTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGGTGTTTTCCGGTGCGCTGACGCCAGGGATGCGGtga                                                                                                                                                        >  2:66121/1‑247 (MQ=255)
                                                                                               gCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGGGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGt                                                                                                                                                     <  1:180555/248‑1 (MQ=255)
                                                                                                 ggTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATTCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCGGACTCCAGGGATGCGGTTAAGTACCTc                                                                                                                                                <  1:101799/251‑1 (MQ=255)
                                                                                                          cgcaccgcacCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCC‑tggg                                                                                                                                      >  1:224205/1‑250 (MQ=255)
                                                                                                            caccgcgccACCACCGCGTGCCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGGGAAAGCGTCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGGCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑tggtgg                                                                                                                                    <  2:432555/251‑1 (MQ=255)
                                                                                                                          cgcggtccACGTCAATGTCATCCATGTTGCCGATCGCAACGGGAGCCAGTGAAAGCTGCGGCTGACCGGAAGGCGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGTCGTTGTTTTTCGGCGTGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGGTGCGCCGGTTTTg                                                                                                                      <  1:472458/245‑1 (MQ=255)
                                                                                                                              tACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGCCCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACCCAGTGG‑GGTTGCGCCGGTTTGGGCgg                                                                                                                  >  2:326374/1‑250 (MQ=255)
                                                                                                                                ccACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCCGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGGTGCCGAAGTTCATCGCCCG‑GACGTGGCTGCCGGTGCCGTTGTTTTCCGGCCCGCTGGCGCCCGGGTTGGGGTTACGTCCCTCCCGCC‑TTGTGGTTGCGCCCGTTTTTGCGaac                                                                                                                >  2:183157/1‑250 (MQ=255)
                                                                                                                                  aCGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATc                                                                                                              <  1:207372/251‑1 (MQ=255)
                                                                                                                                      cAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTc                                                                                                           <  1:233891/250‑1 (MQ=255)
                                                                                                                                      cAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCa                                                                                                          <  1:167928/251‑1 (MQ=255)
                                                                                                                                       aaTGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCGGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGGTGACGCCCGGGGTGCGGGTATGGACCTCACGGA‑TGGTGGGTGCGACGGTTTTGGCGAAATcccccc                                                                                                         >  2:304096/1‑247 (MQ=255)
                                                                                                                                         gggCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCAcgg                                                                                                       <  2:250448/250‑1 (MQ=255)
                                                                                                                                                       ccGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGGAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGCTGCCAAAGTTCATCGCCAG‑GTAGTGGCTTCCGGTGCCTTTGTTTTCCGGCGCACTTGCGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGGGGCTGCGCCGGTTTTGGCGAAATCCCCCcggcggtaccacgtccc                                                                                         >  2:269524/1‑240 (MQ=255)
                                                                                                                                                        cGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGGTGCGGTTAAGTACCTCACGCC‑TGGTGGTTGCGCCGGTTTTGGCGGAATCCTCACGGCGGAttcagtcaccc                                                                                        >  1:293750/1‑241 (MQ=255)
                                                                                                                                                               aaCGGGAGCCAGTGCAAGCTGCGGCGGACCGTAAGGGGCGAAGGGGACGGGGATGGCGGCCATCATGACGGTCGAGCGGCTGGTGAGGCTCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTTCGGCGCGCTGACGTCAGGGATGTGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAg                                                                                 <  1:242603/251‑1 (MQ=255)
                                                                                                                                                                aCGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGc                                                                                >  1:370345/1‑251 (MQ=255)
                                                                                                                                                                 cGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCACAGTTCATCGCCCG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCCCGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGGTCTCGACACGCGGGCCaagc                                                                               >  2:201849/1‑248 (MQ=255)
                                                                                                                                                                       ccAGTGAAAGCAGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGTTGCCGACCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GGCGTGGACGCCGGTGCCGTTGTTTTCCTGCGCGCTGCCGCCAGGGATGCGGTTAAGAGCCTCACGCA‑AGGGGGTTGCGCCGGTTTTGGCGAACTACTCACAGCGTAtcccgtcacccggcgcaacaggtttc                                                                         >  2:350275/1‑234 (MQ=255)
                                                                                                                                                                         aGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGCTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTGTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCaaa                                                                       <  1:394043/251‑1 (MQ=255)
                                                                                                                                                                            gAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATa                                                                     <  1:441422/250‑1 (MQ=255)
                                                                                                                                                                            gAAAGCTGCGGCTGAACGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGAGGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATAc                                                                    <  2:419272/251‑1 (MQ=255)
                                                                                                                                                                                    cGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCG‑tttt                                                           <  1:474743/251‑1 (MQ=255)
                                                                                                                                                                                          gccGTAAGGGGGGAAGGGGACGGGGATGCCGTCCAGCTGGACGGTCGGGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGTTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCG‑TTTTCAcgcg                                                     <  2:282232/250‑1 (MQ=255)
                                                                                                                                                                                             gTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAACTACGTCG‑TTTTCACGCGCAt                                                  >  1:345390/1‑251 (MQ=255)
                                                                                                                                                                                                aGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGTCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCG‑TTTTCACGCGCATCAc                                               <  1:56498/251‑1 (MQ=255)
                                                                                                                                                                                                   ggCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCG‑GTTTCCCGCGCATCAccc                                             >  1:336805/1‑250 (MQ=255)
                                                                                                                                                                                                    gCGAAGGGGACGGGGAGGCCGGCCAGCAGGACGGTCGAGAGGCTGGTGAGGCGCGGGTTCAGGCCCCTGATGCCAAAGTTCATCGCCAG‑GGCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCG‑TTTTCACGCGCATCACccag                                           <  2:445751/251‑1 (MQ=255)
                                                                                                                                                                                                     cGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCG‑TTTTCACGCGCATCACccagc                                          >  1:500056/1‑251 (MQ=255)
                                                                                                                                                                                                     cGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGGGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCCGGGATGCGGTTAAGTACCTCCCGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGCCACGCGCGCCCGCCTGTTCCAATAcatctctttccctcgcaccacccccc                                           >  2:2507/1‑226 (MQ=255)
                                                                                                                                                                                                       aaGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATc                                                                                                              >  2:374282/1‑182 (MQ=255)
                                                                                                                                                                                                       aaGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATc                                                                                                              <  1:374282/182‑1 (MQ=255)
                                                                                                                                                                                                           ggACGGGGACGCCGTCCATCTGGTCGGTGGAGGGGGTGGTGAGGCGCCGGTTCAGTCCCCGGAGGCCAAAGTCTATCTCCAGTGT‑GTGGCGCCCGTTGTTGTTGTTTTCCGGCGCGCTGACGCCAGGTATGCTTTTAATTCCCGCACGCA‑TGGTGGTTGCGCCGGTTATGGCGAAATCCTCACGGCGTATCACGTCACGCGCGCCAGCATGTTCAAATACGTCT‑TGTTCACGGGCAGTACCCAGCCAGTCg                                    <  2:452389/251‑1 (MQ=255)
                                                                                                                                                                                                                     gCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCG‑TTTTCACGCGCATCACCCAGCCAGTCGCCGACCACgg                          >  1:493132/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ggCGGCTGGTGAGGCGCGGGGTCAGGCCCCGGATGCCGAAGTTCAGCGCCAG‑GTCGTGGCTTCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCCCGTCACGCGCGCCAGCATGTTCAAATACGTCG‑TTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGt      <  2:317778/250‑1 (MQ=255)
                                                                                                                                                                                                                                            ggCTGGAGAGGCGCGGGATCAGGCCCCGGGTGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGGTTTCCGGCGCGCTGACGCCAGGGCTGCGGTTAAGTACCTCACGCC‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCCCATCACCCCCCCCAGCATGTTCACATACGTGGGTTTTC‑CACGCACCCCCCAGCCCACCGCCGCCCACGGTCAGGGCACCCTCTTTTGGGgcg   >  2:93824/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCG‑TTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTgcgg  <  1:411755/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        
AGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAG‑GTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCA‑TGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCG‑TTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGG  >  NC_000913/4515868‑4516355

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