Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,565,350 T→C 15.6% intergenic (+20/+129) yffS → / ← eutA CPZ‑55 prophage, uncharacterized protein/reactivating factor for ethanolamine ammonia lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,565,3500TC15.6% 155.0 / 2.3 70intergenic (+20/+129)yffS/eutACPZ‑55 prophage, uncharacterized protein/reactivating factor for ethanolamine ammonia lyase
Reads supporting (aligned to +/- strand):  ref base T (32/24);  new base C (3/9);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 5.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGACTGATTGAGCAGGCCAAAGGGAGAGGATTGATTAGCGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCG  >  NC_000913/2565125‑2565568
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                          
tGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGACTGATTGAGCAGGCCAAAGGGAGAGGATTGATTAGCGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGACACCAGATCGGCAAAGTCAAGAAAGGGGACAAGCTGATACCTTTGCTTCCCAACTGGATTGTTAAGCTGGAGGACAAATAAACCAGCCTTCAACCCCCtctc                                                                                                                                                                                                                           >  2:669577/1‑227 (MQ=255)
            gCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGACTGATTGAGCAGGCCAAAGGGAGAGGATTGATTAGCGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGAAAAGGCAAGAAAAGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGCAGAAAAAACAAGCCTTCAACCCCAACCCAccaaacaaaga                                                                                                                                                                                                               >  1:28575/1‑224 (MQ=255)
                  ggAAGCCAAACAGTAACGGCGCGTTTTATCGACTGATTGAGCAGGCCAAAGGGAGAGGATTGATTAGCGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGccc                                                                                                                                                                                                        <  2:76965/228‑1 (MQ=255)
                                    gacggTTTTATCGAATGATTGAGCAGGCCAAAGGGAGAGGATTGATTAGCGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCt                                                                                                                                                                                      <  2:281447/224‑1 (MQ=255)
                                     cgcgTTTTATCGACTGATTGAGCAGGCCAAAGGGAGAGGATTGATTAGCGTTGGACAGGAAATCGAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGGGACAACCTGATAACTTTGCTTCCCAACTGAATTGATGAGAATGGCGACGGAAAAACCAGCCTCTACCCCCCTCCACACAAACAACCACCGGCCCCCTCACCGGGGTTTTGcca                                                                                                                                                                                     >  1:182705/1‑226 (MQ=255)
                                     cgcgTTTTATCGACTGATTGAGCAGGCCAAAGGGAGAGGATTGATTAGCGTTGAACAGGAAATCAAACACAAGAAAGACGGAAACGGCAACCAGATCGGCAAAGGAAAGAAAGGAGACAAGCTGATAACTTTGCTTCCCAACTGAATTGATAAGCTGGAAAACGACAAAACCACGCTTACAACCCAACCCATAAACAAAAACCAGGCCCCACACCCCGGTTTTTGcaa                                                                                                                                                                                     >  1:269055/1‑226 (MQ=255)
                                     cgcgTTTTATCGACTGATTGAGCAGGCCAAAGGGAGAGGATTGATTAGCGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGAAAAGATGATAACTTTGCTTCCCAACTGAATTGATAAGCTGGGCGACGAATAAACCAGCCTTTAACCCCATATCCTCAATCAACTCCCGTCCCCGCTGCCGGGGTTTTGcct                                                                                                                                                                                     >  1:47028/1‑226 (MQ=255)
                                       cgTTTTATCGACTGATTGAGCAGGCCAAAGGGAGAGGATTGATTAGCGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGGGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAAGCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTAt                                                                                                                                                                                    >  2:230638/1‑227 (MQ=255)
                                           ttATCGACTGATTGAACAGGGCAAAGGGAGAGGATTGATTAGCGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGcac                                                                                                                                                                                <  2:203067/227‑1 (MQ=255)
                                             aTCGACTGATTGAGCAGGCCAAAGGGAGAGGATTGATTAGCGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCAGGCCCCGCCGCCGGGCCTTCGCTATGcacaa                                                                                                                                                                              >  2:450456/1‑225 (MQ=255)
                                                   tGATTGAGCAGGCCAAAGGGAGAGGATTGATTAGCGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGCTTGATAAGCTGGGAGACCAATAAACCAGCCTTCTACCCCATCTCCTCAATCACAGCCCGACCCCGCTGCCAGGTTCTTGCAATGCACCACCATAAc                                                                                                                                                                       >  2:321740/1‑228 (MQ=255)
                                                       tggacAGGCCAAAGGGAAGGGATTGATTAGCGGTGAACAAGAAATAAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCa                                                                                                                                                                   <  1:669724/224‑1 (MQ=255)
                                                       tGAGCAGGCAAAAGGGAGAGGATTGATTAGCGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCa                                                                                                                                                                   <  2:28575/228‑1 (MQ=255)
                                                              gCCAAAGGGAGAGGATTGATTAGCGTTGAACAGGAAATCAAACACAAGAAAGACGGAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAACCCAGCCTTCAACCCCCTCTCATCAATCAAAGCCCGGCCCCGCTTCCGGGTTTTTGCAATGCACCCCAATTACCCCACCCgggt                                                                                                                                                            >  1:396588/1‑226 (MQ=255)
                                                               ccAAAGGGAGAGGATTGATTAGCGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACACCAATTACCCCAACCGGATa                                                                                                                                                           >  2:218879/1‑228 (MQ=255)
                                                               ccAAAGGGAGAGGATTGATTAGCGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCCCTCCCGGGTTTTTCCTATGCACCAAAATAACCCCCACCGGaaa                                                                                                                                                           >  1:140370/1‑226 (MQ=255)
                                                                 aaggtgaggGGATTGAATAGGCTTGAAAAGGAAAGCAAAACAAAGAAAGACCAAAACGGCAACCAGATCGACAAAGGCAATAAAGGTTACATGTTGATAACTGTGCTTCCCAACTGGATTGATAAGCTGGGAGGCGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATaca                                                                                                                                                         <  1:618443/223‑1 (MQ=255)
                                                                        gagGATTGATTAGCGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCCCAATTACCCCAACCGGATACACAGCCgg                                                                                                                                                  >  2:253438/1‑228 (MQ=255)
                                                                              tgattAGCGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCCCAATTACCCCAACCGGATACACCCCCGGCTACa                                                                                                                                             >  2:188024/1‑227 (MQ=255)
                                                                               tattAGAGGTTAACAGGAAATCAAAAACAAGAAACACGAAAAGGGCAAACAGATCGGCAAAGGAAAGAAGGGGGACAAGCTTACAACATTGCTTCCCAACTGGATTGATAAGCTGGGCGACGAATAAACCAGCCTTCATCCCCATCTCAACAATCAACGCCCCGCCCCGCCGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAAt                                                                                                                                           <  1:321740/227‑1 (MQ=255)
                                                                               gattAGCGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACACTTACCCCAACCGGACACACAACCGGATACaa                                                                                                                                            >  2:676601/1‑227 (MQ=255)
                                                                                attAGCGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGATGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATAAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATAACCCCAACCCGCAACACAGCCGGATACAAt                                                                                                                                           >  1:297549/1‑227 (MQ=255)
                                                                                     cGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGCCGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCCCAATTACCCCACCCGGATACACAGCCGGATACAATTcccc                                                                                                                                      >  2:345454/1‑225 (MQ=255)
                                                                                     cGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCCACCGGATACACAGCCGGATACAATTccacc                                                                                                                                     >  2:148447/1‑228 (MQ=255)
                                                                                     cGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAAt                                                                                                                                           <  2:165634/222‑1 (MQ=255)
                                                                                     cGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAAt                                                                                                                                           >  1:165634/1‑222 (MQ=255)
                                                                                     cGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTccacc                                                                                                                                     >  2:553859/1‑228 (MQ=255)
                                                                                     cGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTccac                                                                                                                                      >  2:637204/1‑227 (MQ=255)
                                                                                     cGTTGAACAAGAAGAAAAACAAAAGAAAAACGAAAAAGGCAACCAGGTCGGAAAAGGAAAGAAAGGTGAAAAGCTGATAACTTTGCTGCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTccac                                                                                                                                      <  1:12854/227‑1 (MQ=255)
                                                                                     cGTTGAAAAAGAAATCAAACAAAAGAAAGGCGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTGCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTccac                                                                                                                                      <  1:72033/227‑1 (MQ=255)
                                                                                        taaaCAGGAAATCAAACAAAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAg                                                                                                                                   <  2:105038/225‑1 (MQ=255)
                                                                                        tGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAg                                                                                                                                   <  2:30033/227‑1 (MQ=255)
                                                                                          aaCAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCCCCAGcac                                                                                                                                >  2:619757/1‑228 (MQ=255)
                                                                                               gAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGGTTTTGCTATGCACCCCAATTACCCCAACCGGATACACAACCGGATACAAATCCACCAGcacccaac                                                                                                                           >  1:111316/1‑226 (MQ=255)
                                                                                                             gaaagaCGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGAGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACAAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCACCCACCCAGCGCAAc                                                                                                              >  2:243831/1‑227 (MQ=255)
                                                                                                              aaataCCAAAACGGTCACCATATCGTGAAAGGCAAGAATGGGGACAGCCGAAAAACTTTGCTTCCCAACCGGATTGTGAAGCTGGGAGACGAATAAGCCATCCTTCAACCCCATCTCAACAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCg                                                                                                            <  1:275117/228‑1 (MQ=255)
                                                                                                              aaagaGGAAAACGGCAACCAGATCGGCAAAGGGAAGAAAGGTGACAAGCCGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCg                                                                                                            <  2:269055/228‑1 (MQ=255)
                                                                                                              aaagaCGAAAACGGAAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCg                                                                                                            <  2:253512/228‑1 (MQ=255)
                                                                                                              aaaagcGAAAACGGCAACCAGAACGGCAAAGGGAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCg                                                                                                            <  2:182705/223‑1 (MQ=255)
                                                                                                                agccaaaaaCGGCAACCAGATAGGGAAGGGAAGAAAAGGTGACGAGGTGATAACTTTGCTATCACACTGTAATGCAAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAAACAACGCCTGGCCCCGCTGGCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCg                                                                                                            <  1:692869/221‑1 (MQ=255)
                                                                                                                 gaagAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGAACAGCTGATAACTTTGCTTCCCAACGGGATTGATAAGCTGGGAGACGAATAACCCAGCCTTCAACCCCATCCCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTCCCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTg                                                                                                         <  2:140370/225‑1 (MQ=255)
                                                                                                                     aaaaCGGCAACCAGATCGGAAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACACCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACaa                                                                                                                                            >  2:502356/1‑189 (MQ=255)
                                                                                                                      aaaccggatacAGATCGGCAAAAGGAAGAAAGGGGACAAGATGAAAACTTTGCTTCCCAACTGGATTGAGAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCCCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCGaa                                                                                                    <  1:188024/218‑1 (MQ=255)
                                                                                                                       aaccgtaaCCAGTGCGGCAAAGGCAAAAAAGGTGACAAGCTTGATACTTTGCTTCCAGACTTGACTGACAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCCAGCCACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGc                                                                                                       <  1:129111/218‑1 (MQ=255)
                                                                                                                          ggCAACCAGATCGGAAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCGAATa                                                                                                  <  2:459510/226‑1 (MQ=255)
                                                                                                                          cgcaacctcgccGGCAAAGGGAAGAAAGGTGACAAGCTGATAACGCTGCTTCCCAACTGGATTGATAAGCCGGGAGCCGAATAAACCAGCCTTCAACCCCATCCCATCAAGCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCGAATa                                                                                                  <  1:669577/215‑1 (MQ=255)
                                                                                                                             aacccaaTCGGCAAAGGCAAGAAAGAGTACAAGCAGATAACTTTGCTTCCCACCTGGATTGATAAGCGGGGAGACGAACAAACCAGCCTTCAACCCCATCACACCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTCCCCCAACCGGATCCACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCGAATACCGCa                                                                                             <  1:218879/222‑1 (MQ=255)
                                                                                                                             aaCCAGAACGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCGAATACCGCa                                                                                             <  2:297549/228‑1 (MQ=255)
                                                                                                                               cccaagcGGCAAAGGCAAGAAAGGGGACAAGCTGATAACTTTGCCTCCCAACTGGATTGATAAGCTGGGAAACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCGAATACCGCa                                                                                             <  1:539490/220‑1 (MQ=255)
                                                                                                                                 aggtCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCGCCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCGAATACCGCATTCa                                                                                         <  2:47028/225‑1 (MQ=255)
                                                                                                                                             gCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCCACCCCATCTCATCAATCAACGCCAGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATAACCCCCACCGACTACACAGCCGCATACAATTCCACCAGCACCCACCCACCCAGCGCCCCCCCTAGCGAAAACCCCATCCAGGAAGGaaaga                                                                              >  1:346089/1‑225 (MQ=255)
                                                                                                                                                        gACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTCCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCAACCAAGCGCCACCGCTGGCGCATACAGCATTAAGGAAGGATATGCGCGTGAtttt                                                                   >  2:465206/1‑226 (MQ=255)
                                                                                                                                                                cgATAACTTTCCTTCCCCACTGGATTAGTAAACTGGGAGACGAAAAACCCAGCCTTCAACCCCACCACATCAACCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTtcaccgtcac                                                           <  1:553859/227‑1 (MQ=255)
                                                                                                                                                                  aTAACTTTGCTTCCCAACTGCATTGATAAGCTGGGAGACGAATAACCCAGCCTTCCACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATCCACCACAATTACCCCAACCGGTTCCACAACCGGCACCAATTCCACCAGCACCCAGCCCCCCAGCACCCACGCTGACACAAACCGCATTCCGCACAGAAAACCGAGTCCTTtccccgtcacc                                                          >  1:504182/1‑227 (MQ=255)
                                                                                                                                                                         tGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCGATTACCGCATTCAGGAAGGAAATGCGAGGGATTTCACCGTCCCCGGCACAAc                                                  >  2:576468/1‑228 (MQ=255)
                                                                                                                                                                         tGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCCATGCACCACAATTACCCCAACCGGATACACACCCGGCTACAATTCCACCAGCACCAAGCCACCCAGCGCAACCGCCGGCGAAAACCCCATCCACAAAGCAAATACGACAGATTtcaccgtcaaccgacaaac                                                  >  1:229128/1‑218 (MQ=255)
                                                                                                                                                                         tGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCCGCCGGGTTTTTGCTATGCACCCCAATTACCCCAACCGGCTACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCCAATACCCCACTCACAACAGCAACGCGCGAGACTTCACCGCCACCGGCACAAc                                                  >  1:3426/1‑228 (MQ=255)
                                                                                                                                                                                ccaccggcttattAAGCTGCGAGACGGATTACCCCGCCTTTAACCCCATCCACACCATCAACGCCCCGCCCCGCTGCCGGGTTCTTGCTATCCACCCCAATGACCCCAACCGGAGACCCAGCCGGATACAATTCCAACAGCACCCAGCCACCAAGCGTCACCGCTGGCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCg                                           <  1:293624/216‑1 (MQ=255)
                                                                                                                                                                                      gATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCGAATACCGCATTCAGGAAGGAAATGCCAGTGATTTCACCGTCACCGGCACAACCGGTCCGCCcaaa                                     >  2:540746/1‑228 (MQ=255)
                                                                                                                                                                                          gaaaaGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCCCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCCCCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAagag                                 <  2:500791/225‑1 (MQ=255)
                                                                                                                                                                                           ttagcgcggggggcGAATAAACCGCCCTTCAACCCCACCCCAGCAAACAACGCCCAGGCCCGAGGCGGAGGTTTTTCTATGCACCACACTTACCCCACCCGGAGACACAGCCGGAGACAATTCCACCAGCACCCATCCACCCAGCGCCACCGCTGGCGAATACCTCAGTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAgg                                <  1:219374/215‑1 (MQ=255)
                                                                                                                                                                                                 tcggAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCACCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGAATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCCCCACCGCTGGCGAATACCGCATTCAGCAAGCAAATGCGCGAGATTTCACCGTCCCCCGCCCCACCGATCAGCCAAAAAGAggatga                            >  1:265269/3‑222 (MQ=255)
                                                                                                                                                                                                 cggggGACCAATAAACCAGCCCTCACCCCCACCCCACACATCCACCCCCGGCCCCCCTGCCGGGGTTTTTCTATGCCCCACAATTCCCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCCCCGCTGGCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTAcc                          <  1:676601/227‑1 (MQ=255)
                                                                                                                                                                                                  gggAGACGAATAAACCAGCCTTCAACCCCCTCCCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCCCAATTCCCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAagag                                 <  2:265269/220‑1 (MQ=255)
                                                                                                                                                                                                                   gCCTTCAACCCCATCTCCTCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCGAATACCGCATTCAGGAAGGAAATGCGAGTGCTTTCCCCGTCCCCGGCCCAACCCCACCGCCAAAAAGAGGCGGACCCATGTCGAAAACGTcccc        >  1:458035/1‑228 (MQ=255)
                                                                                                                                                                                                                    ccTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAAGCGTACCAATGTCGCTCTAGCcaccc       >  1:149178/1‑228 (MQ=255)
                                                                                                                                                                                                                         ccccccATCTCACCAACCACCCCCGGGCCCCGCGGCCGGGTTTTTGCTATGCCCCCCAATTCCCCCAACCGGATACACCGCCGGATACAATTCCACCAGCACCCAGCCCCCCAGCGCCTCCGCTGGCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCcgcgcg  <  1:576468/226‑1 (MQ=255)
                                                                                                                                                                                                                         acccccAACCCACCAAACAACCCCCGGCCCCGCTGCCGGGTTTTTGCCATGCACCACACTCACCCCAACCGGATCCACGGCCGGATACAAGCCCAACAGCACCCAGCCACCCAGCGCCCCCGCTGGCGAATACCGCAGTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCcgcg    <  1:243831/224‑1 (MQ=255)
                                                                                                                                                                                                                         aaCCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCCGGCGAATACCGCATCACGGAAGGAAATCCGAGGGATGTCCCCGCCACAGGCCCACCCCTCCAGCCAAAAAAAAGCGAACAAAAGTCAAAAAAATCCcacgcacg  >  1:324898/1‑228 (MQ=255)
                                                                                                                                                                                                                              cATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACCAtt                                                                                                                                          >  1:421596/1‑87 (MQ=255)
                                                                                                                                                                                                                               aTCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTc                                                                                                                                         <  2:421596/87‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                          
TGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGACTGATTGAGCAGGCCAAAGGGAGAGGATTGATTAGCGTTGAACAGGAAATCAAACACAAGAAAGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCG  >  NC_000913/2565125‑2565568

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