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

CATCTTCGTGGAGATCAAACCAAATCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATC‑AGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCACCGACGATACCGCCGAGAATCG  >  NC_000913/4440525‑4441010
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
ccttctcggggagacaaacccAAATCCCGTTTCCACGCGGTGCTATTAAAGACGGAACGGTACGGGATTCCAGCTAAAAGACGTTTTCAATC‑ACCTCGTGTTCCGGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTg                                                                                                                                                                                                                                              <  2:343763/232‑1 (MQ=255)
                 acccAAACCACTTTTTCACGCGGTGTCATTGAAGACGGAACGTTACGGGATTCCAGCACAAAGACGTTTTCAATC‑AGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGaa                                                                                                                                                                                                                            <  1:413613/249‑1 (MQ=255)
                   ccaaatccagtttccccggggTGTCATTGAAAAGGGACCGGTACGGACTTCCAGCTCAAAGACTTTTTCAATC‑AGCTCGTGTTCCTGTTCTCGTAACCCGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAgc                                                                                                                                                                                                                          <  2:459095/232‑1 (MQ=255)
                    cAAATCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATC‑AGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAgcg                                                                                                                                                                                                                         >  1:145744/1‑251 (MQ=255)
                      aaTCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATC‑AGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAgcgc                                                                                                                                                                                                                        >  1:277559/1‑250 (MQ=255)
                          cgtttttcccgcggtgatttggAAGACGGAACGTTACGGGATTCTAGCACAAAGGCGGTTTCA‑TCTAGCTCTGGTTCCTGTTTACGTAACCCGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTGCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCgg                                                                                                                                                                                                                   <  2:421030/230‑1 (MQ=255)
                          ccGTTTTCACGCGGTGTAAGTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATC‑AGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGGCATCTTTACGTACCATTGGCAGGTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCA‑AACGTACAGGCAGAAGCGCATCgg                                                                                                                                                                                                                   <  1:332768/250‑1 (MQ=255)
                          actttctCACGCGGTGTCATTGAAGACGAAACGGCACGGGCTTCCACCTCAAAGCCGTTTTCAATC‑AGCCTGTGTTCCTGTTGACGTAACCCGCCCGCCAGCGCACCGGCTTCCACTGCCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCgg                                                                                                                                                                                                                   <  1:193708/245‑1 (MQ=255)
                                    gcgGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATC‑AGCTCGTGTTCCTGTTTACGTAACACGCACGCCAGCGCACCGGCTTCCACTACCGCGTAGATGACATCAGAAGTGATGTCATCTTTACGTACGATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AAAGTACAGGCAGAAGCGCATCGGGTTgatgat                                                                                                                                                                                                          >  1:312981/1‑250 (MQ=255)
                                                 aGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATC‑AGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCcac                                                                                                                                                                                              >  1:307432/1‑249 (MQ=255)
                                                   aCGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATC‑AGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCt                                                                                                                                                                                          <  1:337343/251‑1 (MQ=255)
                                                    cgaaaCGGTTCGGGATTCCAGCTCAAAGACTTTTTCAATC‑AGCTCGGGTTCCTGTTTACGTAACACCCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATTTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGTCAGAAGCGCATCGGTTTGATGATACGCAAAGCCACAGCtt                                                                                                                                                                                         <  1:131918/248‑1 (MQ=255)
                                                     ggAACGGTACGGGTTCCCAGCCCAAAATCGTTTCCAACC‑AGCTCGTGTTCCTGTTTCCGTACCCCCCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTc                                                                                                                                                                                        <  2:349165/251‑1 (MQ=255)
                                                      gAACGGTCCGGGATTCCAGCTCAAAGACGTTTTCAATC‑ATCTCGAGTTCCTGTTTACGTAACACGCCCCCCAGCGCACCGGCTTCCACTACCACGTAGATGCCCTCAAAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAATATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAAAGTGGT‑CACAACGTAAAGGCAGAAGCGAATCGGGTTGATGATCCGCAAAGCCACAGCGTCt                                                                                                                                                                                       >  1:290552/1‑251 (MQ=255)
                                                           gTACGGGATTCCAGCTCAAAGACGTTTTCAATC‑AGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGACGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGAACCGAAAAGCCACAGCTTCTGgcg                                                                                                                                                                                   >  2:772/1‑250 (MQ=255)
                                                            tACGGGATTCCAGCTCAAAGACGTTTTCAATC‑AGCTCGTGTTCCTGTTTACGTCACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGAACAGGCAGAAGCGCATCGGGTTGATGATACGCCAAGCCACAGCTTCTGgcgc                                                                                                                                                                                  >  1:104529/1‑250 (MQ=255)
                                                                    ccacgccaaacgccGTTTTCAACA‑AGCTCGTGGTCCTTTTTACTGAACACGCCCGCCAGCGCACCGGCTTCCACTGCCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAACCATAcc                                                                                                                                                                         <  2:212060/238‑1 (MQ=255)
                                                                               gACGTTTTCAATC‑AGCTCGTGTTCCTGTTTACGTAACACGACCGCCAGAGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGAACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGAGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGTTGAGACGCAAAGCAACGGCTTCTGGAGCACTTATACCgtggagtttcg                                                                                                                                                              >  2:86494/1‑241 (MQ=255)
                                                                                 ccggttccaaccacctcggtttcgtGTTTACGGAACCCGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGGTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCgg                                                                                                                                                             <  2:19903/227‑1 (MQ=255)
                                                                                       cAATC‑AGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCAACGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTaa                                                                                                                                                       <  1:82795/250‑1 (MQ=255)
                                                                                                   tgtTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACCGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAAc                                                                                                                                             >  1:486596/1‑249 (MQ=255)
                                                                                                       ccTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAATCCACCGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATACGCAAACAGGAt                                                                                                                                        >  2:248004/1‑250 (MQ=255)
                                                                                                         tGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATaa                                                                                                                                      >  1:394726/1‑250 (MQ=255)
                                                                                                           tttACGTAACCCGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAAc                                                                                                                                    <  2:31403/250‑1 (MQ=255)
                                                                                                               cGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGcc                                                                                                                               >  1:127116/1‑251 (MQ=255)
                                                                                                                cagtcacctacgccaggcaaCCGGCTTCCACTACCGCGTAGATGTCATAAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTACGGTCGCCAGGCCGTAGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTGTCGGGGTGAAAGCCGCAAACAGGACAATCATGCCa                                                                                                                              <  2:67322/232‑1 (MQ=255)
                                                                                                                    cacGCCCGCCCGCGCACCGGCTTCCACTACCGCGTCGATGTCATCAGCAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGGGCA‑AACGTAAAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCACTCATACCGAGGCGTTTCGGGGTTAAATCCGCACACAGGATAAACATGCCAGTc                                                                                                                           >  2:239177/1‑250 (MQ=255)
                                                                                                                           gccAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGAGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACg                                                                                                                    >  1:366715/1‑250 (MQ=255)
                                                                                                                            ccAGCGCACCGGCTTCCCCTACCGGGTAGATGTCAGTAGAAGGGATTCCATCTTTACGTACCATTGGCAGTTTGAAAAGACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGAGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACCTGCCAGTCCCTAACgag                                                                                                                  <  2:401895/251‑1 (MQ=255)
                                                                                                                                       gCTTACACTACCACGTAGCTGTCATCAGAAGTGATGTCATCATAACGTACCATTGGCAGTTTAAAAATACGGAAGAATATGTTCGCCAGGCCGTTGAAGAACCACACCAAAGGGGGGAA‑AACGTACAAGCAGAATCGCATCGGGTTGAGGATACGCAAAGACAAAGCTTCTGGCACAACCATAACGATACTTTTCGGGGTTAAATCCGCAAACCGGATAAACATGCCAGTAACTAACgataatgcgagcat                                                                                                       >  2:111660/1‑239 (MQ=255)
                                                                                                                                           ccgcTACCGCGTATATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAg                                                                                                   <  1:386304/248‑1 (MQ=255)
                                                                                                                                                cccGCGTAGATGTCATCAGAAGTGAGTTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAg                                                                                              <  1:109484/250‑1 (MQ=255)
                                                                                                                                                  cgcgTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGtt                                                                                            <  1:415172/251‑1 (MQ=255)
                                                                                                                                                    cgTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCCAACAGGATAAACATGCCAGTCACTACCGCGAAAGAGAGAATAAAGCTCAGTTGc                                                                                          >  1:296282/1‑251 (MQ=255)
                                                                                                                                                     gTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCt                                                                                         <  1:248004/251‑1 (MQ=255)
                                                                                                                                                           gTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGACGATTATGTTCGCCAGGCCGTTGAAGAAACACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAGGCCACAGCTTCTGGCGCAATCACACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATTCCAGTCACTAACGAGAAAGAGAGAATCAAGCTCGGTTGCTCag                                                                                      >  1:496509/1‑248 (MQ=255)
                                                                                                                                                                     gTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACCCCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGa                                                                         <  2:307432/251‑1 (MQ=255)
                                                                                                                                                                     gTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCg                                                                          >  1:21779/1‑250 (MQ=255)
                                                                                                                                                                       gATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTACCGGGAAAGAGAGCATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGAc                                                                        >  2:14813/1‑250 (MQ=255)
                                                                                                                                                                       gATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGAc                                                                        <  1:186988/250‑1 (MQ=255)
                                                                                                                                                                                    cGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCCCTAACGAGAAAGAGAGCATa                                                                                                      >  2:214489/1‑207 (MQ=255)
                                                                                                                                                                                    cGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATa                                                                                                      <  1:214489/207‑1 (MQ=255)
                                                                                                                                                                                        ccATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGgaga                                                           >  2:262271/1‑246 (MQ=255)
                                                                                                                                                                                        ccATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGgaga                                                           <  1:262271/246‑1 (MQ=255)
                                                                                                                                                                                        ccATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAgg                                                      <  1:448809/251‑1 (MQ=255)
                                                                                                                                                                                              gCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTa                                                                                                                       <  1:3441/180‑1 (MQ=255)
                                                                                                                                                                                              gCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTa                                                                                                                       >  2:3441/1‑180 (MQ=255)
                                                                                                                                                                                              gCAGTTTGAAAATACGGAAGATGATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGa                                                <  1:136216/251‑1 (MQ=255)
                                                                                                                                                                                                 gTTTGAAAATACGGAAGATTATGTTCGCAAGGCCGTTGAAGAACCAAAACAACGAGATGCC‑AACGAACAGGAAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGGGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGaaa                                              >  1:299559/1‑250 (MQ=255)
                                                                                                                                                                                                  tttGAAAAAACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGCGCA‑AACTTACAGGCAGCAGCGCATCGGGTTGATGAAACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGACCGTCACTAACGAGAAAGAGAGACTAAAGCTCAGTTGCTCAGAGAGCTCTGCCGCCATGTAGAGGGAGAACAGGCTGTGaaa                                              >  2:259319/1‑249 (MQ=255)
                                                                                                                                                                                                      aaaaTACTGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCGGGCGCAATCAGACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGAGAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCGCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTgg                                        <  2:305965/251‑1 (MQ=255)
                                                                                                                                                                                                       aactccGGAAGATTAGGTTCGCCAGGCCGTGGAAGAACCACCCCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAACCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATATCCATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCGCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGTGAACAGGCTGTGAAAAGCTGga                                       <  2:93762/246‑1 (MQ=255)
                                                                                                                                                                                                              gAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAAAATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGTACGGGCTGTGAAAAGCTGGAGAAAAtt                                >  2:337198/1‑250 (MQ=255)
                                                                                                                                                                                                               aaGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAAAACAGGATGTGAAAAGCTGGAGAAAATGc                               >  2:172418/1‑251 (MQ=255)
                                                                                                                                                                                                               aaGATGATGTTCGCCAGGCCGTAGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGc                               <  2:104469/251‑1 (MQ=255)
                                                                                                                                                                                                                aGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCa                                                                                                                              >  2:105921/1‑155 (MQ=255)
                                                                                                                                                                                                                aGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCa                                                                                                                              <  1:105921/155‑1 (MQ=255)
                                                                                                                                                                                                                     aTGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCAt                          <  1:301567/250‑1 (MQ=255)
                                                                                                                                                                                                                                   gTTGAAGAACCACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCACCGACGATAccgc           >  1:121516/1‑251 (MQ=255)
                                                                                                                                                                                                                                            ccACACCAACGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCACCGACGATACCGCCGCGAATCg  >  1:502389/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
CATCTTCGTGGAGATCAAACCAAATCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATC‑AGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCA‑AACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCACCGACGATACCGCCGAGAATCG  >  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: ‑