Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,517,558 C→A 5.6% E54* (GAA→TAA)  fecR ← anti‑sigma transmembrane signal transducer for ferric citrate transport; periplasmic FecA‑bound ferric citrate sensor and cytoplasmic FecI ECF sigma factor activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,517,5580CA 5.6% 288.2 / 3.2 83E54* (GAA→TAA) fecRanti‑sigma transmembrane signal transducer for ferric citrate transport; periplasmic FecA‑bound ferric citrate sensor and cytoplasmic FecI ECF sigma factor activator
Reads supporting (aligned to +/- strand):  ref base C (38/40);  new base A (3/2);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAA  >  NC_000913/4517410‑4517705
                                                                                                                                                    |                                                                                                                                                   
gTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTc                                                                                                                                                     <  2:178788/149‑1 (MQ=255)
       aCTGCCAGCGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCACCctgc                                                                                                                                              >  6:120487/1‑149 (MQ=255)
        cTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCCCGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCACCctgct                                                                                                                                             >  4:193147/1‑149 (MQ=255)
           ccAGAGCTGTCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCACGCTGGTTGCGCAGGTTGAAAACctgc                                                                                                                                              >  5:52799/1‑145 (MQ=255)
           ccAGAGCTGTCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTGAAAACctgc                                                                                                                                              <  6:52799/145‑1 (MQ=255)
             agagCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTc                                                                                                                                                     <  8:387950/136‑1 (MQ=255)
             agagCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTGTTc                                                                                                                                                     >  7:387950/1‑136 (MQ=255)
                gCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGcc                                                                                                                                     >  4:137631/1‑149 (MQ=255)
                       gCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGt                                                                                                                              <  7:269644/149‑1 (MQ=255)
                         cTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTa                                                                                                                            >  1:206137/1‑149 (MQ=255)
                             gCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCAc                                                                                                                                  >  6:14595/1‑139 (MQ=255)
                             gCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCAc                                                                                                                                  <  5:14595/139‑1 (MQ=255)
                               cAGCTCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctga                                                                                                                      <  3:259185/149‑1 (MQ=255)
                               cAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctga                                                                                                                      <  3:83397/149‑1 (MQ=255)
                               cAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctga                                                                                                                      <  3:137631/149‑1 (MQ=255)
                                 gcgcCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctgatc                                                                                                                    <  8:324807/149‑1 (MQ=255)
                                        gCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCa                                                                                                             <  6:294153/149‑1 (MQ=255)
                                         cAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCa                                                                                                             >  7:211477/1‑148 (MQ=255)
                                                  tAATCCTTTCATCAGGTGACGGCGGTTGAGGCGGTTATCGTTCAACGCCCTGCTGGAAACGTAGCCAGGCACACCGCCAAGCTGTTTGCGAAGGTTTTAAACCTTCTGCCAGGCCCACTGGTTATCCTGATCCTGTTAATACCACTGtt                                                                                                   <  6:326987/149‑1 (MQ=255)
                                                  tAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGtt                                                                                                   <  7:214644/149‑1 (MQ=255)
                                                    aTCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGc                                                                                                 >  5:367881/1‑149 (MQ=255)
                                                       cTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCa                                                                                               <  8:211477/148‑1 (MQ=255)
                                                       cTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAg                                                                                              >  4:285949/1‑149 (MQ=255)
                                                        tttCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAAc                                                                                                                                                  >  5:245030/1‑96 (MQ=255)
                                                        tttCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAAc                                                                                                                                                  <  6:245030/96‑1 (MQ=255)
                                                         ttCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGCTCCTGTTCATACAACTGTTGCCAgcg                                                                                            >  7:137601/1‑149 (MQ=255)
                                                               aCGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTc                                                                                      <  1:213063/149‑1 (MQ=255)
                                                                cGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTcc                                                                                     <  7:134252/149‑1 (MQ=255)
                                                                 gTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCt                                                                                    >  8:339312/1‑149 (MQ=255)
                                                                 gTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCt                                                                                    >  3:356170/1‑149 (MQ=255)
                                                                     cggcggGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCtgttg                                                                                >  7:59750/1‑149 (MQ=255)
                                                                      ggcggGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCtgttgt                                                                               <  2:347211/149‑1 (MQ=255)
                                                                       gcggGTGAGGCGGGTATCGGGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTtgtg                                                                              <  2:107792/149‑1 (MQ=255)
                                                                        cggGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGCTCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTgg                                                                             >  7:126232/1‑149 (MQ=255)
                                                                         ggggGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGGCGCCAGGCACACCGCCTAGCTGGGTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCTCTGTTGCCAGCGCGCTGCCTGTTGTggg                                                                            <  8:351042/149‑1 (MQ=255)
                                                                         gggTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTACCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTggg                                                                            >  8:281469/1‑149 (MQ=255)
                                                                         gggTAAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTggg                                                                            <  2:255654/149‑1 (MQ=255)
                                                                          ggTGAGGCGGGTATCGAGCCACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGGTGCGCAGGGTTTCAACCTGCTGCCAGGCCCACTGGTTAACCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGc                                                                           <  6:199616/149‑1 (MQ=255)
                                                                              aGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGAc                                                                       <  3:193147/149‑1 (MQ=255)
                                                                                   ggTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCt                                                                  >  1:228715/1‑149 (MQ=255)
                                                                                         gTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCtgtt                                                                                 <  3:129677/128‑1 (MQ=255)
                                                                                         gTGCAACGCCCGGATGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCtgtt                                                                                 >  4:129677/1‑128 (MQ=255)
                                                                                          tGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCACCCTGCTGCCAGGCCCACTGGATATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCt                                                           >  1:293116/1‑149 (MQ=255)
                                                                                          tGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTGCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCt                                                           <  5:120487/149‑1 (MQ=255)
                                                                                            cAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTa                                                         <  5:75864/149‑1 (MQ=255)
                                                                                               cGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCa                                                      <  4:356170/149‑1 (MQ=255)
                                                                                                gCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCAc                                                     <  7:339312/149‑1 (MQ=255)
                                                                                                 cccGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGCTCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACg                                                    <  3:165211/149‑1 (MQ=255)
                                                                                                  ccGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACgg                                                   <  2:206137/149‑1 (MQ=255)
                                                                                                   cGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACgg                                                   >  8:227597/1‑148 (MQ=255)
                                                                                                     gCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCa                                                 >  1:240705/1‑148 (MQ=255)
                                                                                                         gCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTtgtg                                                                              <  1:184142/115‑1 (MQ=255)
                                                                                                         gCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTtgtg                                                                              >  2:184142/1‑115 (MQ=255)
                                                                                                         gCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATAc                                              >  6:107196/1‑147 (MQ=255)
                                                                                                           aaCGTCGGCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGt                                          <  6:367881/149‑1 (MQ=255)
                                                                                                           aaCGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGt                                          >  5:138047/1‑149 (MQ=255)
                                                                                                           aaCGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGt                                          <  4:293320/149‑1 (MQ=255)
                                                                                                           aaCGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGt                                          >  7:320994/1‑149 (MQ=255)
                                                                                                            aCGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTg                                         <  3:56514/149‑1 (MQ=255)
                                                                                                            aCGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTg                                         <  3:210221/149‑1 (MQ=255)
                                                                                                             cGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTgg                                        <  8:59750/149‑1 (MQ=255)
                                                                                                               tCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTggg                                       >  7:2802/1‑148 (MQ=255)
                                                                                                                 gCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAg                                    >  2:2586/1‑149 (MQ=255)
                                                                                                                  ccAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATGCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGc                                   <  5:245465/149‑1 (MQ=255)
                                                                                                                   cAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCt                                  <  2:228715/149‑1 (MQ=255)
                                                                                                                       cacaCCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAAc                              >  4:84256/1‑149 (MQ=255)
                                                                                                                        acacCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACg                             >  5:9055/1‑149 (MQ=255)
                                                                                                                         cacCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGc                            <  3:89765/149‑1 (MQ=255)
                                                                                                                          acCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACg                             >  6:67671/1‑147 (MQ=255)
                                                                                                                             gCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAgcg                        >  3:190791/1‑149 (MQ=255)
                                                                                                                              ccAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCa                                                 <  8:24129/123‑1 (MQ=255)
                                                                                                                              ccAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCa                                                 >  7:24129/1‑123 (MQ=255)
                                                                                                                              ccAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAgcgc                       >  7:345717/1‑149 (MQ=255)
                                                                                                                                aaGCTGGTTGCGCATGGTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGCTCCTGTGCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTcgccgc                                                            <  3:25363/110‑1 (MQ=255)
                                                                                                                                    tGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAgcgc                       <  6:373461/143‑1 (MQ=255)
                                                                                                                                    tGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAgcgc                       >  5:373461/1‑143 (MQ=255)
                                                                                                                                       ttGCGCAGGTTTTCCACCTGCCGCCAGGCCCACTGGTTATCCTGCTCCTGTTCCTACCACTGTTGCCACCGCGCTTCCTGTTGCGGGCTGACGCGCGCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCa                           <  4:297238/136‑1 (MQ=255)
                                                                                                                                       ttGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCa                           >  3:297238/1‑136 (MQ=255)
                                                                                                                                       ttGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGAcggcg              >  6:20243/1‑149 (MQ=255)
                                                                                                                                         cagcAGGTTGAAAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGAcggcgg             >  2:378633/3‑148 (MQ=255)
                                                                                                                                          agcAGGTTGAAAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCggg            >  3:76688/2‑148 (MQ=255)
                                                                                                                                           gcAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGaa          <  4:208813/149‑1 (MQ=255)
                                                                                                                                           gcAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGaa          >  5:80515/1‑149 (MQ=255)
                                                                                                                                                ttttCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAgcgc                       >  2:19599/1‑131 (MQ=255)
                                                                                                                                                ttttCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAgcgc                       <  1:19599/131‑1 (MQ=255)
                                                                                                                                                 tttCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGt     >  8:210643/1‑148 (MQ=255)
                                                                                                                                                    cAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTaa  <  3:285949/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAA  >  NC_000913/4517410‑4517705

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