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→CAAfecRferric 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

TCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTC  >  NC_000913/4517110‑4517597
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
tccccagcggtctgctgcTCGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTCCCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACTTCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGAc                                                                                                                                                                                                                                                <  2:325300/239‑1 (MQ=255)
  cacaGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGACCAATCTTCCAGTTGTTGACGGc                                                                                                                                                                                                                                             >  2:125048/1‑251 (MQ=255)
       cGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCCGTTGTTGACGGCTGACt                                                                                                                                                                                                                                        >  2:398387/1‑251 (MQ=255)
       cGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACt                                                                                                                                                                                                                                        >  1:356670/1‑251 (MQ=255)
               tgcACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCttt                                                                                                                                                                                                                                >  1:165388/1‑251 (MQ=255)
                   cGTCAAGCGGCGTGAAATTATCCTGCTGGCGGACGGGAAATTCGTTCCCTAAAGCAGTGAGCTGTCCCGTACGGGTCAGGACACGGAAGGGGCGTGGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGGACCAGAGCCGGACGGTGGGCGGAGGCGCATCAAAACGCACTTCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACGGTGCCTTTGgcg                                                                                                                                                                                                                            <  2:352567/251‑1 (MQ=255)
                       aaGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGACGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGCTTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGAATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcggtgc                                                                                                                                                                                                                        >  1:28051/1‑251 (MQ=255)
                       aaGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTTAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcggtgc                                                                                                                                                                                                                        <  1:251594/251‑1 (MQ=255)
                                aaaaTTACCCTGGTGGCGGGCGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCATAGCCGGACGGTGCGCGGATGCGCATCAAAACGAACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTTCGGTAATCTg                                                                                                                                                                                                               <  2:51246/250‑1 (MQ=255)
                                  aaTTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGGCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGTTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGcc                                                                                                                                                                                                             <  1:280288/251‑1 (MQ=255)
                                                 ggACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCg                                                                                                                                                                                              >  2:104559/1‑251 (MQ=255)
                                                     ggTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGttt                                                                                                                                                                                          >  1:398748/1‑251 (MQ=255)
                                                       gaaaTTCTGTCCCTAAAGCAGGGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCGTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCGCGCAGACCTTCGCCGGTTTcc                                                                                                                                                                                        <  1:274910/250‑1 (MQ=255)
                                                              tGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGAGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGACGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTTATCTGCCCTCAGACCTTCTCCGGTTTCCGACttcc                                                                                                                                                                                 >  1:240863/1‑248 (MQ=255)
                                                                   cTAAAGCAGTGAGCTAGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAACACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCCGACCTTCGCCGGGTTCCGACTGCCCGAGc                                                                                                                                                                            >  2:334873/1‑251 (MQ=255)
                                                                      atgCAGTGAGCTGGCCCTGCCGGGTCAGGACACCGTAGGGGCGTGGCAGGGCATCTTCCGCGGGAGTAATGTCGATGTCACCGTACCAGAGCCGGACGTTGCGCTGATGCCCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGTCAGAGCTGc                                                                                                                                                                         <  2:311351/249‑1 (MQ=255)
                                                                      aaGCAGTGAGCTGGCCCGGACGGGTCAGGACACGGAAGGGGGGTTGCAGGGCATCTTGCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTTCGCTGATGCGCATCAAAACGCACCTCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCAGGTTCCAGTTGTTGACTGCTGACTGTGCCTTTGGCGGTGCGTTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGc                                                                                                                                                                         <  1:356947/251‑1 (MQ=255)
                                                                       aGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTGCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGTGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGcc                                                                                                                                                                        <  1:30189/251‑1 (MQ=255)
                                                                                 gggCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGcc                                                                                                                                                              <  2:332385/250‑1 (MQ=255)
                                                                                  ggCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGcc                                                                                                                                                              >  1:207574/1‑250 (MQ=255)
                                                                                    cccTGACGGGTCAGGACACGGAAGGGGCGTGGAAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTAGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAgc                                                                                                                                                           <  1:319182/251‑1 (MQ=255)
                                                                                         aCGGGTCAGGACAAGGAAGGGGCGTTGAAGGGCATATTTAGCGGTAGTAAGAGCGATTTCACCGTACCAGAGCGGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGTCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCGTCCGCCAGAGCCgg                                                                                                                                                      >  1:89823/1‑250 (MQ=255)
                                                                                           gggTCAGGACACGGAAGGGGGGTGGCAGGGCACCTTGCGCGGTAGTAATGGCGATTGCCCCGTACCAGAGCCGGACGGTGCGCGGAGGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACGGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGc                                                                                                                                                    <  2:83105/251‑1 (MQ=255)
                                                                                              tCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAAc                                                                                                                                                 >  1:71096/1‑251 (MQ=255)
                                                                                                   gcacGGAAGGGGCGTGGCAGGGCATCTTTCGCGGTAGTAATGGCGCTTGCCCCGTACCAGAGCCGGTCGGTGCGCTGATGCGCATCAAAACCCCCATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAAcagcag                                                                                                                                            <  2:76603/250‑1 (MQ=255)
                                                                                                          gggggCGTTGCAGGGCTGCTTCCGCGTTCGTAAGGGCGATTTCACCGTACCGGAGCCGGACGGTGCGATGGTGCCCATCAAAACTCCCATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTGTGGCGGGGCGGTAATCTGCGCGCAGACCTGCGCCGGTTTCCGACTTCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCt                                                                                                                                     <  2:271641/250‑1 (MQ=255)
                                                                                                           ggggCGTTGCAGGGCATCTTTCGCGGCAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGAGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGCCTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACCGCAGTAATCCtt                                                                                                                                    >  2:192221/1‑251 (MQ=255)
                                                                                                                 ttGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTtcatcc                                                                                                                              >  1:252550/1‑250 (MQ=255)
                                                                                                                  tGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCAc                                                                                                                             >  1:83205/1‑251 (MQ=255)
                                                                                                                  ggCAGGGCACCTGTCGGGGTAGTAATGGCGATTTCCCCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCTCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGTCCTTCGCCGGTTTCCGACTGCCATAGCGGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCAc                                                                                                                             <  2:369587/250‑1 (MQ=255)
                                                                                                                     ggggCATCTTTCGCGGTAGTAATGGCGATTTCCCCGTCCCAGGGCCGGACGTTGCGCTGATGCGCATCAAAACGCCACTCCGCCGCTCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTTCCTGTGGCGGTTCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTg                                                                                                                          <  2:51631/250‑1 (MQ=255)
                                                                                                                         cactttccgcgGTAGTAAGGGCGTTTTCCCCTGACCAGAGTCGGACGTGGCGCTGATGCGCATCACAACTCACATCCGCCGCGCTTTGGGGATTCAGCGTGAGCAGGGAGTCATCTTCCAGTTGTGGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTGCGCCGGTTTCCGACTGCCAGTGCTGCCATCCTCCGCCATCGCCGAGCAACAGCAGTAATCCTTTCATCGCGTGAcgg                                                                                                                      <  2:277247/244‑1 (MQ=255)
                                                                                                                              ttCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGAcggc                                                                                                                     <  1:398387/247‑1 (MQ=255)
                                                                                                                                         aaTGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTa                                                                                                       <  1:168171/250‑1 (MQ=255)
                                                                                                                                          aTGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCGGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATc                                                                                                     <  1:5251/251‑1 (MQ=255)
                                                                                                                                           tGGCGATTTCACCGTACCAGAGACGGACTGGGCGCTGATGCGCATCAAAGCGCACATCCGCCGAGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTACAGTTGTTGACGGCTGACTGTGCCTTTGGCGGGGCGGTAATCTGCCCGCAGACCTTCGCCGGATTCCGACTGCCAGAGCTGCCAGACTCCGCCAGCTCCGAGCACCAGCACGCATCCTTTCATCACGTggcggcacggtaggtcggggatg                                                                                                    >  2:16972/1‑234 (MQ=255)
                                                                                                                                           gggCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCCCATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCg                                                                                                    <  2:207574/250‑1 (MQ=255)
                                                                                                                                               gATTTCACCGTACCAGGGCCGGACGGTGCGATGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCa                                                                                                <  1:246737/251‑1 (MQ=255)
                                                                                                                                                    ccccGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGc                                                                                            <  1:261956/248‑1 (MQ=255)
                                                                                                                                                      ccGTCCCAGCGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTTACGGCGGGTGAGGCGGGTATTGTGCAACGCCCg                                                                                         <  1:327005/251‑1 (MQ=255)
                                                                                                                                                            ccaggCCGGACGGCGCCCTGATCCGAAGCAAACCGCACACCCGCCGCGCTTTGTGTATTCAGCGTGATCATGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTTCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTCACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGc                                                                                   <  2:231834/247‑1 (MQ=255)
                                                                                                                                                             agagCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCa                                                                                  <  1:256810/251‑1 (MQ=255)
                                                                                                                                                             agagCCGAGCGGTGCGGTGAGGCGCAGCAAACCGCGCATCCGCCGCGCTGTGGGTATTCAGCGGGAGCGGGGATCCATCTTCCAGGTGTTGACGGCGGACGTTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCGGCCAGCCTCCGGCAGCGCCGAGCAACAGCAGTAATCCGTGCATCACATGACGGCGGGTGAGGCGGGTATCTTGCAACGCCCGGCTGGCa                                                                                  <  2:158627/251‑1 (MQ=255)
                                                                                                                                                             agagCAGGACGGTGCGCTGATGCGCCTCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGGGAGGCgggaagcgggggaccacacgccgggc                                                                                  >  2:66596/1‑228 (MQ=255)
                                                                                                                                                                  cGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTc                                                                             <  1:364853/251‑1 (MQ=255)
                                                                                                                                                                           cgcTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGcac                                                                    <  1:20883/251‑1 (MQ=255)
                                                                                                                                                                                         aaaccGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTGGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCGGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTAGCGTGCAACGCCCGGCTGGCAACCTCGCCAGGCTCACCGCCAAGCTGGt                                                      <  2:398748/247‑1 (MQ=255)
                                                                                                                                                                                                            gcTTGGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACct                                    <  2:239065/251‑1 (MQ=255)
                                                                                                                                                                                                                tGGGTATTCAGCTGTAGCTGTGAGCCAGTTTCCAGTTTTTGACGGCGGACTGTGCCCTTGGCGGGGCGGTAACCTGCCCGCAGACCTTCGCCGGTTGCCGACTGCCAGAGCGGCTAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTTATCACGTGACGGCGGGGGAGTCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACGGCCAAGCTGGTTTCGCAGGTTTTccacaggct                                 <  2:129898/250‑9 (MQ=255)
                                                                                                                                                                                                                 gggTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACctgctg                                <  1:139596/250‑1 (MQ=255)
                                                                                                                                                                                                                 gggTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACctgctg                                <  1:104559/250‑1 (MQ=255)
                                                                                                                                                                                                                                            ccTTCCAGTTGTTGTCGGCTGACTGTGCCGTTGGCGGTGCGGTATGCTGCCCGCCGTCCTTCGCCGGTTTCCGACTGCCAGAGTTGCCAGTCCCGGCCAGCGTCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTTGCAACGTGCCCAGGCACTCCGCCAATCTGTTTGCTCATGTTTTCAACCTGCTGCCTGTCCCTCTGGTTatcctgatcct      <  2:227471/249‑1 (MQ=255)
                                                                                                                                                                                                                                                acAGTTGTTGGCGGCTGTCTGTGCCGTTGGCGGTGCGGTAATCTGCCCCGAGTCCTTCGCCTGTTTCCGTCTGCCCGAGCTGCCAGCCCCCGCCAGCGCAGCGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGTTATCGGGCAACGCCCGGCTGGCAACTTCGCCAGGCTCACCGTTTAGTTGGTTGCGCAGGTTTTCAGCCTGCTGCCAGGCCCACTGGTTATTTTGATCCTGTTc  <  2:388427/249‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
TCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTC  >  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: