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

AGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGG  >  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)
                                             gTCCTGAGGAATGGCACGGGTAACAACGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGGGGGGCGCGGGTTCCGGCCCCGGGTGCCAAAGTTCATCGCCAGGGCGTGGcggc                                                                                                                                                                                                    >  2:233891/1‑248 (MQ=255)
                                                 ggCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGt                                                                                                                                                                                                <  1:499275/250‑1 (MQ=255)
                                                         ggCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGTGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGGGCCgttgt                                                                                                                                                                                        >  2:411755/1‑251 (MQ=255)
                                                         ggCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCgttgt                                                                                                                                                                                        >  1:421828/1‑251 (MQ=255)
                                                                gTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCgg                                                                                                                                                                                 >  1:516/1‑251 (MQ=255)
                                                                  aaCAAAGTGCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCCCCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCTGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGgc                                                                                                                                                                                <  2:118326/250‑1 (MQ=255)
                                                                          tCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGCTGCCGTCCATCAGGACGGTCGAGCGGCTGGTGCGGCGCGGGTTCAGGCCCCGGATGCCACAGTTCATCGCCCGGTCGGGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTggcg                                                                                                                                                                       >  2:26132/1‑248 (MQ=255)
                                                                          tCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGAc                                                                                                                                                                        <  1:448693/250‑1 (MQ=255)
                                                                            accacGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGcc                                                                                                                                                                     >  1:436333/1‑251 (MQ=255)
                                                                             ccacGCCGCCCACGCTCGGCTGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGTCCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCa                                                                                                                                                                    <  1:68654/251‑1 (MQ=255)
                                                                              cacGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCGTCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGACTGCCAAAGTTCATCGCCCGGTCGCGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGcc                                                                                                                                                                     >  1:433785/1‑249 (MQ=255)
                                                                              cacGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAg                                                                                                                                                                   <  1:472415/251‑1 (MQ=255)
                                                                                  ccgccCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTTCCGTTGTTTTCCGGGGCGCTGACGCCAGGGAt                                                                                                                                                               <  1:154516/251‑1 (MQ=255)
                                                                                  ccacccccgcccggCGGGCCGTACCGCGCCGCACCCGCACCGCGTTCTACGTCAAGGGCCTCCATGTTGCCGCGCGAAACGGGAGCCAGGGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGTACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGAGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGAt                                                                                                                                                               <  2:286115/238‑1 (MQ=255)
                                                                                   cgccCACGCTCGGCGGTCCGTACCGCCCCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATg                                                                                                                                                              <  2:153610/251‑1 (MQ=255)
                                                                                        ccGCTCTGCGGTCCGTAACGCCCCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGTTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGtt                                                                                                                                                         <  1:255650/250‑1 (MQ=255)
                                                                                           ctctGCGGTCCGTAACGCACCGCACCAACACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACAGTAAGGGGCGAAGGGGACGGGGATGCAGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGGTCATCGCCAGGTCGTGGCTGCCGGTGCCGGTGTTTTCCGGTGCGCTGACGCCAGGGATGCGGtga                                                                                                                                                        >  2:66121/1‑247 (MQ=255)
                                                                                               gCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGGGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGt                                                                                                                                                     <  1:180555/248‑1 (MQ=255)
                                                                                                 ggTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATTCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCGGACTCCAGGGATGCGGTTAAGTACCTc                                                                                                                                                <  1:101799/251‑1 (MQ=255)
                                                                                                          cgcaccgcacCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCCtggg                                                                                                                                      >  1:224205/1‑250 (MQ=255)
                                                                                                            caccgcgccACCACCGCGTGCCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGGGAAAGCGTCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGGCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCAtggtgg                                                                                                                                    <  2:432555/251‑1 (MQ=255)
                                                                                                                          cgcggtccACGTCAATGTCATCCATGTTGCCGATCGCAACGGGAGCCAGTGAAAGCTGCGGCTGACCGGAAGGCGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGTCGTTGTTTTTCGGCGTGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGGTGCGCCGGTTTTg                                                                                                                      <  1:472458/245‑1 (MQ=255)
                                                                                                                              tACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGCCCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACCCAGTGGGGTTGCGCCGGTTTGGGCgg                                                                                                                  >  2:326374/1‑250 (MQ=255)
                                                                                                                                ccACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCCGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGGTGCCGAAGTTCATCGCCCGGACGTGGCTGCCGGTGCCGTTGTTTTCCGGCCCGCTGGCGCCCGGGTTGGGGTTACGTCCCTCCCGCCTTGTGGTTGCGCCCGTTTTTGCGaac                                                                                                                >  2:183157/1‑250 (MQ=255)
                                                                                                                                  aCGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATc                                                                                                              <  1:207372/251‑1 (MQ=255)
                                                                                                                                      cAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTc                                                                                                           <  1:233891/250‑1 (MQ=255)
                                                                                                                                      cAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCa                                                                                                          <  1:167928/251‑1 (MQ=255)
                                                                                                                                       aaTGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCGGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGGTGACGCCCGGGGTGCGGGTATGGACCTCACGGATGGTGGGTGCGACGGTTTTGGCGAAATcccccc                                                                                                         >  2:304096/1‑247 (MQ=255)
                                                                                                                                         gggCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCAcgg                                                                                                       <  2:250448/250‑1 (MQ=255)
                                                                                                                                                       ccGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGGAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGCTGCCAAAGTTCATCGCCAGGTAGTGGCTTCCGGTGCCTTTGTTTTCCGGCGCACTTGCGCCAGGGATGCGGTTAAGTACCTCACGCATGGGGGCTGCGCCGGTTTTGGCGAAATCCCCCcggcggtaccacgtccc                                                                                         >  2:269524/1‑240 (MQ=255)
                                                                                                                                                        cGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGGTGCGGTTAAGTACCTCACGCCTGGTGGTTGCGCCGGTTTTGGCGGAATCCTCACGGCGGAttcagtcaccc                                                                                        >  1:293750/1‑241 (MQ=255)
                                                                                                                                                               aaCGGGAGCCAGTGCAAGCTGCGGCGGACCGTAAGGGGCGAAGGGGACGGGGATGGCGGCCATCATGACGGTCGAGCGGCTGGTGAGGCTCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTTCGGCGCGCTGACGTCAGGGATGTGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAg                                                                                 <  1:242603/251‑1 (MQ=255)
                                                                                                                                                                aCGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGc                                                                                >  1:370345/1‑251 (MQ=255)
                                                                                                                                                                 cGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCACAGTTCATCGCCCGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCCCGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGGTCTCGACACGCGGGCCaagc                                                                               >  2:201849/1‑248 (MQ=255)
                                                                                                                                                                       ccAGTGAAAGCAGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGTTGCCGACCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGGCGTGGACGCCGGTGCCGTTGTTTTCCTGCGCGCTGCCGCCAGGGATGCGGTTAAGAGCCTCACGCAAGGGGGTTGCGCCGGTTTTGGCGAACTACTCACAGCGTAtcccgtcacccggcgcaacaggtttc                                                                         >  2:350275/1‑234 (MQ=255)
                                                                                                                                                                         aGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGCTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTGTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCaaa                                                                       <  1:394043/251‑1 (MQ=255)
                                                                                                                                                                            gAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATa                                                                     <  1:441422/250‑1 (MQ=255)
                                                                                                                                                                            gAAAGCTGCGGCTGAACGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGAGGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATAc                                                                    <  2:419272/251‑1 (MQ=255)
                                                                                                                                                                                    cGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGtttt                                                           <  1:474743/251‑1 (MQ=255)
                                                                                                                                                                                          gccGTAAGGGGGGAAGGGGACGGGGATGCCGTCCAGCTGGACGGTCGGGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGTTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCAcgcg                                                     <  2:282232/250‑1 (MQ=255)
                                                                                                                                                                                             gTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAACTACGTCGTTTTCACGCGCAt                                                  >  1:345390/1‑251 (MQ=255)
                                                                                                                                                                                                aGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGTCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCAc                                               <  1:56498/251‑1 (MQ=255)
                                                                                                                                                                                                   ggCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGGTTTCCCGCGCATCAccc                                             >  1:336805/1‑250 (MQ=255)
                                                                                                                                                                                                    gCGAAGGGGACGGGGAGGCCGGCCAGCAGGACGGTCGAGAGGCTGGTGAGGCGCGGGTTCAGGCCCCTGATGCCAAAGTTCATCGCCAGGGCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACccag                                           <  2:445751/251‑1 (MQ=255)
                                                                                                                                                                                                     cGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACccagc                                          >  1:500056/1‑251 (MQ=255)
                                                                                                                                                                                                     cGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGGGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCCGGGATGCGGTTAAGTACCTCCCGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGCCACGCGCGCCCGCCTGTTCCAATAcatctctttccctcgcaccacccccc                                           >  2:2507/1‑226 (MQ=255)
                                                                                                                                                                                                       aaGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATc                                                                                                              >  2:374282/1‑182 (MQ=255)
                                                                                                                                                                                                       aaGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATc                                                                                                              <  1:374282/182‑1 (MQ=255)
                                                                                                                                                                                                           ggACGGGGACGCCGTCCATCTGGTCGGTGGAGGGGGTGGTGAGGCGCCGGTTCAGTCCCCGGAGGCCAAAGTCTATCTCCAGTGTGTGGCGCCCGTTGTTGTTGTTTTCCGGCGCGCTGACGCCAGGTATGCTTTTAATTCCCGCACGCATGGTGGTTGCGCCGGTTATGGCGAAATCCTCACGGCGTATCACGTCACGCGCGCCAGCATGTTCAAATACGTCTTGTTCACGGGCAGTACCCAGCCAGTCg                                    <  2:452389/251‑1 (MQ=255)
                                                                                                                                                                                                                     gCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACgg                          >  1:493132/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ggCGGCTGGTGAGGCGCGGGGTCAGGCCCCGGATGCCGAAGTTCAGCGCCAGGTCGTGGCTTCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCCCGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGt      <  2:317778/250‑1 (MQ=255)
                                                                                                                                                                                                                                            ggCTGGAGAGGCGCGGGATCAGGCCCCGGGTGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGGTTTCCGGCGCGCTGACGCCAGGGCTGCGGTTAAGTACCTCACGCCTGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCCCATCACCCCCCCCAGCATGTTCACATACGTGGGTTTTCCACGCACCCCCCAGCCCACCGCCGCCCACGGTCAGGGCACCCTCTTTTGGGgcg   >  2:93824/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTgcgg  <  1:411755/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        
AGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGG  >  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: