Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,516,378 G→A N100N (AAC→AAT fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,516,3780GA96.5% 161.6 / ‑3.4 60N100N (AAC→AATfecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (27/29);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

GATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTC  >  NC_000913/4516132‑4516621
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
gATGCCAAAGTTCAGGCCCAGGGCGTGGCTGCCGGTGCCGTTGTTGTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCtattat                                                                                                                                                                                                                                                 <  2:516/251‑1 (MQ=255)
          ttCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGCTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATAAAGTCACACGCGCCAGCCTGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCCCCGACCACGGTCGCGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCCGCTATTAGTTCCCAGctg                                                                                                                                                                                                                                       >  1:235171/1‑249 (MQ=255)
           tCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGATAAGCACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGCACCCGACACGCGCGCCAGCAAGTTCAAATACGTCGTTTTCACGCCCATTACCCCGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGGGCAGGCGCCGGCTCCCGCGCCCAGCTATTTTTTCCCGGCGGt                                                                                                                                                                                                                                      >  2:389443/1‑251 (MQ=255)
           ccagcgccggggcgcggCTGCCGGTGCCGTTGTTTTCCGGGGCGCTGACGCCAGGGATGCGGTTAAGGACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCTCCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCTGGCGCCAGCGTCCAGCTATTATTTCCCAGCGGt                                                                                                                                                                                                                                      <  2:289264/236‑1 (MQ=255)
            cATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGt                                                                                                                                                                                                                                      <  1:436068/250‑1 (MQ=255)
            cATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCAGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGACAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTCTTATTTCCCAGCGGtt                                                                                                                                                                                                                                     >  1:460756/1‑251 (MQ=255)
               cGCCAGGTCGTGGCTGCCGGTGCCTTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCCCGCATGGTGGTTGCGCCGGTGTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTGCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTGTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTa                                                                                                                                                                                                                                  <  2:436333/251‑1 (MQ=255)
               cGCCAGGTCGTGGCTGCCGGTGCCTGTGTTTTCCGGCGCGCGGGCGCCAGGGATGCGGTTAAGTTCCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAGTCCGCACGGCGGAGCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACTGTCAGGGCATCTGCTTTTGGTGCGGGCGCGGGCACCAGCGTCCAGCTATTATGTCCCAGCGGTTTTa                                                                                                                                                                                                                                  <  2:498510/251‑1 (MQ=255)
                      ccttGGCTGCCGGTGCCGTTGTTTTCCGGGTAGCTGACGCCCGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGAGCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAg                                                                                                                                                                                                                           <  2:143136/248‑1 (MQ=255)
                          ggCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCg                                                                                                                                                                                                                       <  1:66121/251‑1 (MQ=255)
                                      cGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTGTTGGCGAAATCCTCACGGCGGATCCCGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGGCCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCcagc                                                                                                                                                                                                           <  2:9643/251‑1 (MQ=255)
                                        ttgttTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATAAGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCCCGGTCAAGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGGCCCGCACTTATTTCCCAGCGGTTGTACCTTCAGTCCGCTGCCGTCcagc                                                                                                                                                                                                           >  2:341812/1‑249 (MQ=255)
                                                     cgcTGACGCCAGGGATGCGGTTAAGTAGCTCACGCATGGTGGTTGCGCCGGTTTTGGAGAAATAATCACAGCGGATCACGTCACGCGGGCCAGCATTTTCAACACCGTCGTTGTCACGCGCATCACCCAGCCAGTCGCCGACCAAGGTAATGGCATCTTCCTTTGGTGCGGGCGCGAGCTCCGGCGTACAGCTATTATTTCCCCTAGGTTTTACCTGCAGTCCGCTACCGGACAGCAGTTGTTGCCGGccc                                                                                                                                                                                            >  2:504886/1‑250 (MQ=255)
                                                     cgcTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCg                                                                                                                                                                                            <  1:2507/251‑1 (MQ=255)
                                                      gcTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCCGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCCGCCAGTCGCCGACCACGGACAGGGCATCTTCTTTTGGTGCGGGCCCGGGCCCCAGCGGCCCGCCATTTTTTCCCAGCGGGTTTACCTTCAGTCCCCTGCCGGCCCGCAGTTGTTTGAGGCCGc                                                                                                                                                                                           >  2:428805/1‑251 (MQ=255)
                                                                   aTGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCGTCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGAcgtcgt                                                                                                                                                                               >  1:64747/1‑250 (MQ=255)
                                                                       ggTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCCAATACGTCGTTTTCACGCGCATCACCCCGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCCGCGGTTTTACCTTCAGTCCGCTGCCGTCCCGCAGTTTTTGCAGGCCGGTTTCCACGTCGTAATcc                                                                                                                                                                          >  2:305804/1‑250 (MQ=255)
                                                                         ttAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCAAGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGcc                                                                                                                                                                        >  1:156451/1‑251 (MQ=255)
                                                                               aCCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCCGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCCGTCCGCTTCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGCCGTCGTAATCGCCGTGc                                                                                                                                                                    >  1:57357/1‑249 (MQ=255)
                                                                                        aTGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCCTCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGct                                                                                                                                                         >  1:64264/1‑251 (MQ=255)
                                                                                                  cgcCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCcg                                                                                                                                               <  1:341812/251‑1 (MQ=255)
                                                                                                   gcCGGTTTTGGCGCAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTACCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCCGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGtt                                                                                                                                                            >  2:319829/1‑237 (MQ=255)
                                                                                                   gcCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGtt                                                                                                                                                            <  1:319829/237‑1 (MQ=255)
                                                                                                       gTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGt                                                                                                                                          <  1:349251/251‑1 (MQ=255)
                                                                                                          ttGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCa                                                                                                                                        <  1:26132/250‑1 (MQ=255)
                                                                                                             gCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCGCCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTGACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCt                                                                                                                                    <  1:183157/251‑1 (MQ=255)
                                                                                                                aaaaccgccgggcggATCACGTCACGGGCGCCAGCATGTTCAAATCCGGCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCGTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTAGTATTTCCCAGCGGTTTGACCGGCAGTCCGCTGCCGTCCAGCAGTTGTGGCAGGCCGCTCTCGACGGCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCTGCGCGTCAGGCTg                                                                                                                                   <  2:57357/239‑1 (MQ=255)
                                                                                                                    ccTCACGGCGGATCACGGCACGCGCGCCAGCATGTTCAAATACGTCGTGTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCa                                                                                                                             <  1:366689/251‑1 (MQ=255)
                                                                                                                          ggcggATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCACCCgcg                                                                                                                       >  1:350379/1‑249 (MQ=255)
                                                                                                                             ggagcacctcacgCGCCCCGGCATGTTCAAATACGTCGTTTTCACGCGCCTCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAgg                                                                                                                    <  2:428450/243‑1 (MQ=255)
                                                                                                                                 cacgtcacgCGCGCCAGCATGTTCAAATACGTCGTTTTAACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAgggaa                                                                                                                 >  2:455045/1‑248 (MQ=255)
                                                                                                                                   cttcacgTGCGCCAGCATGTTCAAATACGGCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTAATTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGt                                                                                                                  <  1:57310/245‑1 (MQ=255)
                                                                                                                                    gtcacgCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCCCCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCGGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATc                                                                                                             <  2:64747/251‑1 (MQ=255)
                                                                                                                                    gtcacgCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATc                                                                                                             >  1:215647/1‑251 (MQ=255)
                                                                                                                                          cgcgcCAGCATGTTCAAATACGTCGTTATCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCtgt                                                                                                       <  2:336805/251‑1 (MQ=255)
                                                                                                                                                   aTGTTCAAATACGTCGTTTCCACGCGCACCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGGGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCAt                                                                                              <  2:224205/251‑1 (MQ=255)
                                                                                                                                                       tCAAATCCGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGGGAGGTGAAATCCGCTGTGGGCGGCATACTg                                                                                          >  2:217465/1‑251 (MQ=255)
                                                                                                                                                          aattcgccggtttCACGCGTATCTCCCAGCCAGGCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGGAAATCCGCTGTGTGCGGCATACTGAtt                                                                                       <  2:292502/241‑1 (MQ=255)
                                                                                                                                                                  gtTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGAGCtt                                                                               >  1:502283/1‑251 (MQ=255)
                                                                                                                                                                     ttCACGCGCATCTCCCAGGCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCATTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCATGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTa                                                                             <  1:68933/250‑1 (MQ=255)
                                                                                                                                                                        ccccgcATCCCCCAGCAGGTCGCCGACCACGGTCAGGGCATCTTCTGTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGa                                                                         <  2:370345/248‑1 (MQ=255)
                                                                                                                                                                        aCGCGCAGCACCCAGCCAGTCACCGACCATGGTCAGGGCTTCGTCTGTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTGCCCAGCGTTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTTTGCGGCATACTGATTGAGCGCTTTATCg                                                                          <  1:278161/250‑1 (MQ=255)
                                                                                                                                                                                cACCCAGCCAGTCTCCGCCCACGGTCAGGGCATCTTTTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATccc                                                                 <  2:156451/251‑1 (MQ=255)
                                                                                                                                                                                 aCCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCg                                                                >  1:445165/1‑251 (MQ=255)
                                                                                                                                                                                    cagccagTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTg                                                             >  1:463543/1‑251 (MQ=255)
                                                                                                                                                                                      gccagTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCg                                                           >  1:133635/1‑251 (MQ=255)
                                                                                                                                                                                      gccagTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCg                                                           <  1:197834/251‑1 (MQ=255)
                                                                                                                                                                                         agTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGACGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATATCGATTAAGCGATGTATCGAGCGATCCCGGTGCGat                                                         >  2:321348/1‑250 (MQ=255)
                                                                                                                                                                                               cGACCACGGACAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTACTCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAAcc                                                  >  1:341203/1‑251 (MQ=255)
                                                                                                                                                                                                 aCCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCtg                                                >  1:87803/1‑251 (MQ=255)
                                                                                                                                                                                                   cACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGACTGGCGTCAACCGCGAGGGTAAATCCGCTGTTGGCGGCATACTGATTGATGGCTTTATCGTGCGATCCCGGTGCGGTATTAACCtgtg                                              >  2:147865/1‑251 (MQ=255)
                                                                                                                                                                                                          aGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCaa                                        >  1:106368/1‑250 (MQ=255)
                                                                                                                                                                                                            ggCATCTGCTTTTGGTGCGGGCTCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAg                                     <  2:120494/251‑1 (MQ=255)
                                                                                                                                                                                                               aTCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGCCGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGCGAGGGTAAATCTGCTGTGTGCGGCATACTGGTTGGGCGCTTTATCGAGCGCTCCCGGTGCGATATTAAACTGTGCAGCAACAGCgg                                  >  2:489930/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAgg     <  1:201849/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCt   <  1:43914/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCt   >  1:8547/1‑251 (MQ=255)
                                                                                                                                                                                                                                                ccAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCt   >  1:82729/1‑249 (MQ=255)
                                                                                                                                                                                                                                                ccAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCt   <  2:82729/249‑1 (MQ=255)
                                                                                                                                                                                                                                                 cAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGTAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTc  >  1:109474/1‑249 (MQ=255)
                                                                                                                                                                                                                                                 cAGCTATTATTTCCCAGCGGTTTTACCTGCCGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGGGAGGGGAAATCCGCTGTGTGCGGCATACTGCTTGAGCGCTTTATCGCGCGATCCCGGTGCGATATTAACCTGTGCCGCAAAAGCGGAAAAAGAGCAACCGGCCAGCGGCAGCAGGccc  >  2:81844/1‑247 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
GATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTC  >  NC_000913/4516132‑4516621

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