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,3550CT95.7% 217.0 / NA 71Q121Q (CAG→CAAfecRferric citrate regulator FecR
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (26/40);  minor base G (0/3);  total (28/43)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTG  >  NC_000913/4517110‑4517594
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                
tCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACg                                                                                                                                                                                                                                             >  1:17894/1‑251 (MQ=255)
   acaGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCt                                                                                                                                                                                                                                          >  1:378734/1‑251 (MQ=255)
   acaGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCt                                                                                                                                                                                                                                          <  1:192465/251‑1 (MQ=255)
    caGCGTGCTGCTGCACGACAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTACATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGAGCATCTTTCACGGTAGTAATGGCGATTTCACCGTACCAGTGCCGGACGGTGCGCTGATGCAGATCAAAACGCACATCCGCCACGCCTTGGGTATTCCGCGTGTGCAGGGAGCCCTCTTCCTGTTTCTTACGGCTg                                                                                                                                                                                                                                         >  2:526778/1‑251 (MQ=255)
         ggctgctgaACGTCAAGCGGCGGGAAATGATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTTAGGACACGGAAGGGGCGTGGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACtgt                                                                                                                                                                                                                                    <  1:137618/250‑1 (MQ=255)
            tgctgcACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGcc                                                                                                                                                                                                                                 >  1:295126/1‑251 (MQ=255)
             gctgcACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTAAGCGTGAGAAGGGAGCCATCATCCCGTTGTTGACGGCTGACTTTGCCt                                                                                                                                                                                                                                >  2:381333/1‑251 (MQ=255)
                       aaGCTGCGTGAACTTATCCTGCTGGCGGACGGGAAATTCTGTCACTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGAAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGAAGGTGCGCTGATGCGCATCAAAACGACCATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGATGGCGGATGACAGTGCCTTTGgcggtgg                                                                                                                                                                                                                      >  2:330364/1‑250 (MQ=255)
                       aaGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgtggtgc                                                                                                                                                                                                                      >  1:233803/1‑251 (MQ=255)
                         gCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcggtgcgg                                                                                                                                                                                                                    <  1:13016/251‑1 (MQ=255)
                            gCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTaa                                                                                                                                                                                                                 <  1:56327/251‑1 (MQ=255)
                            gCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCAAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTa                                                                                                                                                                                                                  <  1:518521/250‑1 (MQ=255)
                              gTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATc                                                                                                                                                                                                               <  1:13465/251‑1 (MQ=255)
                               tGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGTGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCt                                                                                                                                                                                                              >  1:203494/1‑251 (MQ=255)
                                  aaTTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGGACCAGAGCCGGCCGGTGCACTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGCATTCAGCGTGGGCAGGGCGCCAACTTCCCGTTGTTGACTGCTGACTGTGCCTTTGGCGGTGGGGTAATATGcc                                                                                                                                                                                                           >  2:385797/1‑251 (MQ=255)
                                   agtacccTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCGGGCCCTGACGGGTCAGGACACGGAGGGGGCGTTGCGGGGCATCTTTCGCGGTAGGAATGGCGATTTCAGCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGccc                                                                                                                                                                                                          <  1:331659/246‑1 (MQ=255)
                                     tATCCTGCGGGCGGACGGTAAATTCTTTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGTTTTCCCCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGTCATCTTCCAGTTGTTGACGGCTTACTGTGCCTTTGGCGGTGCGGTAATCTGCCCg                                                                                                                                                                                                         <  1:296592/250‑1 (MQ=255)
                                      aTCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGcc                                                                                                                                                                                                       >  1:238645/1‑250 (MQ=255)
                                              ggCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTCATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGAATTCAGAGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGGGCGGTAATTTGCCCGCATAccctcc                                                                                                                                                                                               >  2:394891/1‑247 (MQ=255)
                                                 ggACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGCTGCGCATCAAAACGCACATCCGCCGCGCTTTTGGTATTCAGCGTGAGCAGGGCGCCCTCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGGGGGAATCTGCCCGCAGACCTTCGCCg                                                                                                                                                                                            >  2:322696/1‑251 (MQ=255)
                                                     ggaaatgtCTGTCCCTAAAGCAGTGCGCTGGCCCCGACGGGCCAGGACACGGAAGGGTCGTTGCGGGGCATCTTCCGCGGTTGTAATGGCGATGCCCCCGTACCAGAGCCGGACGTTGCGCGAATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTGGGTGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGttt                                                                                                                                                                                        <  2:233803/244‑1 (MQ=255)
                                                         aaTTCTGTCCCTAACGCAGTGAGCTTGACCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCCTCTTTCGCGGTAGTAATGGCGATTTCAACGGACCACCGACGGACGGCGCGCAGGTGCGCATCAAAACGCACATCCGCCGCGCTTTGGGAATTCACCTTGAGCAGGGACGCCTCGTCCCGGTGTTGACGGCTGACTGGGACTTTGGCGGTGCGGTGATCTGGCCTCAGACCATCGCCGGTTTCCgt                                                                                                                                                                                    >  2:270707/1‑250 (MQ=255)
                                                                      aagcggggggcgggCCCTGACGTGTCAGGACACGGAAGGGGTGTTGCAGGGCTGCTGTCGCGGTAGTAATGGCGATTTCACCTTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGGGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGc                                                                                                                                                                       <  2:27443/239‑1 (MQ=255)
                                                                         cAGTGAGCTGGCCCTGCCGGGTCAGGACACGGAAGGGGTGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGTTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAg                                                                                                                                                                    <  1:207567/251‑1 (MQ=255)
                                                                          aGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTATTCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCCGAGCTGCCAGc                                                                                                                                                                   >  1:138261/1‑251 (MQ=255)
                                                                          aGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCCGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGc                                                                                                                                                                   >  1:463609/1‑251 (MQ=255)
                                                                          aGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGc                                                                                                                                                                   >  1:265434/1‑251 (MQ=255)
                                                                             gAGCTGGCCCGGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTc                                                                                                                                                                <  1:528936/251‑1 (MQ=255)
                                                                                cTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGGTTCCGACTGCCAGAGCTGCCAGCCTCCGc                                                                                                                                                             >  1:353110/1‑251 (MQ=255)
                                                                                        gACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGcc                                                                                                                                                      <  1:223920/250‑1 (MQ=255)
                                                                                          cGGGTCAGGACACGGAAGGGGCGTTGCAGGGCACCTTCCGCGGTAGTAAGGGCGATTTCCCCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCTCATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGTTTTCCGACTTCCAGAGCTGCCAGCCTCCGCCAGCGCCGAg                                                                                                                                                   <  1:131202/251‑1 (MQ=255)
                                                                                             gTCAGAACACGGATGGGGCGATGCAGGGCATCTTTCGCGGTAGTAAGGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAACACGCACATCCGCCGCGCTTTGAGTATTCAGCGTGAGCCGGGAGCCATCTTCCTGTAGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCCGACCGTTGACGGTTTCCGGCTGACAGAGCTGCCAGACTccccacgcaccgagca                                                                                                                                                 >  2:106402/1‑236 (MQ=255)
                                                                                                aGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAAcag                                                                                                                                             <  1:322696/251‑1 (MQ=255)
                                                                                                 ggACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCTATTTCACCGTACCAGTGCCGGACGGTGCGCTGATGCGCATCAAACCGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAAcagc                                                                                                                                           <  1:319449/251‑1 (MQ=255)
                                                                                                  gACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAAcagca                                                                                                                                          <  1:291149/251‑1 (MQ=255)
                                                                                                       ggAAGGGGCGTTGCAGTGCACTTTTCGCGATAGCACTGGCGATTTCACCGTACCAGAGCCGGAGGGTGCGCGGAGGCGCATCAAAACGCCCATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACAGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAAt                                                                                                                                     <  2:264402/251‑1 (MQ=255)
                                                                                                         aaGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGAGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGAGCCGAGCAACAGCAGTAATcc                                                                                                                                   >  1:357314/1‑251 (MQ=255)
                                                                                                             ggCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTtc                                                                                                                               <  1:457522/251‑1 (MQ=255)
                                                                                                                gTTGCTGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCCGCGCCGAGCCACAGCAGTAATCCTTtcatc                                                                                                                            >  1:398422/1‑251 (MQ=255)
                                                                                                                gTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCCACAGCAGTAATCCTTtcaac                                                                                                                            >  1:495189/1‑249 (MQ=255)
                                                                                                                  tGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTCCCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTtcatca                                                                                                                           <  1:140673/250‑1 (MQ=255)
                                                                                                                  tGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTtcatca                                                                                                                           <  1:206610/250‑1 (MQ=255)
                                                                                                                       ggCATCTTTCGCGGTAGCAATGGCGATTTCACCGTACCAGAGTCGGACGGTGCTCTGATGCGCATCATACCTCCCCTCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGGCGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCGGAGCTGCCGGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGAc                                                                                                                     <  1:32058/251‑1 (MQ=255)
                                                                                                                         cATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGAcgg                                                                                                                   >  1:430276/1‑251 (MQ=255)
                                                                                                                          aTCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGCCGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTTCGACTGCCcgag                                                                                                                                                                           >  2:67168/1‑195 (MQ=255)
                                                                                                                          aTCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCagag                                                                                                                                                                           <  1:67168/195‑1 (MQ=255)
                                                                                                                            cgttggcgGTAGTAATGGCTATTTCCCCGTCCCGGAGCCGGTCGGTGCGCTGTGTCGCCTAAAACCCAAATTCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACGGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGAcggc                                                                                                                  <  2:521975/244‑1 (MQ=255)
                                                                                                                                cgcgGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGGAATCCTTTCATCACGTGACGGCgggg                                                                                                              >  1:46776/1‑248 (MQ=255)
                                                                                                                                 gcgGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGCGCAACAGCAGTAATCCTTTCATCACGTGACGGCggggg                                                                                                             >  1:98612/1‑247 (MQ=255)
                                                                                                                                 gcgGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAg                                                                                                           <  1:418568/251‑1 (MQ=255)
                                                                                                                                       gtaATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGt                                                                                                     <  1:330364/251‑1 (MQ=255)
                                                                                                                                         aaTGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTAt                                                                                                   <  1:178846/251‑1 (MQ=255)
                                                                                                                                          aTGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCCGCGTGAGCAGGGAGCCATCTGCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGGCCTTCTCCGGGTTCCGCCTGCCAGCGCTGCCCGCCTCCCGCAGCCCCGAACAACCGCAGAAATCCCTTCATCAccgtaccgcgggcgagggcggtgac                                                                                                  >  2:149269/1‑227 (MQ=255)
                                                                                                                                                      ccGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCCGCCCCGAGCCACCACAGTAAACCTTTCATCCCGTGCCGGCGGGGGGGGCgggggtgtagggaatacc                                                                                       >  2:413897/1‑235 (MQ=255)
                                                                                                                                                      ccGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGt                                                                                                     <  1:413897/236‑1 (MQ=255)
                                                                                                                                                          aCCAGAGCCGGCCGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGAATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTATGCCTTTGGCGGTGCGGGAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAACCTCCGCCAACCCCGAGCAACAGACAGAATCATTTCCTCACGTTACGGAGGGAGCGGCggggagtctgagaaccccacgac                                                                                  >  2:500226/1‑231 (MQ=255)
                                                                                                                                                                gCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGCGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACg                                                                            <  1:432323/251‑1 (MQ=255)
                                                                                                                                                                 ccGGACGGTGCGCTGATGCGCATCAAACCGCCCATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGt                                                                           <  1:75723/251‑1 (MQ=255)
                                                                                                                                                                  cGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTc                                                                          <  1:276836/251‑1 (MQ=255)
                                                                                                                                                                          gcgttGATGCGCATCAAACCGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACACCCGGCTGGTAACGTCGCCAGGca                                                                  <  1:295647/247‑1 (MQ=255)
                                                                                                                                                                                  ggcgacataaaccgcccaTCCGCCTCGATTTGGGTATTCAGTGGGCGCAGGGAGCCATCTTCCAGTTGGTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGcc                                                           <  2:311182/235‑1 (MQ=255)
                                                                                                                                                                                  gcgcATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGGTGACGTCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGcc                                                           <  1:267281/250‑1 (MQ=255)
                                                                                                                                                                                                 caaTCCGCCGCGGTTTGGTTATTCAGCCTGAGGAGGGAGCCATCTTACAGTTGGGGTCGGCTGACGTTGCCGTTGGCGGTGCGTGAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCGGCCAGCCTCCGCCAGCGCAGAGCAAAAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGTTGGCAACGTCGCCAGGCACACCGCCAAGCTTGGTGCGCAg                                            <  2:396181/249‑1 (MQ=255)
                                                                                                                                                                                                        gctctgcTTTGGGTATCACGTGTGAGCAGGGTGCCTTTCTCCAGTTGTTGACGGCTGTCTGTGCCCTTGGCGGTGCGGGAATCGTCCCGCAGACCGTCGCCGGTTCCCTACTTCCAGTGCTGCCAGCCTCCGCCAGCTCCGAGCGACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCAGGCTGGCAACGTCGCCAGGTCCACCGCCAAGCTGGTTGCGCAGGtt                                         <  2:517396/243‑1 (MQ=255)
                                                                                                                                                                                                                         cAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGcc                                                                                                                                                                                             <  1:516736/84‑1 (MQ=255)
                                                                                                                                                                                                                         cAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGccc                     <  1:394891/251‑1 (MQ=255)
                                                                                                                                                                                                                              gggggCAGGGAGCCATCTTCCAGTTTTTGACGTCTGACTGTGCCTTTGGGGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTg                 <  2:110495/247‑1 (MQ=255)
                                                                                                                                                                                                                                     gggAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcc          <  1:381333/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gAGCCATCTTCCAGTTGTTGACGGCTGCCTGTGCCTTGGGCGGTGCGGTAATCTGCCCGCCGCCCTTCGCCGGTTCCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctg        <  2:489157/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aGCCATCTTCCAGTTGTGGACGGCTGGCTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctga       <  1:300604/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 tCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctgatcctg  <  1:444225/248‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                
TCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTG  >  NC_000913/4517110‑4517594

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