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,1270CA93.5% 183.3 / 1.0 63V197V (GTG→GTT) fecRferric citrate regulator FecR
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (28/30);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

AGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCC  >  NC_000913/4516885‑4517372
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
aGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGGAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAg                                                                                                                                                                                                                                               <  2:365743/251‑1 (MQ=255)
     ggCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCg                                                                                                                                                                                                                                          >  1:46663/1‑251 (MQ=255)
              gggATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTAt                                                                                                                                                                                                                                 <  1:96847/251‑1 (MQ=255)
                   cgcagcgcagCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCATGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCctg                                                                                                                                                                                                                             <  2:218393/250‑1 (MQ=255)
                     cagcgcagcACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGCAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCCCCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGCAATTATCctgctg                                                                                                                                                                                                                          >  2:137885/1‑251 (MQ=255)
                                gCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTaa                                                                                                                                                                                                               <  2:378132/251‑1 (MQ=255)
                                       cGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGt                                                                                                                                                                                                        <  1:22981/251‑1 (MQ=255)
                                              gggTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTaa                                                                                                                                                                                                  >  1:4032/1‑250 (MQ=255)
                                                      gTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTg                                                                                                                                                                                           >  1:47931/1‑249 (MQ=255)
                                                        ggCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCt                                                                                                                                                                                       <  1:248580/251‑1 (MQ=255)
                                                            aTCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGcc                                                                                                                                                                                   >  1:160393/1‑251 (MQ=255)
                                                                         gtggTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAgg                                                                                                                                                                      <  2:136357/249‑1 (MQ=255)
                                                                             tttATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGCCGGGTCAGGgcac                                                                                                                                                                  >  2:98037/1‑251 (MQ=255)
                                                                                       aaGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGTCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAgggggg                                                                                                                                                        >  2:283594/1‑249 (MQ=255)
                                                                                          cTCAGGAAGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCTCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCCCTTCCACCGCGTGCTGCTGCAC‑GTCAAGCTGAGTGAAATTATCCTGCTGGCGGACCGTAAATTATGTACCTAAAGCCGTGCTCTGGCCCTGACGGGTCCGGACACGGAAGGGGCggtg                                                                                                                                                     >  2:8541/1‑248 (MQ=255)
                                                                                            ccggattcacttggTCCAGCTTGCACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGAGTTCTCCAGCGTTCACGATGTGTTTTTGCGCGGGGGCACTTGCGAGAATCACGTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGATACGGAAGGGGTGTTGc                                                                                                                                                    <  1:494763/237‑1 (MQ=255)
                                                                                            cAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCt                                                                                                                                                                                 >  1:208861/1‑221 (MQ=255)
                                                                                            cAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTg                                                                                                                                                     <  2:143960/249‑1 (MQ=255)
                                                                                            cAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTg                                                                                                                                                     <  1:152161/249‑1 (MQ=255)
                                                                                            cAGGATGTCCTTCGGCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCt                                                                                                                                                                                 <  2:208861/221‑1 (MQ=255)
                                                                                               gATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGGAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGCAAGGA‑ACGGAAGGGGCGGTGccaggg                                                                                                                                               >  2:451864/1‑246 (MQ=255)
                                                                                                aTGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGc                                                                                                                                               >  1:126802/1‑251 (MQ=255)
                                                                                                    ccTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCAC‑GTCAAGCGGCGTGAAATTATACTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTAAGGACACGGAAGGGGCGTTGCAGGGCATc                                                                                                                                            <  1:24332/250‑1 (MQ=255)
                                                                                                    ccTTCGTCCAGCTTGTACTCCCGTCCTCCCGCGGTTTCCCGGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCt                                                                                                                                           <  2:402186/251‑1 (MQ=255)
                                                                                                       ccGTCCAGCTTGTACTCCCGTCAGCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCtt                                                                                                                                          <  2:128188/248‑1 (MQ=255)
                                                                                                         gTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGCCGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTcgc                                                                                                                                      >  2:425544/1‑251 (MQ=255)
                                                                                                                         cGTCATCCAGCGGTTTCACTGCGCCAATATCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCg                                                                                                                         <  1:374369/248‑1 (MQ=255)
                                                                                                                         cGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGAtt                                                                                                                      <  1:284644/251‑1 (MQ=255)
                                                                                                                          gTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGAttt                                                                                                                     <  2:46663/251‑1 (MQ=255)
                                                                                                                            cATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTcc                                                                                                                   >  1:14759/1‑250 (MQ=255)
                                                                                                                             aTCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTcgcgaggtcggaagagcgtcggg                                                                                                                  >  2:239056/1‑232 (MQ=255)
                                                                                                                             aTCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTcgcg                                                                                                                                     <  1:239056/232‑1 (MQ=255)
                                                                                                                               ccAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCGCCAGCGTTCGCGAGGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGGCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCg                                                                                                                <  2:47931/251‑1 (MQ=255)
                                                                                                                                cAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGt                                                                                                               >  1:224618/1‑251 (MQ=255)
                                                                                                                                    ggTTTCACTGCGCCAGACTCAGTGGCGCTGAACTGCATGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCCCTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGATGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCa                                                                                                           <  2:115031/251‑1 (MQ=255)
                                                                                                                                     gTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGCGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGCGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGtccc                                                                                                            >  2:355920/1‑246 (MQ=255)
                                                                                                                                     gTTTCAATGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCAATGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGATGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTTCCag                                                                                                          <  2:353517/251‑1 (MQ=255)
                                                                                                                                       ttCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCaga                                                                                                         <  1:464096/250‑1 (MQ=255)
                                                                                                                                         cACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTTAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGcc                                                                                                      <  1:81844/251‑1 (MQ=255)
                                                                                                                                             gcgcCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGa                                                                                                   <  2:160393/250‑1 (MQ=255)
                                                                                                                                               gcCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCCGGGAATCTTTCGCGGTAGTAATGGCGATTTACCCGTAACAAGGCCGGACgg                                                                                                >  2:432091/1‑251 (MQ=255)
                                                                                                                                                  aaaCTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTACCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTgc                                                                                             <  1:385868/251‑1 (MQ=255)
                                                                                                                                                         gagGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCCGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGAtt                                                                                      >  2:110712/1‑250 (MQ=255)
                                                                                                                                                         gagGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATGCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATg                                                                                      <  1:321348/251‑1 (MQ=255)
                                                                                                                                                           ggCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGCTgcg                                                                                    >  2:416391/1‑251 (MQ=255)
                                                                                                                                                                      gCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAACTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGc                                                                         <  1:95782/251‑1 (MQ=255)
                                                                                                                                                                                 ccAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGCGCCGGACGGTGCGCTGATgc                                                                                     >  2:260205/1‑228 (MQ=255)
                                                                                                                                                                                 ccAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATgc                                                                                     <  1:260205/228‑1 (MQ=255)
                                                                                                                                                                                       ttCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGTATCAAAACGCACATCCGCCGCGCTTTg                                                        >  2:337984/1‑251 (MQ=255)
                                                                                                                                                                                                 gTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGc                                              >  1:313453/1‑251 (MQ=255)
                                                                                                                                                                                                     gtGCGGGGGGGCACTGGCGAGAAGCACTTCCACAGCGGGCTGCTGAAC‑GTCAAGTGGCGGGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAGGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGAGGGTGCGCTGATGCGCATCAAAACGCCCATCCGCCGCGCTTTGGGTATTCAGCGTg                                           <  2:390201/249‑1 (MQ=255)
                                                                                                                                                                                                          cgGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAgg                                     >  1:230721/1‑251 (MQ=255)
                                                                                                                                                                                                          cgGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGGTTTCAGCGTGAGCAgg                                     >  2:442045/1‑251 (MQ=255)
                                                                                                                                                                                                              ggCACTGGCGAGAAGCACTTCCCCAGCGGGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGc                                 <  2:224618/251‑1 (MQ=255)
                                                                                                                                                                                                              ggCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGc                                 >  1:269/1‑251 (MQ=255)
                                                                                                                                                                                                                   tGGCGAGAAGCACTTCCAAAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCt                            >  1:112973/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ccACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACgg             >  1:233069/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ccACAGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACg              <  1:416391/250‑1 (MQ=255)
                                                                                                                                                                                                                                   cacaGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGc            >  1:24643/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cacaGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGGCAGGACACGGAAGGGGCGTTGCAGGGAATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGCTGCGCATCAAAACGCACATCCGCCGCGCTTTGGGAATTCAGCGTGAGCAGGGAGCCATCCTCCAGTTGCTGACGGc            >  2:289816/1‑251 (MQ=255)
                                                                                                                                                                                                                                    acaGCGTGCTGCTGAAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCCCCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAAAGCACATCCGCCGCGCTTTGGGAATTCAGCTTGAGCAGGGAGCCATCTTCCAGTTCCTGACGGCt           >  2:304104/1‑251 (MQ=255)
                                                                                                                                                                                                                                     caGCGCGCTGCTGACC‑GTAAAGCTGCGTGAAATTATACTGTTGGCGGACGGAAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCGGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGAGGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCAGCTTCCAGTTGCTGACGGCTg          <  2:134025/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gTGCTGCTG‑ACAGTCAAGCTGCGTGAAATTATCCTGCTGCCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGAGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCCAAACGCACCTCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCCTCTTCCAGTTGCTGCCGGCTGACtg      >  2:457168/1‑251 (MQ=255)
                                                                                                                                                                                                                                             tgctgaAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGcc  >  1:20870/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
AGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCAC‑GTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCC  >  NC_000913/4516885‑4517372

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑