Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,358,859 T→C Y110C (TAT→TGT)  puuP ← putrescine:H(+) symporter PuuP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,358,8590TC95.2% 122.6 / NA 43Y110C (TAT→TGT) puuPputrescine:H(+) symporter PuuP
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base C (13/27);  minor base A (2/0);  total (15/28)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

GACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATCA  >  NC_000913/1358615‑1359105
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
gACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCca                                                                                                                                                                                                                                                  <  1:279655/251‑1 (MQ=255)
      gCCAACGCCTTCTACTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCagcagc                                                                                                                                                                                                                                             <  1:234907/250‑1 (MQ=255)
          aCGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTAACACCCACAACC‑ACGGCGGCACTTCCGGGGAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGAAA‑ATCCAGCAGCGATg                                                                                                                                                                                                                                         >  1:45217/1‑250 (MQ=255)
           cGCCTTCTCCTGTATGCATTCCCTGAACTACCAGGAAGATAAACACCCCCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGGTGACGCTCTTCAGATTCGCGGCGGTTAAAAGGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCTCTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCACCGGCAAAAAGAGACA‑ATCCagcagctattac                                                                                                                                                                                                                                       <  2:90958/251‑7 (MQ=255)
            gCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGAcc                                                                                                                                                                                                                                      <  2:143229/251‑1 (MQ=255)
                       tATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCAt                                                                                                                                                                                                                            <  1:181451/250‑1 (MQ=255)
                                  tGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGG‑GGAGAGATAGATTTTTGCCAACACGACGGTAATCATCGGCCAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATCAAACCGActgg                                                                                                                                                                                                                >  1:409791/1‑248 (MQ=255)
                                  tGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGAATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCCCCCACACCC‑CCGGCGGCACTTCCGGGGAGAGGGG‑GGGGAGATAGATTTTCTCCAAAAAGACGGTAATTATCGGAAAAAAGAGACA‑AAACAGCAACGAAGACCAAGCGACCAAAAATCCCACGTg                                                                                                                                                                                                                >  2:477010/1‑251 (MQ=255)
                                        aCCAGGAAGATAAACACCACCATGATGGAGAGTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCAGCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGt                                                                                                                                                                                                           <  2:82226/250‑1 (MQ=255)
                                             gAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCg                                                                                                                                                                                                     <  2:146015/251‑1 (MQ=255)
                                                gATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGGAGGTTTCCACCCACACCC‑CCGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCACCAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGAAAAGCCGACCATAAATCCCACGTGCGGGGTAATCGACt                                                                                                                                                                                                  >  1:114824/1‑251 (MQ=255)
                                                    aaCACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGAACTTCCCGGCAGGGGGC‑GGAGAGATAGATTTTCCCCAACCAGACGGTAATCATCGGCAAAAAGAGACA‑ATCCAGCCGCGATGACCCGCCGACCCTAAATCCGACGTGCGGGTTAATCGCCttttt                                                                                                                                                                                              >  1:453335/1‑250 (MQ=255)
                                                      caccacCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTgcg                                                                                                                                                                                            <  2:39597/251‑1 (MQ=255)
                                                            caTGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAgg                                                                                                                                                                                      <  1:431139/251‑1 (MQ=255)
                                                                  ggAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑GCGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAAGCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGTGTTAATCGACTGTTGCGCGTAGGTATAgg                                                                                                                                                                                <  2:251319/251‑1 (MQ=255)
                                                                   gagaTTTGCACCAGTACAAACAGGGTATTGAAGATAGCGACCAGGTTGACGCTCTGCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGTGCGTAGGTCTAGGc                                                                                                                                                                               <  2:91032/251‑1 (MQ=255)
                                                                    agaTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCAACCACAACC‑ACGGCGGCACTTCCGGGGAGAGGGG‑GGAGAGAAAGAATTTCGCCAACAAGACGTTAATCATCGGCAAAAAAAGACA‑ATCCAGCAGCGGTGAACAGCCGAACATAAATCCGACGGGCGGGTTAATGGACATTTGCGCGTAGGTATAggac                                                                                                                                                                              >  1:292202/1‑249 (MQ=255)
                                                                        ttGCACCAGTACAAACAAGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTGGCGCGTAGGTATAGGCCGaa                                                                                                                                                                           <  2:326193/250‑1 (MQ=255)
                                                                           cACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTAGAGGCCGAACCgg                                                                                                                                                                       <  1:404757/251‑1 (MQ=255)
                                                                                      aCAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGATAAAATGGCGACGAAGGTTACAACCCACACCCCACG‑CCGCACTTCCGGGAAGAGGGCCGGA‑AGATAGATTTTTGCCAACAAGACGGTAAAAAATGGCAAAAAAAAACA‑ATCCAGCAGCGATGACCAGCCCACCCAAAATCCGACGTGCGGGTTAATCGCCTTTTTCACGTAGGTATAGGCCGAAcccggctccagtaa                                                                                                                                                             >  1:73215/1‑238 (MQ=255)
                                                                                       cAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTAACACCCCCACCC‑AAGGCGGCACTTCCGGGACGAGGGC‑GGAGAGATAGATTTTTGCCAACAAGACGCTAATCATCGGCACAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGGCCCTAAATCCGACGTGCGGGGTAATCGGCTTTTGCGCGTAGGTATAAGCCGAACCGGCCTCCGGGAAcc                                                                                                                                                           >  1:34657/1‑250 (MQ=255)
                                                                                         gggTAATGAAGTGAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTgg                                                                                                                                                         <  2:453335/251‑1 (MQ=255)
                                                                                          ggTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTAACACCCCCACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGGTTTTCGCCAACAAGACGTTAAACCTCGGCAAAAAGAGACA‑AACCAGCAGCGATGACCAGCCGACCATAAAACCGACGTGGGGGGTAATCGACTTTTTCGGGTAGGAATAGGACGAACCGGCCCCCCGAAAATGGc                                                                                                                                                        >  2:379851/1‑251 (MQ=255)
                                                                                                     tAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGttt                                                                                                                                              <  2:45217/250‑1 (MQ=255)
                                                                                                      aGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGc                                                                                                                                            <  2:38701/251‑1 (MQ=255)
                                                                                                       gCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGcc                                                                                                                                           <  2:9630/251‑1 (MQ=255)
                                                                                                          aCCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCGCTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTa                                                                                                                                        <  1:166614/251‑1 (MQ=255)
                                                                                                          aCCAGGTTGACGCTCGTCAGATTCGCGGCGGTTAAAAGGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTa                                                                                                                                        <  2:171925/251‑1 (MQ=255)
                                                                                                                   accagctTCACATTCGCGTCGGTCAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GTAGAGATAGATTGTCGCCAACGAGACCTTAATCATCGGCCAAAAGAGACC‑ATCCAGCAGCGATGCCCAGCCGATCATCGATCCGCCGTGCGGGTTAATCGACGTTTGCGCGTAGGTATAGGCCTAACCGGCGTCCGGAACCGGGCTCACCAGTTTTCCGTAGCTGATacc                                                                                                                               <  2:50524/246‑3 (MQ=255)
                                                                                                                       tctTCGGATTCGCGGCGGTTAAAATGGAGACGAAGGTTACCACCCACACAC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTTGCCCAAAAGACGGTAATCATCGGGAAAAAGAGAAA‑ATc                                                                                                                                                                                                                                                    >  2:409905/1‑130 (MQ=255)
                                                                                                                       tctTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGTCAAAAAGAGACA‑ATc                                                                                                                                                                                                                                                    <  1:409905/130‑1 (MQ=255)
                                                                                                                                                                                 gcACGTCCGGGAAGGGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAACTAGGACGCCGGAACGTGACCGTCGCTAATg                                                                  >  2:415016/1‑251 (MQ=255)
                                                                                                                                                                                   ccTTCCGGGAAGAGGGC‑GGAGAGATAGATTTCCCCCAACAATCCGTTAATCATCGGCAAAAAGAGACA‑ACCCAGCAGCTATTACCAGCCGACCATAAACCCGCCGTGCGGGTTAATCGACTGTTGCGCGTAGGTATAGGCCGAACCGCCCGCCGGAAACGGGCAAACCAGTTGGCCGTTGCTGATAGCGGTAAACGGTACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGcc                                                                <  2:34657/250‑1 (MQ=255)
                                                                                                                                                                                      tCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCt                                                                      >  1:179808/1‑242 (MQ=255)
                                                                                                                                                                                      tCCGGAAAGAGGGG‑GGAGAGATAGATTTTCGCAAACAAGACGTTAACCAGCGGCAAAAAGAGACA‑AACCAGCAGCGAGGACCAGCCGACCAAAAATCCGACGGGCGGGTTAACCGTCTTTTGCGCGTAGGTATAGGCCGAACGGCCCCCCGGAATCTGGCGACCCAGTTGGCCATAGTTGATAGCGGTAAACAGCACCCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCt                                                                      <  2:179808/242‑1 (MQ=255)
                                                                                                                                                                                       ccGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGACCGTGACCGTCGCTAATGCCgg                                                              >  2:382450/1‑249 (MQ=255)
                                                                                                                                                                                        cGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGCCGTGCGGGTTACTCGACTTTTGCGCGTAGGTATAGGCCGACCCGACCTCCTGAAACTGGCGCACCAGTTTGCCGTAGCTGATATCGGAAAACAGCACACCCACCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTCATGCCgttg                                                            >  1:90457/1‑247 (MQ=255)
                                                                                                                                                                                                         atagatatTTGCGCAAACAAGACTGTAAACAGCGGCAAAAAGAG‑TATGTCCAGCAGCGAGGACCTGCGGATCAAAAATCAGGCGCGCGGGGTAATCGATTTGTGCGTGTGGGTGTAGGCGGAACCGTTCTCGGGAAACTGGCAAAGGAGTTTGCGGGTGCGGATCGCGGTAATCAGGGCGCCCGCCATCGCCAGCAAATAGGACGCCGGAGCGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATc                                           <  2:455096/245‑1 (MQ=255)
                                                                                                                                                                                                           atgaattgtCGCAAACAAGACGTAAATCATCGAAAAAAGGAGACA‑ATCCAGCAGCGATGACCAGCGGACCATAAATCCGACGTCCGGGTAAATCGCCTGTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCAGTAGCTGAGCGCGGAAAACAGCACACCCGCCAGCGCCAGCAAATAGAACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCaa                                         <  2:211521/243‑1 (MQ=255)
                                                                                                                                                                                                                    cGCCAACAAGACGTGAATCAACGGCAAAAAGAGACA‑ATCCACCACGGATGCCCAGCCGACCATAAATCCGACGTGCGGTTTAATCTACTTTTGCGCGTAGGTATAGGCCGAACCGGCCCCCGGAAACTGGCGACCCAGTTTGCCGTAGCTGATAGCGTTAAACACCACCCCCGCCAGCTCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCa                                <  2:49372/251‑1 (MQ=255)
                                                                                                                                                                                                                          aaaGACTTTAATCATCGGCAAAAAGAGCCA‑ATCCAGCAGCGATGACCAGCCGACCATAAACCCGGCGTGCGGGTTAAGCGACCTTTGCGCGGGGGGATAGGCCGAACCGGCCCCCGGAACCTGGCGAACCAGTTTGCCGTAGGTGATGGCGGTAAACAGCACTCCCGCCAGCGCCAGCAAATAGTACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCg                          <  1:232448/250‑1 (MQ=255)
                                                                                                                                                                                                                                   aaTCATCGGCAAAAAGAGACA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCCGCAAATAGGACGCCGGCACGTCACCGTCGCTACTGCCGGAGACAATACCAAAAGTATCAAATAAGGACATCGGCGTGAGATAgg                 >  1:59495/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  gagaCA‑ATCCAGCAGCGATGACCGGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTATCGGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCAtcatca  <  1:359707/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
GACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCC‑ACGGCGGCACTTCCGGGAAGAGGGC‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATA‑ATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATCA  >  NC_000913/1358615‑1359105

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: ‑