Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,516,632 T→C T16A (ACC→GCC)  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,6320TC91.4% 186.5 / 6.6 74T16A (ACC→GCC) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base C (26/38);  total (31/43)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.96e-01

CCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTC  >  NC_000913/4516385‑4516872
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                
cccAGCGGTTTTGCCTGCAGTCCGCTGCCCTCCAGCAGTTGTTGCGGGCCGCTCTCGCCGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGtt                                                                                                                                                                                                                                               <  2:8984/251‑1 (MQ=255)
  cAGCGGTTTTTCCTTCCGGCCGCTCCCGTCCAGCAGTTGTTGCGGTCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCTCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTaa                                                                                                                                                                                                                                             <  2:425704/251‑1 (MQ=255)
                gCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGAGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATGTTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGAGAATGGGGTTAACCAAAGGTgatgtt                                                                                                                                                                                                                               >  2:413925/1‑251 (MQ=255)
                 cAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAACGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGGTAAACAAAGGAGTTGttt                                                                                                                                                                                                                              >  2:34806/1‑251 (MQ=255)
                 cAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGttt                                                                                                                                                                                                                              >  1:231109/1‑251 (MQ=255)
                 cAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGttt                                                                                                                                                                                                                              >  1:115031/1‑251 (MQ=255)
                   gTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTa                                                                                                                                                                                                                            <  1:489930/251‑1 (MQ=255)
                                   cagTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACg                                                                                                                                                                                                            <  1:93824/251‑1 (MQ=255)
                                           gCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGATAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCAt                                                                                                                                                                                                     <  1:217465/250‑1 (MQ=255)
                                            cAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGGGTTGTTTTACGACAAACGCGTAACGGCGTCATAc                                                                                                                                                                                                   >  1:223075/1‑251 (MQ=255)
                                               gCCGCTCCCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCtt                                                                                                                                                                                                <  2:493132/251‑1 (MQ=255)
                                                 cGCTCCCGGCGTCGTAACCGCCGGGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGGCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCGTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCTGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTcc                                                                                                                                                                                              <  2:89232/251‑1 (MQ=255)
                                                       ggACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCcat                                                                                                                                                                                          <  1:284640/248‑1 (MQ=255)
                                                       cggcgtcgtAATCGCCGGGCAGGCCGTTGCTCTGCTGGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCcatca                                                                                                                                                                                        <  2:500056/248‑1 (MQ=255)
                                                       cGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGGGAAGGTAAATCCGCTGTGAGCGGCATACTGATTGAGCACTTTATCGAGCGATCCCGGCGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGAGCGGCCAGCGGCAGCAGGCTCAGGCGACTGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCATAACGGCGTCATACCTTccccctct                                                                                                                                                                                        >  1:313154/1‑247 (MQ=255)
                                                           gtcgtAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCcatcatcat                                                                                                                                                                                    <  1:326374/251‑1 (MQ=255)
                                                           gtcgtAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCcatcatcat                                                                                                                                                                                    >  1:366458/1‑251 (MQ=255)
                                                             cgtAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTAGCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCAGCAttt                                                                                                                                                                                  <  2:463543/251‑1 (MQ=255)
                                                                                gTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGTGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGGGCGCTTTATCGGGCGCTCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGTAAACGCGAGCCCGGCCAGCGGCAGCATGCTCAGGCGATTGGCGTTAACCAAAGGTGTTGTTTTACAACAGACGTGTAACGGCGTCCTACCTTCCCCATCATCCTTTTTGTTGTGTTCAGCTATgc                                                                                                                                                               >  1:20400/1‑250 (MQ=255)
                                                                                   gctttGCTTGCCGCGCGTCAGGCGGGCGTCAACCGAGAGGGTAAATCCGCTGTGGGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCGGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTc                                                                                                                                                            <  2:421828/247‑1 (MQ=255)
                                                                                   gCTCTGCGTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTc                                                                                                                                                            <  1:269524/251‑1 (MQ=255)
                                                                                        gctggCCGGGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGGGGGAGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCAGCATTTTTGTTGTGTTCAGATATGAGTCGAACg                                                                                                                                                       <  2:350379/247‑1 (MQ=255)
                                                                                        gCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAAc                                                                                                                                                        <  1:305804/250‑1 (MQ=255)
                                                                                             ccGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACgaga                                                                                                                                                    >  1:2942/1‑249 (MQ=255)
                                                                                                   gTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTaa                                                                                                                                             >  1:175706/1‑250 (MQ=255)
                                                                                                    gcAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTaaaa                                                                                                                                           <  2:215647/250‑1 (MQ=255)
                                                                                                     cgcgCTGGCGTCAACCTCAGGGGTAAATCCGCTGTGTTGGGCATACTGATTGAGCTCGTTATCAAGCGATCCGCGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGTACGGCCAGCGTCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTCCGAAAATCGCTTAACGGCGTCATACCTTCCCCATCATCATTTTTGGTGTTTTCAGCTATGATTCGAACTAGAAGGGTaaaaa                                                                                                                                          <  2:293750/248‑1 (MQ=255)
                                                                                                        gCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGa                                                                                                                                       <  2:445165/251‑1 (MQ=255)
                                                                                                         cTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGAGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCCGCTATGAGTCGAACGAGCAGGGTAAAAAGGAc                                                                                                                                      >  2:474673/1‑251 (MQ=255)
                                                                                                          tGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGTGTCGAACGAGAAGGGTAAAAAGGAc                                                                                                                                      >  2:376740/1‑250 (MQ=255)
                                                                                                           ggCGTCAACCGAGATGGTAAATCCGCTGTGGGCGGCAGACTGATCGAGCGCTTTATCGAGCGATCTCGGTGCGATACTAACCTGTGCAGCAAAAGCGGAAACCGAGAGACCGGCAAGAGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACCAAGAACGCGTAACGGCGTCATACCTTCCCCATCACCATTGTTGTTGTGTTCAGTTATTAGTCGAACGAGAAGGGTAAAAAGGAc                                                                                                                                      <  1:226530/249‑1 (MQ=255)
                                                                                                           ggCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCTGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGAc                                                                                                                                      <  1:350275/249‑1 (MQ=255)
                                                                                                              gTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGGGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCg                                                                                                                                 <  1:492622/251‑1 (MQ=255)
                                                                                                                 aaCCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCTATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGaaa                                                                                                                              <  1:428805/251‑1 (MQ=255)
                                                                                                                      gagaGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATa                                                                                                                            >  2:377386/1‑249 (MQ=255)
                                                                                                                          gggTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCACAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAatt                                                                                                                      >  2:248580/1‑251 (MQ=255)
                                                                                                                              aaaTCAGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCAAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGatttattt                                                                                                                  >  2:410303/1‑251 (MQ=255)
                                                                                                                                 cccGCGGTGTGCGTCATACTGATTGAGCGCTTTAGCGAGCGGGCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGGAACGGCGGCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATATGAATTATTTTcc                                                                                                               <  2:284703/250‑1 (MQ=255)
                                                                                                                                     cTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAACAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTAc                                                                                                           >  2:152161/1‑251 (MQ=255)
                                                                                                                                              gCATATTGATTGAGCGCTGTAGCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGATTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGaa                                                                                                  <  2:476496/251‑1 (MQ=255)
                                                                                                                                                        tGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATc                                                                                         <  1:377386/250‑1 (MQ=255)
                                                                                                                                                            cgcTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTa                                                                                    <  1:286146/251‑1 (MQ=255)
                                                                                                                                                             gcTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTc                                                                                                                                                            >  1:232648/1‑178 (MQ=255)
                                                                                                                                                             gcTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTc                                                                                                                                                            <  2:232648/178‑1 (MQ=255)
                                                                                                                                                                  aacgagcgggcccGGTGCGATATGAACCAGGGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTGGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACcgcg                                                                               <  2:160159/240‑1 (MQ=255)
                                                                                                                                                                      aGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCGTACAGGGGTGAAATGTTTATCCAGTACCGCGTAATa                                                                          <  1:191040/251‑1 (MQ=255)
                                                                                                                                                                        cGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGa                                                                        <  1:410303/251‑1 (MQ=255)
                                                                                                                                                                         gATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGAc                                                                       >  1:134025/1‑251 (MQ=255)
                                                                                                                                                                            cccGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGATAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGTTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATGTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTg                                                                     <  1:489108/250‑1 (MQ=255)
                                                                                                                                                                            cccGGTGCGATATTAACCATTGCAGCAAAAGTGGAAATCGTAAGACCGGCCAGCGGAAGCAGGCTAAAGCGAATGGAGGTAACCAAAGGTGTTGGTTTACGAAAAACGTGTAACGGCGTCATACCTTCCCAAACATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTAGTATCCTTACAGTGGTGAAATGTTTATCCTGTACCGCGTACTAGACTg                                                                     <  2:433785/250‑1 (MQ=255)
                                                                                                                                                                              cGGTGCGATATTAACCTGTGCAGCAAAAGCGGCAAACGCGAGACCGGCCCGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGATGTGTTCAGCTATGAGTCGAACGCGAAGGGTAAAAAGGACAATCGAAATAAGACTTATTTTCGTTACCGTGGTGAAATGTTTATCCCGTACCGCGTAATAGACTGAAt                                                                  >  2:239895/1‑251 (MQ=255)
                                                                                                                                                                                       atTAACATGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGCAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGAATTAGACTGAATTTTAACggg                                                         >  2:153022/1‑251 (MQ=255)
                                                                                                                                                                                       agtAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCTGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACgg                                                          <  2:395854/248‑1 (MQ=255)
                                                                                                                                                                                            ccTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAAACAAAGGAGTTGTTTTACGAAAAACGCGTAACGGAGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGGGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCg                                                    >  2:95782/1‑251 (MQ=255)
                                                                                                                                                                                              tgtgCAGCAAAAGCGGAACACGAGAGACCGGACAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGGAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGtt                                                  >  2:58599/1‑251 (MQ=255)
                                                                                                                                                                                               gtgCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGtt                                                  <  1:417952/250‑1 (MQ=255)
                                                                                                                                                                                                 gcagcaAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTGTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTg                                                <  1:343145/250‑1 (MQ=255)
                                                                                                                                                                                                     caAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGcgcg                                            >  1:128188/1‑250 (MQ=255)
                                                                                                                                                                                                            ggAAACCGAGAGACCGGCGAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGCCAATCGCATTAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGTCTGACTTTTAACGGGAAGCGTTTGCGCGATAACGtt                                    >  2:385868/1‑251 (MQ=255)
                                                                                                                                                                                                                  cGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGAt                              >  1:61905/1‑251 (MQ=255)
                                                                                                                                                                                                                   gagagaCCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATc                             <  1:389443/251‑1 (MQ=255)
                                                                                                                                                                                                                        aCCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATc                        <  2:231109/251‑1 (MQ=255)
                                                                                                                                                                                                                         ccGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCg                       >  1:353517/1‑251 (MQ=255)
                                                                                                                                                                                                                          cGGCCAGCGGCAGCAGGCTTAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGGGATAACGTTCAGGagcgaagagg                      <  1:16293/251‑10 (MQ=255)
                                                                                                                                                                                                                               aGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTAttt                 <  1:455045/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ggCAGCAGGCGCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGACATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATACCGTTCAGGATCGCATCGGTATTTTTc              >  2:482330/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ccgcagGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACTGCGTCACACCTTCCCTTTCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAg            <  2:109474/249‑1 (MQ=255)
                                                                                                                                                                                                                                         ggcgcAGGCGAATGGCGTTAACCAAAGCTGTTGTGTTACGAAAAACGCGTAACGTCGGCATACCTTCCCCATCATCATTTTTGTTGTGTTCATCTATGATTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGa       <  2:345390/247‑1 (MQ=255)
                                                                                                                                                                                                                                          gCTCAGGCGAATGGCGTTAACCCAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGa       >  1:378132/1‑250 (MQ=255)
                                                                                                                                                                                                                                           cTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTCCGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAg                                                                                                                                                              >  2:503066/1‑98 (MQ=255)
                                                                                                                                                                                                                                           cTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAg                                                                                                                                                              <  1:503066/98‑1 (MQ=255)
                                                                                                                                                                                                                                            tCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACg    >  1:365743/1‑251 (MQ=255)
                                                                                                                                                                                                                                              aGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAATGGCGTCATTCCTTCCCCATCATCATTTTTGTTGTGGTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAAGTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAAGAGACTGAATTTTAACGGGAAGCTTTTGCGCGATAAGGTTCAGGATGGCAACGGGATTTTTCAGCGGGAACGTc  <  2:502283/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTc  >  2:380355/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                
CCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTC  >  NC_000913/4516385‑4516872

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