Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,514,139 A→C I57S (ATC→AGC)  fecB ← ferric citrate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,514,1390AC97.9% 140.0 / NA 47I57S (ATC→AGC) fecBferric citrate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (23/23);  minor base G (0/1);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATAT  >  NC_000913/4513900‑4514378
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
aaTTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCAcg                                                                                                                                                                                                                                      <  1:454994/251‑1 (MQ=255)
     tCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACggcggc                                                                                                                                                                                                                                 >  1:306727/1‑251 (MQ=255)
              gtttcgtgtcgtgACTTAAGCAGAAGTACCGGCGCGATTTGCGGCAAGGCGATGGAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTGTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCAGCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGTCCAGCGCATc                                                                                                                                                                                                                        <  2:318100/239‑1 (MQ=255)
                       cGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAAACCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACgg                                                                                                                                                                                                               >  1:316864/1‑250 (MQ=255)
                              tAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAATCCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTGCCAGGCTCGGCTGCGCGCGCGTTCCGACGGTCTGCCACGGTTTCAGGTGCGCACGCTCTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTc                                                                                                                                                                                                         <  1:230435/250‑1 (MQ=255)
                                 gcagcaGTACCGGCGCGATTTGCTGCAAGGCGATGTCAACCCCCGCATGGCGCCTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCACTGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGCCGGACTGCTACGGTTTCAGGTGCGCACGCGCTTCGGGCAGGATGCGTTTTGCCTCGTTATCGTCGGCAATACCGATAGGGCTGACGTCCACGGCGGCCCGCGCATCGGCGCACGagagatacagc                                                                                                                                                                                                     >  2:159888/1‑244 (MQ=255)
                                           cgtcgcgATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATcc                                                                                                                                                                                           <  1:495493/248‑1 (MQ=255)
                                           cGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATcc                                                                                                                                                                                           <  2:104572/251‑1 (MQ=255)
                                            ggCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGCACGAGAGTTCCAGCACAACAAtacg                                                                                                                                                                                          >  1:339581/1‑248 (MQ=255)
                                                       ctgcAAGGCGATGTAAACCCCCGCATAGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCGTCGGCGAACGAGAGTTCCAGCACCACAATCCGTAGTGGCGttt                                                                                                                                                                               >  1:127861/1‑251 (MQ=255)
                                                           aaGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTc                                                                                                                                                                           >  1:399543/1‑251 (MQ=255)
                                                                              cATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCACTGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGCTCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGGGGCGTTTTTTCGCGTGTAAACGTGCCgtgt                                                                                                                                                        >  2:493884/1‑251 (MQ=255)
                                                                                 ggCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTc                                                                                                                                                      <  1:474690/250‑1 (MQ=255)
                                                                                      ctgctgTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGAGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTAGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGCGTTACAGCACCACAATCCGTTGTGGCGTTTTTTCGAGGGGAAACGTGCCGTGTTCGTCCTg                                                                                                                                                >  1:69273/1‑251 (MQ=255)
                                                                                           gTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGGTCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGGGAGTTCCAGCCCCCCCATCCGTTGTGGCGTTTTTTCGAGTGGAAACCTGGCGGGTTCGTCCTGAACCg                                                                                                                                           >  2:32222/1‑251 (MQ=255)
                                                                                             cGGCAATGATCAGGCCTGGTTTCAGAGCGGCAATGGCTTCCAGGCGCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGCTTCAGGTGCGCACGCACTTCGTGCAGGATGCGTTTTTCATCGTTATCGTCGGCAATACCGATCGGGCTGGCGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGGTCCAGCCGCAAAATCCGGTGTGGAGTTTTTTCTAGTTTAAACGTGCCGTGTTCGTCGTGAAcatcg                                                                                                                                         >  2:483519/1‑247 (MQ=255)
                                                                                                       cAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGTGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAgg                                                                                                                               <  1:109709/251‑1 (MQ=255)
                                                                                                        aGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGc                                                                                                                              >  1:407113/1‑251 (MQ=255)
                                                                                                               gTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTg                                                                                                                       <  1:221709/251‑1 (MQ=255)
                                                                                                                       gCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACcag                                                                                                               <  1:32222/251‑1 (MQ=255)
                                                                                                                         ggCAATGGCTTCCGGGCCCGGCCGCGCACGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCCCTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACcagca                                                                                                             <  2:348891/251‑1 (MQ=255)
                                                                                                                            aaTGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACcagcagca                                                                                                          <  2:424662/251‑1 (MQ=255)
                                                                                                                             aggggttCCAGGCCCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACGTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAGCGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTGTCGAGTGTAAACGTGCCGTGGTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGAGCACcagcagcag                                                                                                         <  2:25974/246‑1 (MQ=255)
                                                                                                                                      aGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGGGGCTGATCACCAGCAGCAGGCCTGCaa                                                                                                 >  1:391895/1‑250 (MQ=255)
                                                                                                                                        gCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCaaaaa                                                                                              <  1:121666/251‑1 (MQ=255)
                                                                                                                                          ccGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCaaaaaca                                                                                            <  1:449651/250‑3 (MQ=255)
                                                                                                                                           cGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCCCTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGGTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTTCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCaaaaacaa                                                                                           <  2:266979/251‑4 (MQ=255)
                                                                                                                                            ggCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGACGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGCTCACCAGCAGCAGGCCTGCaaaaacaa                                                                                           >  1:141935/1‑247 (MQ=255)
                                                                                                                                                 cgcgcgcgGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCTCCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGa                                                                                     <  1:266432/251‑1 (MQ=255)
                                                                                                                                                            cGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCCAAGGCGTGGCTGATCACCCGCAGCAGGCCTGCAAAAACAAAACGGATaaaagacaac                                                                          >  2:117482/1‑244 (MQ=255)
                                                                                                                                                                ggACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATa                                                                       >  1:224097/1‑250 (MQ=255)
                                                                                                                                                                  gcggCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATc                                                                    <  2:107494/248‑1 (MQ=255)
                                                                                                                                                                   cTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATc                                                                    >  1:332041/1‑250 (MQ=255)
                                                                                                                                                                     gCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCTCCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATcaca                                                                 <  1:402014/251‑1 (MQ=255)
                                                                                                                                                                     gCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATca                                                                   <  2:332041/249‑1 (MQ=255)
                                                                                                                                                                            tttCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTcccgc                                                          >  1:94142/1‑248 (MQ=255)
                                                                                                                                                                              tCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCt                                                         <  1:157935/250‑1 (MQ=255)
                                                                                                                                                                                                    cAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGAGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGATAAAAGCCCGGCAAGCCGGGCGtt                                  >  1:411315/1‑251 (MQ=255)
                                                                                                                                                                                                           cGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAacaaa                            >  2:195386/1‑248 (MQ=255)
                                                                                                                                                                                                            gTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCCGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCCGCTAAAAGCCCAGCAAGCCGGGCGTTAACCCATc                          >  2:136227/1‑251 (MQ=255)
                                                                                                                                                                                                               ttGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAg                        <  1:6944/250‑1 (MQ=255)
                                                                                                                                                                                                                 gCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAAc                     >  1:330150/1‑251 (MQ=255)
                                                                                                                                                                                                                  caacGTTATCGGCGGCAATACCGAGCGGGGTTACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACACTCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCTCCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACt                    <  2:399543/248‑1 (MQ=255)
                                                                                                                                                                                                                        tATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAAc               <  1:493884/250‑1 (MQ=255)
                                                                                                                                                                                                                            gtcgGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCAGGCAAGCCGGGCGTTAACACATCAGAACTTCAACAAccc          >  2:166524/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCatat  >  2:158673/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aCCGATCGGGCTGACGGCCACGGCGGCCAGCGCATCTGCGAACGAGAGTTCCAGCACCACAATCCGTGGGGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGAGCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCagat  <  2:339581/250‑4 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
AATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATAT  >  NC_000913/4513900‑4514378

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