Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,517,355 C→T Q121Q (CAG→CAA)  fecR ← ferric citrate regulator FecR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,517,3550CT98.1% 152.5 / ‑3.7 53Q121Q (CAG→CAA) fecRferric citrate regulator FecR
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (17/35);  total (17/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

TCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTC  >  NC_000913/4517110‑4517597
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
tccccagcggtctgctgcTCGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTCCCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACT‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGAc                                                                                                                                                                                                                                                <  2:325300/239‑1 (MQ=255)
  cacaGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGACCAATCTTCCAGTTGTTGACGGc                                                                                                                                                                                                                                             >  2:125048/1‑251 (MQ=255)
       cGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCCGTTGTTGACGGCTGACt                                                                                                                                                                                                                                        >  2:398387/1‑251 (MQ=255)
       cGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACt                                                                                                                                                                                                                                        >  1:356670/1‑251 (MQ=255)
               tgcACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCttt                                                                                                                                                                                                                                >  1:165388/1‑251 (MQ=255)
                   cGTCAAGCGGCGTGAAATTATCCTGCTGGCGGACGGGAAATTCGTTCCCTAAAGCAGTGAGCTGTCCCGTACGGGTCAGGACACGGAAGGGGCGTGGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGGACCAGAGCCGGACGGTGGGCGGAGGCGCATCAAAACGCACT‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACGGTGCCTTTGgcg                                                                                                                                                                                                                            <  2:352567/251‑1 (MQ=255)
                       aaGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGACGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGCTTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGAATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcggtgc                                                                                                                                                                                                                        >  1:28051/1‑251 (MQ=255)
                       aaGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTT‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcggtgc                                                                                                                                                                                                                        <  1:251594/251‑1 (MQ=255)
                                aaaaTTACCCTGGTGGCGGGCGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCATAGCCGGACGGTGCGCGGATGCGCATCAAAACGAACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTTCGGTAATCTg                                                                                                                                                                                                               <  2:51246/250‑1 (MQ=255)
                                  aaTTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGGCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGTTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGcc                                                                                                                                                                                                             <  1:280288/251‑1 (MQ=255)
                                                 ggACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCg                                                                                                                                                                                              >  2:104559/1‑251 (MQ=255)
                                                     ggTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGttt                                                                                                                                                                                          >  1:398748/1‑251 (MQ=255)
                                                       gaaaTTCTGTCCCTAAAGCAGGGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCGTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCGCGCAGACCTTCGCCGGTTTcc                                                                                                                                                                                        <  1:274910/250‑1 (MQ=255)
                                                              tGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGAGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGACGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTTATCTGCCCTCAGACCTTCTCCGGTTTCCGACttcc                                                                                                                                                                                 >  1:240863/1‑248 (MQ=255)
                                                                   cTAAAGCAGTGAGCTAGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAACACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCCGACCTTCGCCGGGTTCCGACTGCCCGAGc                                                                                                                                                                            >  2:334873/1‑251 (MQ=255)
                                                                      atgCAGTGAGCTGGCCCTGCCGGGTCAGGACACCGTAGGGGCGTGGCAGGGCATCTTCCGCGGGAGTAATGTCGATGTCACCGTACCAGAGCCGGACGTTGCGCTGATGCCCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGTCAGAGCTGc                                                                                                                                                                         <  2:311351/249‑1 (MQ=255)
                                                                      aaGCAGTGAGCTGGCCCGGACGGGTCAGGACACGGAAGGGGGGTTGCAGGGCATCTTGCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTTCGCTGATGCGCATCAAAACGCACC‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCAGGTTCCAGTTGTTGACTGCTGACTGTGCCTTTGGCGGTGCGTTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGc                                                                                                                                                                         <  1:356947/251‑1 (MQ=255)
                                                                       aGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTGCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGTGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGcc                                                                                                                                                                        <  1:30189/251‑1 (MQ=255)
                                                                                 gggCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGcc                                                                                                                                                              <  2:332385/250‑1 (MQ=255)
                                                                                  ggCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGcc                                                                                                                                                              >  1:207574/1‑250 (MQ=255)
                                                                                    cccTGACGGGTCAGGACACGGAAGGGGCGTGGAAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTAGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAgc                                                                                                                                                           <  1:319182/251‑1 (MQ=255)
                                                                                         aCGGGTCAGGACAAGGAAGGGGCGTTGAAGGGCATATTTAGCGGTAGTAAGAGCGATTTCACCGTACCAGAGCGGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGTCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCGTCCGCCAGAGCCgg                                                                                                                                                      >  1:89823/1‑250 (MQ=255)
                                                                                           gggTCAGGACACGGAAGGGGGGTGGCAGGGCACCTTGCGCGGTAGTAATGGCGATTGCCCCGTACCAGAGCCGGACGGTGCGCGGAGGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACGGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGc                                                                                                                                                    <  2:83105/251‑1 (MQ=255)
                                                                                              tCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAAc                                                                                                                                                 >  1:71096/1‑251 (MQ=255)
                                                                                                   gcacGGAAGGGGCGTGGCAGGGCATCTTTCGCGGTAGTAATGGCGCTTGCCCCGTACCAGAGCCGGTCGGTGCGCTGATGCGCATCAAAACCCCCA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAAcagcag                                                                                                                                            <  2:76603/250‑1 (MQ=255)
                                                                                                          gggggCGTTGCAGGGCTGCTTCCGCGTTCGTAAGGGCGATTTCACCGTACCGGAGCCGGACGGTGCGATGGTGCCCATCAAAACTCCCA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTGTGGCGGGGCGGTAATCTGCGCGCAGACCTGCGCCGGTTTCCGACTTCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCt                                                                                                                                     <  2:271641/250‑1 (MQ=255)
                                                                                                           ggggCGTTGCAGGGCATCTTTCGCGGCAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGAGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGCCTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACCGCAGTAATCCtt                                                                                                                                    >  2:192221/1‑251 (MQ=255)
                                                                                                                 ttGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTtcatcc                                                                                                                              >  1:252550/1‑250 (MQ=255)
                                                                                                                  tGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCAc                                                                                                                             >  1:83205/1‑251 (MQ=255)
                                                                                                                  ggCAGGGCACCTGTCGGGGTAGTAATGGCGATTTCCCCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCTCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGTCCTTCGCCGGTTTCCGACTGCCATAGCGGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCAc                                                                                                                             <  2:369587/250‑1 (MQ=255)
                                                                                                                     ggggCATCTTTCGCGGTAGTAATGGCGATTTCCCCGTCCCAGGGCCGGACGTTGCGCTGATGCGCATCAAAACGC‑CACTCCGCCGCTCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTTCCTGTGGCGGTTCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTg                                                                                                                          <  2:51631/250‑1 (MQ=255)
                                                                                                                         cactttccgcgGTAGTAAGGGCGTTTTCCCCTGACCAGAGTCGGACGTGGCGCTGATGCGCATCACAACTCACA‑TCCGCCGCGCTTTGGGGATTCAGCGTG‑AGCAGGGAGTCATCTTCCAGTTGTGGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTGCGCCGGTTTCCGACTGCCAGTGCTGCCATCCTCCGCCATCGCCGAGCAACAGCAGTAATCCTTTCATCGCGTGAcgg                                                                                                                      <  2:277247/244‑1 (MQ=255)
                                                                                                                              ttCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGAcggc                                                                                                                     <  1:398387/247‑1 (MQ=255)
                                                                                                                                         aaTGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTa                                                                                                       <  1:168171/250‑1 (MQ=255)
                                                                                                                                          aTGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCGGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATc                                                                                                     <  1:5251/251‑1 (MQ=255)
                                                                                                                                           tGGCGATTTCACCGTACCAGAGACGGACTGGGCGCTGATGCGCATCAAAGCGCACA‑TCCGCCGAGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTACAGTTGTTGACGGCTGACTGTGCCTTTGGCGGGGCGGTAATCTGCCCGCAGACCTTCGCCGGATTCCGACTGCCAGAGCTGCCAGACTCCGCCAGCTCCGAGCACCAGCACGCATCCTTTCATCACGTggcggcacggtaggtcggggatg                                                                                                    >  2:16972/1‑234 (MQ=255)
                                                                                                                                           gggCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCCCA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCg                                                                                                    <  2:207574/250‑1 (MQ=255)
                                                                                                                                               gATTTCACCGTACCAGGGCCGGACGGTGCGATGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCa                                                                                                <  1:246737/251‑1 (MQ=255)
                                                                                                                                                    ccccGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGc                                                                                            <  1:261956/248‑1 (MQ=255)
                                                                                                                                                      ccGTCCCAGCGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTTACGGCGGGTGAGGCGGGTATTGTGCAACGCCCg                                                                                         <  1:327005/251‑1 (MQ=255)
                                                                                                                                                            ccaggCCGGACGGCGCCCTGATCCGAAGCAAACCGCACA‑CCCGCCGCGCTTTGTGTATTCAGCGTG‑ATCATGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTTCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTCACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGc                                                                                   <  2:231834/247‑1 (MQ=255)
                                                                                                                                                             agagCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCa                                                                                  <  1:256810/251‑1 (MQ=255)
                                                                                                                                                             agagCCGAGCGGTGCGGTGAGGCGCAGCAAACCGCGCA‑TCCGCCGCGCTGTGGGTATTCAGCGGG‑AGCGGGGATCCATCTTCCAGGTGTTGACGGCGGACGTTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCGGCCAGCCTCCGGCAGCGCCGAGCAACAGCAGTAATCCGTGCATCACATGACGGCGGGTGAGGCGGGTATCTTGCAACGCCCGGCTGGCa                                                                                  <  2:158627/251‑1 (MQ=255)
                                                                                                                                                             agagCAGGACGGTGCGCTGATGCGCCTCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGGGAGGCgggaagcgggggaccacacgccgggc                                                                                  >  2:66596/1‑228 (MQ=255)
                                                                                                                                                                  cGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTc                                                                             <  1:364853/251‑1 (MQ=255)
                                                                                                                                                                           cgcTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGcac                                                                    <  1:20883/251‑1 (MQ=255)
                                                                                                                                                                                         aaaccGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTGGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCGGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTAGCGTGCAACGCCCGGCTGGCAACCTCGCCAGGCTCACCGCCAAGCTGGt                                                      <  2:398748/247‑1 (MQ=255)
                                                                                                                                                                                                            gcTTGGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACct                                    <  2:239065/251‑1 (MQ=255)
                                                                                                                                                                                                                tGGGTATTCAGC‑TGTAGCTGTGAGCCAGTTTCCAGTTTTTGACGGCGGACTGTGCCCTTGGCGGGGCGGTAACCTGCCCGCAGACCTTCGCCGGTTGCCGACTGCCAGAGCGGCTAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTTATCACGTGACGGCGGGGGAGTCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACGGCCAAGCTGGTTTCGCAGGTTTTccacaggct                                 <  2:129898/250‑9 (MQ=255)
                                                                                                                                                                                                                 gggTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACctgctg                                <  1:139596/250‑1 (MQ=255)
                                                                                                                                                                                                                 gggTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACctgctg                                <  1:104559/250‑1 (MQ=255)
                                                                                                                                                                                                                                            ccTTCCAGTTGTTGTCGGCTGACTGTGCCGTTGGCGGTGCGGTATGCTGCCCGCCGTCCTTCGCCGGTTTCCGACTGCCAGAGTTGCCAGTCCCGGCCAGCGTCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTTGCAACGTGCCCAGGCACTCCGCCAATCTGTTTGCTCATGTTTTCAACCTGCTGCCTGTCCCTCTGGTTatcctgatcct      <  2:227471/249‑1 (MQ=255)
                                                                                                                                                                                                                                                acAGTTGTTGGCGGCTGTCTGTGCCGTTGGCGGTGCGGTAATCTGCCCCGAGTCCTTCGCCTGTTTCCGTCTGCCCGAGCTGCCAGCCCCCGCCAGCGCAGCGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGTTATCGGGCAACGCCCGGCTGGCAACTTCGCCAGGCTCACCGTTTAGTTGGTTGCGCAGGTTTTCAGCCTGCTGCCAGGCCCACTGGTTATTTTGATCCTGTTc  <  2:388427/249‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
TCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACA‑TCCGCCGCGCTTTGGGTATTCAGCGTG‑AGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTC  >  NC_000913/4517110‑4517597

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