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,6320TC98.1% 155.6 / ‑5.0 52T16A (ACC→GCC) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (19/32);  total (19/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.32e-01

TACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCC  >  NC_000913/4516396‑4516874
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
tACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGgtg                                                                                                                                                                                                                                       >  1:384023/1‑251 (MQ=255)
  ccTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGtgtt                                                                                                                                                                                                                                     <  1:83541/251‑1 (MQ=255)
     gctgtccgcgtccGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTGGCTCGGCTTGCCGCGCGTCAGGCTGACGTCAACCGAGAGGGTAAATCCGCGGTGTGCGGCATACTGATTGAGCGCTGTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGtgttg                                                                                                                                                                                                                                    <  2:72409/238‑1 (MQ=255)
     gCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGtgttg                                                                                                                                                                                                                                    <  1:252432/249‑1 (MQ=255)
     gCAGTCCGCTGCCGTCCAGCAGTTGGTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGtgttg                                                                                                                                                                                                                                    <  1:46840/249‑1 (MQ=255)
       aGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAAACTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGtttt                                                                                                                                                                                                                                >  1:226291/1‑251 (MQ=255)
               gCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGaaaa                                                                                                                                                                                                                         <  1:124366/250‑1 (MQ=255)
               gCCGTCCAGCAGTTGTTGCAGGCCGCTCCCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGaaaaa                                                                                                                                                                                                                        <  1:363475/251‑1 (MQ=255)
                  gTCCAGCAGTGGTTGCTGGCCGCCTCATACGTCGTAATCTCCGTGCAGTCCGTTGCTCTGCGGGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAAcgc                                                                                                                                                                                                                     <  2:374723/251‑1 (MQ=255)
                            tgtgtcAGGCCGCGCTCTATGTTGTAATCGCCTTGCAGTCTGTCGCTCTGCTTGCCGCGCGTCAGGCGGTCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATCCTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCGTTGCAGCAAAAGCGGAAAACGAGAGACCGGCTAGCGGCAGTAGGCTCAGTCGAATGGCGTTTACCAAAGGTGTTGTTTTACGAAAATCGCGTAACGGCg                                                                                                                                                                                                            <  2:44338/245‑1 (MQ=255)
                                 cAGGCCGCTCTCGAGGTCGTAACCTCCGTGCAGTCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGAGATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCATCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTTCGAAAAACGCGTAACGGCGTCATAc                                                                                                                                                                                                      <  2:129490/251‑1 (MQ=255)
                                   ggCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCt                                                                                                                                                                                                    <  1:89898/251‑1 (MQ=255)
                                        ctctCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGAGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTcccc                                                                                                                                                                                               >  1:225883/1‑251 (MQ=255)
                                          ctcGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTcccc                                                                                                                                                                                               <  1:257375/249‑1 (MQ=255)
                                              gACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCcatcat                                                                                                                                                                                          <  1:333789/251‑1 (MQ=255)
                                                    gtAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAACGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCAtttt                                                                                                                                                                                    >  1:353744/1‑251 (MQ=255)
                                                      aaTCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGTGAGGGTAACTCCGCTGTGTGCGTCATACTGCTTGAGCGCTTTATCGAGCGATCCCGATGAGATATTAACCTGTGCAGCAAAAGCGGAAAACACGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGGTTTGTTTTACGAAAAACGCGTACAGGCGTCATAACTTCCCCATACTCATTTttg                                                                                                                                                                                  >  2:2838/1‑251 (MQ=255)
                                                                    ccGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACAGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCCTCATTTTTGTTGTGTTCAGCTAt                                                                                                                                                                   <  1:12913/251‑1 (MQ=255)
                                                                         gCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTc                                                                                                                                                              <  1:207667/251‑1 (MQ=255)
                                                                                 tGCCGCGCGGAAGGCTGTCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATTGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGTCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACgag                                                                                                                                                       <  1:352703/250‑1 (MQ=255)
                                                                                  gCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACgaga                                                                                                                                                      <  1:65381/250‑1 (MQ=255)
                                                                                   ccttgcgTCAGGCTGGCCTCAACCGTGATGGTAAATCAGCTGTGTGCTTCATTCTGATTGAGCGCTTTAACGGGCGATCCCGGTGCGATACTTACCTTGGCAGCAAAATCGGAATACGATAGACCGGCCAGCGGCAGCATGCTCAGGCGAATTGCGTTAACCAAATGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATGTTGGTTGTGTTCAGCTATGCGTCGAACGAGAAg                                                                                                                                                    <  2:242815/247‑1 (MQ=255)
                                                                                         gTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTaaa                                                                                                                                              <  1:106027/251‑1 (MQ=255)
                                                                                                tgtcGTCAACCGGTAGGGTAAATCCGCTTGGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCTGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCGTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGAc                                                                                                                                        <  2:354997/247‑1 (MQ=255)
                                                                                                 ggCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACaa                                                                                                                                      <  1:287462/251‑1 (MQ=255)
                                                                                                      cAACCGAGAGGGAAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGa                                                                                                                                  >  1:300122/1‑250 (MQ=255)
                                                                                                           gagaGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATa                                                                                                                              >  1:40360/1‑249 (MQ=255)
                                                                                                                      tCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCGGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCCTCATTTTTGTTGTGTTCCGCTATGAGTCGAAAGAGACTGGCAAAAAGGACAATCGAATTAAGAATTATTTTcc                                                                                                                 >  2:107231/1‑251 (MQ=255)
                                                                                                                        cGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCtt                                                                                                               >  1:295297/1‑251 (MQ=255)
                                                                                                                         gCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGCGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCtt                                                                                                               >  1:327124/1‑250 (MQ=255)
                                                                                                                          cTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATCTTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACAAAAGGTGTTGTTTTACGAAAACCGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTa                                                                                                              >  2:107497/1‑250 (MQ=255)
                                                                                                                           tgtgtgCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGGCAATCGAAATAAGAATTATTTTCCTTAc                                                                                                             >  1:59913/1‑250 (MQ=255)
                                                                                                                            gtgtgCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACa                                                                                                            <  1:1628/250‑1 (MQ=255)
                                                                                                                              gtgCGGCATACGGATTGAGCGCTTTACCAAGCTATCCCGGGGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGTCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAgt                                                                                                          <  2:91083/250‑1 (MQ=255)
                                                                                                                                gCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAgtgg                                                                                                        <  1:333562/250‑1 (MQ=255)
                                                                                                                                       aCTGATTGAGCGCTTTAGCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTGCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATg                                                                                                 <  1:374949/250‑1 (MQ=255)
                                                                                                                                               aGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAg                                                                                        <  1:326745/251‑1 (MQ=255)
                                                                                                                                                     ttATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACcg                                                                                   <  1:157015/250‑1 (MQ=255)
                                                                                                                                                                 cccGGTACGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATAACTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGACTTATTTTCCTTATAGTGGTGAAATGTTTATCCAGAACCGCGTAATAAACTg                                                                       >  2:104335/1‑250 (MQ=255)
                                                                                                                                                                  ccGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAg                                                                                                                                                                       <  1:156715/153‑1 (MQ=255)
                                                                                                                                                                  ccGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAg                                                                                                                                                                       >  2:156715/1‑153 (MQ=255)
                                                                                                                                                                       gCGATATTAACCTGGGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAAT‑‑ACTGAATTTTaa                                                               <  2:192834/250‑1 (MQ=255)
                                                                                                                                                                             cgttttacgcttttaaaaGCGGAAAACAAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAAACAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCAACATCATTTTTGTTGTGTTCAGATATGAGTCGAACGAGAAGGGTAACAAGGACAAACGAACTAAGAATTATTATCCTTACACTGGTGAGATGTTTATCCAGTAACTCCTAATAGACTGAATATTacaggg                                                           >  2:363578/15‑245 (MQ=255)
                                                                                                                                                                                                  gAAAACGAGAGACCGACCCGCTGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGCGTTGTTTTCCGCAAAACGCGAAACGGCGTCATACCTTACACATCATCATTTTTGGTGTGTACAGCTATGAGTCGAACGAGACAGGCAAAAAGAAACATAGAAATAAGACTTCTTTTCCTTACAGCGATGAAATGTTTGTCCAGTAAAGCGTAATAGACTGAATTTTAACGGGAAGCGTTTgcgctgtaacattc                                     >  2:175035/1‑241 (MQ=255)
                                                                                                                                                                                                   aaaaCTAGAGACCGGTCAGCGGCAGCAGGTTCAGGCGAATGTCGTGAACCAAAGGTGTTGTTGTACGAAAAACGCGTAACGGCGGCATACCTTCCCCATCAATATTTTTGTTGTGTTCAGCTATGAGGTGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACAGAATTTTAACGGGATGCGTTTGCTCGATAACGTTCa                                    <  2:225883/251‑1 (MQ=255)
                                                                                                                                                                                                        gagagaCCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATc                               <  1:230467/251‑1 (MQ=255)
                                                                                                                                                                                                            gaCCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCAt                           >  1:173227/1‑251 (MQ=255)
                                                                                                                                                                                                              cctgccatcggcagaagtgtgaGGCGCATGGCGTTATCCAAAGGTGTTGTGTTACGAAAAACGCGTAACGGCGTCAGACCATCTCCAACATCATTTTGGTTGTTTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAAGCGAAATAAGAATTATTTTCCTTACAGTGGGTTAATTTTTATCCAGTTACGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCg                         <  2:276060/230‑1 (MQ=255)
                                                                                                                                                                                                                 gCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGt                       >  1:318787/1‑250 (MQ=255)
                                                                                                                                                                                                                     gccgcagcaGGCTCAGGCGAATGGCGTTAACCTAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATTCCCTCCCCATCCTCATTTTTGTTTGGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTGCCTTACAGTGGTGAAATGTTTATCCAGTTCCGCCTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTAtttt                  <  1:27853/248‑1 (MQ=255)
                                                                                                                                                                                                                       gtcagcagtgtCAGGAGAATGGCATTAACCAAATGTGTTGGTTGATGATAAATGCGTTACGGCGTCTTACCCGTCCTATCACCATTTTTGGTTGGTTCTGCTATGAGGCGAACTAGAATGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGGGGTTAAATGTTTATCCAGCACCGCGTAATAGACTGAATTTTCTCGGGAAGCGTTTTGTCTATACCGTTCAGTATCGCAtatgtttttttc                <  2:138655/241‑12 (MQ=255)
                                                                                                                                                                                                                                     gCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTccc  >  1:79073/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
TACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCC  >  NC_000913/4516396‑4516874

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