Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,517,127 C→A V197V (GTG→GTT fecR ← ferric citrate regulator FecR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,517,1270CA97.8% 126.3 / ‑4.3 47V197V (GTG→GTTfecRferric citrate regulator FecR
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (21/24);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

CGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGC  >  NC_000913/4516895‑4517371
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
gggcggGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCAGCCAGCGGTTTCACGGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGaaa                                                                                                                                                                                                                                    <  2:177868/250‑1 (MQ=255)
             gcgcagcACGCCGTTGCGGTAACGGGTTAGCGTGGCTCTCACCTCCCCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGCGAAGCACTTCCACAGCGAGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTTCTGGc                                                                                                                                                                                                                       >  2:165175/1‑251 (MQ=255)
                    aCGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGTCGGACGGt                                                                                                                                                                                                                >  1:1795/1‑251 (MQ=255)
                                 aaccggTTAGCGTGGCTATCTCCGCACCAATTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTGATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGGTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACGTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCGGGCGGACGGTAAATGCTGTCCCt                                                                                                                                                                                                   <  2:79073/247‑1 (MQ=255)
                                              ggCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCt                                                                                                                                                                                      <  1:373001/251‑1 (MQ=255)
                                                         caccCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACTATGCGTTGTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGAc                                                                                                                                                                             <  1:222958/249‑1 (MQ=255)
                                                              agtggTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTc                                                                                                                                                                        >  1:358494/4‑249 (MQ=255)
                                                                  gTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAACGCAGTGAGCTGGCCCTGACGGGTCAGGac                                                                                                                                                                   >  2:90923/1‑250 (MQ=255)
                                                                  gTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGac                                                                                                                                                                   >  1:245672/1‑250 (MQ=255)
                                                                                  cAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCCGCGGTTTCCCTGCGCCACACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGATTTTGCGCGGGGGCACTGGCGCGACGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCACGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGATGGCCCTGCCTGGGCAGGACCCGGactgggcgatcc                                                                                                                                                   >  2:235982/1‑239 (MQ=255)
                                                                                  cAGGAGGTCCTTCGTCCAGCTTGTTCTCGAGTCCACCAGCGTATTCACTGCGCCAAACTCAGAGGCGCTGAACTTCAGGCTGTCCCCAGCGTTCCCGAAGCGTTTTTTCGCGGGGGGACTGGCGAGAAGCACTTCCACAGCGTGCTGCTTAACGTCAAGCTGCGTGAAATTATCCTGCTGGGGGACGGTAAATTCGGTCCCTAAAGCAGTGATCTGTCCCTGACGGGTCAGGACACGGAATGGGCGTTGCa                                                                                                                                                  <  2:300122/251‑1 (MQ=255)
                                                                                                   ggCTTGTACTCACTTCATCCAGCGGTGTCCCTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTACACCAGCGTGCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCGGCGTGAAATTGTCCTGCTGGCGGACGGTAAATTCTGTCCCTAATGCAGTGCGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCTCGgta                                                                                                                                 <  2:187528/250‑1 (MQ=255)
                                                                                                   agcttggacacacGTCATCCAGCGGTGTCACTGCGCCAAACTCAGTGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGAGGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTGCCACGGCGTGCGGCTGAGCGTCAATCTGCTTGAAATTATCCTGCTGGCGGGGGGTAAATTCTGTCTCTAAAGCAGTGAGCTGGCCCt                                                                                                                                                                                <  2:158523/192‑1 (MQ=255)
                                                                                                   aGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTACAGCAGTGAGCTGGCCCt                                                                                                                                                                                >  1:158523/1‑204 (MQ=255)
                                                                                                      ttGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGgtagt                                                                                                                               >  1:354750/1‑250 (MQ=255)
                                                                                                      gtttacgcccttcagccAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCATCGTTCACGATGCGTTTTTGCGCGGGGGCACGGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGGAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTGTCGCGgtagt                                                                                                                               <  2:354750/235‑1 (MQ=255)
                                                                                                       tGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACCCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGTTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGgtagta                                                                                                                              <  1:43984/250‑1 (MQ=255)
                                                                                                             ctcGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCg                                                                                                                        >  1:56612/1‑250 (MQ=255)
                                                                                                                  cATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCa                                                                                                                  <  1:165175/251‑1 (MQ=255)
                                                                                                                        gctgTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTAc                                                                                                            <  1:24907/248‑1 (MQ=255)
                                                                                                                          ggTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGACGCTTTCAACGTAcc                                                                                                           >  2:274028/1‑250 (MQ=255)
                                                                                                                             ttCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGCCGGTAAAt                                                                                                                                                                                                            >  2:195798/1‑150 (MQ=255)
                                                                                                                             ttCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAAt                                                                                                                                                                                                            <  1:195798/150‑1 (MQ=255)
                                                                                                                                   gcgcCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGAc                                                                                                 <  1:235982/251‑1 (MQ=255)
                                                                                                                                       cAAACTCAGAGGCGCTGAACTGAAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCCGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGGAATGGCGATTTTCCCGGACCAGAGCCGGACGGTg                                                                                             >  2:252770/1‑251 (MQ=255)
                                                                                                                                        aaaCTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTgc                                                                                            >  1:20533/1‑251 (MQ=255)
                                                                                                                                           cTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCt                                                                                         >  1:152092/1‑251 (MQ=255)
                                                                                                                                                       gAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCaa                                                                             <  1:104335/251‑1 (MQ=255)
                                                                                                                                                        aaCTGCAGGCTGTCCCCAGTGTTCACGATGCTTTTTTGCGCGGGGGCACTGGCGAGAATCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTTAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCaaa                                                                            <  1:142701/251‑1 (MQ=255)
                                                                                                                                                         aCTGCGGGCTGTCACCAGCGTTTACGATGCGTTTTTGCTCGGGGGCACTGGCGATAAGCACTTCCACAGCGTGCTGCTGAACTTCAAGCTGCGTGAAATTATCCTGCGGGCGGACTGTAAATGCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCaaaa                                                                           <  2:40360/251‑1 (MQ=255)
                                                                                                                                                              aGGCTCTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGcac                                                                      >  1:349798/1‑251 (MQ=255)
                                                                                                                                                                gCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCATTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACAt                                                                    <  1:327670/251‑1 (MQ=255)
                                                                                                                                                                gCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACGGGCGTTAAGAACTTACAAAGCATGCTGATGAACGTCAAGCTACGTGGAATTATACTGCTGGCGGACGGTAAATTCTGTCCCTAAAGAAGTGCGATGGCACTGACGGGTCCGGAAACGGAGGGAGCGTTGCAGGGCGTCTTTCGAGGTAGTACTGGCGATTTCACTGTAACCGACCCGGACGGTGCGCGGAGACGCATCAATACGcac                                                                      >  2:374567/1‑249 (MQ=255)
                                                                                                                                                                      cccAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATccgccg                                                              <  1:252770/250‑1 (MQ=255)
                                                                                                                                                                       ccAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCcgccgc                                                             >  1:307214/1‑251 (MQ=255)
                                                                                                                                                                         agcgtgaaCGATGAGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAACGAGCTGCTGTACGTAAAGCTGCGTGAAATTATACTGCTGGCGGCCGGTCAATTCTGTCCCTAAAGAAGTGAGCGTGCCCTGACGGGTCAGGAAACGGAAGGGGCGTTGCAGGGGCTCTCTCTCGGTAGGAAGGGCGATGTCACCGTACCCTAACCGGACGGGGCGCTGAGGCGCAACAGAACgaaaaacaacagcgc                                                           >  2:341050/8‑237 (MQ=255)
                                                                                                                                                                         aGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCCAAACGCACATCCGCCgcgc                                                           >  1:65069/1‑251 (MQ=255)
                                                                                                                                                                               cACGAGGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTTCTGCTGAACGTCAAGCTGCTGGAAATTATCCGGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATGTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATTAAAACGCACATCCGCCGCGCTTTggg                                                     <  1:360879/251‑1 (MQ=255)
                                                                                                                                                                                      cGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAg                                              <  1:363578/251‑1 (MQ=255)
                                                                                                                                                                                             gcgcGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGGCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGCCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTCCGCGGTATTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGc                                       <  1:175035/251‑1 (MQ=255)
                                                                                                                                                                                                cgGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAgg                                    <  1:20875/251‑1 (MQ=255)
                                                                                                                                                                                                  gtggCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAGGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCATAGCCGGACGGTGCGCTGGTGCGCATCAAAACGCACCTCCGCCGCGCTTTGGGTATTCAGCGTGCGCAgggg                                  >  2:115076/3‑250 (MQ=255)
                                                                                                                                                                                                      cACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCa                              <  1:373097/251‑1 (MQ=255)
                                                                                                                                                                                                           gCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTc                         <  1:217546/251‑1 (MQ=255)
                                                                                                                                                                                                                aaGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGCCGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGACGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCACTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGAGAGCCATCTTCCAGtt                    >  1:301115/1‑251 (MQ=255)
                                                                                                                                                                                                                       tCCACAGCGTGCTGCTGAACGTCAAGCTGCGGGAAATTTGCCTGTGGGCGGACGGGAAATTATGTCACTAAAGTAGTGAGCTGGCACTGTCGGGTCAGGACACGGAAGGGGCGTTGCAGGGCAGCTTGCGCGGTATTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGTCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGAc              <  1:305599/250‑1 (MQ=255)
                                                                                                                                                                                                                          acaGCGTGATGCTGAAAGTCAAGCTGCGTGACATTATCCTGCTGGCGGACGGTACATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCCGGACACGGAAGGGGCGTTGCAGGGCATCTTACGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGGATCAAAACGCACCTCAGCCGCGCTTTGGGCATTCAACGGGCGCAGGGAGCCATCTTTCAGTAGCAGACGGCt          >  2:295345/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ctgctgaACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGc  >  1:331671/1‑251 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
CGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGC  >  NC_000913/4516895‑4517371

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