Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,440,768 G→A L150L (CTG→TTG)  paeA ← putative cadaverine/putrescine exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,440,7680GA95.0% 173.5 / ‑1.4 60L150L (CTG→TTG) paeAputative cadaverine/putrescine exporter
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (27/30);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

CATCTTCGTGGAGATCAAACCAAATCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCACCGACGATACCGCCGAGAATCG  >  NC_000913/4440525‑4441010
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
ccttctcggggagacaaacccAAATCCCGTTTCCACGCGGTGCTATTAAAGACGGAACGGTACGGGATTCCAGCTAAAAGACGTTTTCAATCACCTCGTGTTCCGGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTg                                                                                                                                                                                                                                              <  2:343763/232‑1 (MQ=255)
                 acccAAACCACTTTTTCACGCGGTGTCATTGAAGACGGAACGTTACGGGATTCCAGCACAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGaa                                                                                                                                                                                                                            <  1:413613/249‑1 (MQ=255)
                   ccaaatccagtttccccggggTGTCATTGAAAAGGGACCGGTACGGACTTCCAGCTCAAAGACTTTTTCAATCAGCTCGTGTTCCTGTTCTCGTAACCCGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAgc                                                                                                                                                                                                                          <  2:459095/232‑1 (MQ=255)
                    cAAATCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAgcg                                                                                                                                                                                                                         >  1:145744/1‑251 (MQ=255)
                      aaTCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAgcgc                                                                                                                                                                                                                        >  1:277559/1‑250 (MQ=255)
                          cgtttttcccgcggtgatttggAAGACGGAACGTTACGGGATTCTAGCACAAAGGCGGTTTCATCTAGCTCTGGTTCCTGTTTACGTAACCCGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTGCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCgg                                                                                                                                                                                                                   <  2:421030/230‑1 (MQ=255)
                          ccGTTTTCACGCGGTGTAAGTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGGCATCTTTACGTACCATTGGCAGGTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCgg                                                                                                                                                                                                                   <  1:332768/250‑1 (MQ=255)
                          actttctCACGCGGTGTCATTGAAGACGAAACGGCACGGGCTTCCACCTCAAAGCCGTTTTCAATCAGCCTGTGTTCCTGTTGACGTAACCCGCCCGCCAGCGCACCGGCTTCCACTGCCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCgg                                                                                                                                                                                                                   <  1:193708/245‑1 (MQ=255)
                                    gcgGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCACGCCAGCGCACCGGCTTCCACTACCGCGTAGATGACATCAGAAGTGATGTCATCTTTACGTACGATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAAAGTACAGGCAGAAGCGCATCGGGTTgatgat                                                                                                                                                                                                          >  1:312981/1‑250 (MQ=255)
                                                 aGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCcac                                                                                                                                                                                              >  1:307432/1‑249 (MQ=255)
                                                   aCGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCt                                                                                                                                                                                          <  1:337343/251‑1 (MQ=255)
                                                    cgaaaCGGTTCGGGATTCCAGCTCAAAGACTTTTTCAATCAGCTCGGGTTCCTGTTTACGTAACACCCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATTTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGTCAGAAGCGCATCGGTTTGATGATACGCAAAGCCACAGCtt                                                                                                                                                                                         <  1:131918/248‑1 (MQ=255)
                                                     ggAACGGTACGGGTTCCCAGCCCAAAATCGTTTCCAACCAGCTCGTGTTCCTGTTTCCGTACCCCCCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTc                                                                                                                                                                                        <  2:349165/251‑1 (MQ=255)
                                                      gAACGGTCCGGGATTCCAGCTCAAAGACGTTTTCAATCATCTCGAGTTCCTGTTTACGTAACACGCCCCCCAGCGCACCGGCTTCCACTACCACGTAGATGCCCTCAAAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAATATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAAAGTGGTCACAACGTAAAGGCAGAAGCGAATCGGGTTGATGATCCGCAAAGCCACAGCGTCt                                                                                                                                                                                       >  1:290552/1‑251 (MQ=255)
                                                           gTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGACGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGAACCGAAAAGCCACAGCTTCTGgcg                                                                                                                                                                                   >  2:772/1‑250 (MQ=255)
                                                            tACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTCACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGAACAGGCAGAAGCGCATCGGGTTGATGATACGCCAAGCCACAGCTTCTGgcgc                                                                                                                                                                                  >  1:104529/1‑250 (MQ=255)
                                                                    ccacgccaaacgccGTTTTCAACAAGCTCGTGGTCCTTTTTACTGAACACGCCCGCCAGCGCACCGGCTTCCACTGCCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAACCATAcc                                                                                                                                                                         <  2:212060/238‑1 (MQ=255)
                                                                               gACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGACCGCCAGAGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGAACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGAGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGTTGAGACGCAAAGCAACGGCTTCTGGAGCACTTATACCgtggagtttcg                                                                                                                                                              >  2:86494/1‑241 (MQ=255)
                                                                                 ccggttccaaccacctcggtttcgtGTTTACGGAACCCGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGGTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCgg                                                                                                                                                             <  2:19903/227‑1 (MQ=255)
                                                                                       cAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCAACGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTaa                                                                                                                                                       <  1:82795/250‑1 (MQ=255)
                                                                                                   tgtTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACCGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAAc                                                                                                                                             >  1:486596/1‑249 (MQ=255)
                                                                                                       ccTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAATCCACCGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATACGCAAACAGGAt                                                                                                                                        >  2:248004/1‑250 (MQ=255)
                                                                                                         tGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATaa                                                                                                                                      >  1:394726/1‑250 (MQ=255)
                                                                                                           tttACGTAACCCGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAAc                                                                                                                                    <  2:31403/250‑1 (MQ=255)
                                                                                                               cGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGcc                                                                                                                               >  1:127116/1‑251 (MQ=255)
                                                                                                                cagtcacctacgccaggcaaCCGGCTTCCACTACCGCGTAGATGTCATAAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTACGGTCGCCAGGCCGTAGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTGTCGGGGTGAAAGCCGCAAACAGGACAATCATGCCa                                                                                                                              <  2:67322/232‑1 (MQ=255)
                                                                                                                    cacGCCCGCCCGCGCACCGGCTTCCACTACCGCGTCGATGTCATCAGCAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGGGCAAACGTAAAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCACTCATACCGAGGCGTTTCGGGGTTAAATCCGCACACAGGATAAACATGCCAGTc                                                                                                                           >  2:239177/1‑250 (MQ=255)
                                                                                                                           gccAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGAGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACg                                                                                                                    >  1:366715/1‑250 (MQ=255)
                                                                                                                            ccAGCGCACCGGCTTCCCCTACCGGGTAGATGTCAGTAGAAGGGATTCCATCTTTACGTACCATTGGCAGTTTGAAAAGACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGAGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACCTGCCAGTCCCTAACgag                                                                                                                  <  2:401895/251‑1 (MQ=255)
                                                                                                                                       gCTTACACTACCACGTAGCTGTCATCAGAAGTGATGTCATCATAACGTACCATTGGCAGTTTAAAAATACGGAAGAATATGTTCGCCAGGCCGTTGAAGAACCACACCAAAGGGGGGAAAACGTACAAGCAGAATCGCATCGGGTTGAGGATACGCAAAGACAAAGCTTCTGGCACAACCATAACGATACTTTTCGGGGTTAAATCCGCAAACCGGATAAACATGCCAGTAACTAACgataatgcgagcat                                                                                                       >  2:111660/1‑239 (MQ=255)
                                                                                                                                           ccgcTACCGCGTATATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAg                                                                                                   <  1:386304/248‑1 (MQ=255)
                                                                                                                                                cccGCGTAGATGTCATCAGAAGTGAGTTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAg                                                                                              <  1:109484/250‑1 (MQ=255)
                                                                                                                                                  cgcgTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGtt                                                                                            <  1:415172/251‑1 (MQ=255)
                                                                                                                                                    cgTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCCAACAGGATAAACATGCCAGTCACTACCGCGAAAGAGAGAATAAAGCTCAGTTGc                                                                                          >  1:296282/1‑251 (MQ=255)
                                                                                                                                                     gTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCt                                                                                         <  1:248004/251‑1 (MQ=255)
                                                                                                                                                           gTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGACGATTATGTTCGCCAGGCCGTTGAAGAAACACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAGGCCACAGCTTCTGGCGCAATCACACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATTCCAGTCACTAACGAGAAAGAGAGAATCAAGCTCGGTTGCTCag                                                                                      >  1:496509/1‑248 (MQ=255)
                                                                                                                                                                     gTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACCCCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGa                                                                         <  2:307432/251‑1 (MQ=255)
                                                                                                                                                                     gTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCg                                                                          >  1:21779/1‑250 (MQ=255)
                                                                                                                                                                       gATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTACCGGGAAAGAGAGCATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGAc                                                                        >  2:14813/1‑250 (MQ=255)
                                                                                                                                                                       gATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGAc                                                                        <  1:186988/250‑1 (MQ=255)
                                                                                                                                                                                    cGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCCCTAACGAGAAAGAGAGCATa                                                                                                      >  2:214489/1‑207 (MQ=255)
                                                                                                                                                                                    cGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATa                                                                                                      <  1:214489/207‑1 (MQ=255)
                                                                                                                                                                                        ccATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGgaga                                                           >  2:262271/1‑246 (MQ=255)
                                                                                                                                                                                        ccATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGgaga                                                           <  1:262271/246‑1 (MQ=255)
                                                                                                                                                                                        ccATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAgg                                                      <  1:448809/251‑1 (MQ=255)
                                                                                                                                                                                              gCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTa                                                                                                                       <  1:3441/180‑1 (MQ=255)
                                                                                                                                                                                              gCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTa                                                                                                                       >  2:3441/1‑180 (MQ=255)
                                                                                                                                                                                              gCAGTTTGAAAATACGGAAGATGATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGa                                                <  1:136216/251‑1 (MQ=255)
                                                                                                                                                                                                 gTTTGAAAATACGGAAGATTATGTTCGCAAGGCCGTTGAAGAACCAAAACAACGAGATGCCAACGAACAGGAAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGGGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGaaa                                              >  1:299559/1‑250 (MQ=255)
                                                                                                                                                                                                  tttGAAAAAACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGCGCAAACTTACAGGCAGCAGCGCATCGGGTTGATGAAACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGACCGTCACTAACGAGAAAGAGAGACTAAAGCTCAGTTGCTCAGAGAGCTCTGCCGCCATGTAGAGGGAGAACAGGCTGTGaaa                                              >  2:259319/1‑249 (MQ=255)
                                                                                                                                                                                                      aaaaTACTGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCGGGCGCAATCAGACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGAGAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCGCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTgg                                        <  2:305965/251‑1 (MQ=255)
                                                                                                                                                                                                       aactccGGAAGATTAGGTTCGCCAGGCCGTGGAAGAACCACCCCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAACCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATATCCATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCGCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGTGAACAGGCTGTGAAAAGCTGga                                       <  2:93762/246‑1 (MQ=255)
                                                                                                                                                                                                              gAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAAAATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGTACGGGCTGTGAAAAGCTGGAGAAAAtt                                >  2:337198/1‑250 (MQ=255)
                                                                                                                                                                                                               aaGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAAAACAGGATGTGAAAAGCTGGAGAAAATGc                               >  2:172418/1‑251 (MQ=255)
                                                                                                                                                                                                               aaGATGATGTTCGCCAGGCCGTAGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGc                               <  2:104469/251‑1 (MQ=255)
                                                                                                                                                                                                                aGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCa                                                                                                                              >  2:105921/1‑155 (MQ=255)
                                                                                                                                                                                                                aGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCa                                                                                                                              <  1:105921/155‑1 (MQ=255)
                                                                                                                                                                                                                     aTGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCAt                          <  1:301567/250‑1 (MQ=255)
                                                                                                                                                                                                                                   gTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCACCGACGATAccgc           >  1:121516/1‑251 (MQ=255)
                                                                                                                                                                                                                                            ccACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCACCGACGATACCGCCGCGAATCg  >  1:502389/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
CATCTTCGTGGAGATCAAACCAAATCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCACCGACGATACCGCCGAGAATCG  >  NC_000913/4440525‑4441010

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