Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 717,518 G→A 5.8% P250S (CCA→TCA)  potE ← putrescine/proton symporter: putrescine/ornithine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913717,5180GA 5.8% 218.0 / 2.3 71P250S (CCA→TCA) potEputrescine/proton symporter: putrescine/ornithine antiporter
Reads supporting (aligned to +/- strand):  ref base G (33/33);  new base A (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

ACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTCCACTACATCAGTATTCGCACACGCAGACTCCAGACCAAGA  >  NC_000913/717370‑717657
                                                                                                                                                    |                                                                                                                                           
acCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTgg                                                                                                                                             <  2:59742/149‑1 (MQ=255)
        aaTGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGAGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATcc                                                                                                                                     >  8:360780/1‑149 (MQ=255)
          tGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCg                                                                                                                                   <  8:209939/149‑1 (MQ=255)
           ggTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCgg                                                                                                                                  >  1:398677/1‑149 (MQ=255)
            gTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCgg                                                                                                                                  <  4:256516/148‑1 (MQ=255)
                aCTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATc                                                                                                                             <  3:160169/149‑1 (MQ=255)
                        ccaAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTtggt                                                                                                                     <  2:11635/149‑1 (MQ=255)
                        ccaAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTtggt                                                                                                                     >  2:292202/1‑149 (MQ=255)
                         caAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTtggtg                                                                                                                    <  8:418689/149‑1 (MQ=255)
                           aGTAGCGAACAGGATCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGAGTGAACATCTGCGCGAAGGACAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTTGTGga                                                                                                                  <  3:8394/149‑1 (MQ=255)
                               gCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGAc                                                                                                               <  4:199274/148‑1 (MQ=255)
                                     cGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATaa                                                                                                        <  6:402689/149‑1 (MQ=255)
                                      gcagcaGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCAGAATCCCGGCaa                                                                                                                               >  6:185827/1‑125 (MQ=255)
                                      gcagcaGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCaa                                                                                                                               <  5:185827/125‑1 (MQ=255)
                                      gcagcaGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCaa                                                                                                                               <  5:264371/125‑1 (MQ=255)
                                      gcagcaGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCaa                                                                                                                               >  6:264371/1‑125 (MQ=255)
                                           agGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCAcc                                                                                                  >  3:55671/1‑149 (MQ=255)
                                             gACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATccc                                                                                                                                    >  6:301374/1‑113 (MQ=255)
                                             gACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATccc                                                                                                                                    <  5:301374/113‑1 (MQ=255)
                                                aTCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACcgccgc                                                                                             <  4:55671/149‑1 (MQ=255)
                                                   aCCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCAcc                                                                                          <  7:309071/149‑1 (MQ=255)
                                                      aTCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTaa                                                                                       >  7:69892/1‑149 (MQ=255)
                                                      aTCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGAAATAAATAACCGCCGCACCTac                                                                                       >  1:399827/1‑148 (MQ=255)
                                                            gcCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTCGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTAcc                                                                                 >  4:70918/1‑149 (MQ=255)
                                                                tAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCa                                                                                                    <  5:397177/126‑1 (MQ=255)
                                                                tAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCa                                                                                                    >  6:397177/1‑126 (MQ=255)
                                                                  aTGACTTTACCCACTTCCGGCGGGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGa                                                                            <  1:292202/148‑1 (MQ=255)
                                                                  aTGACTTTACCCACTTCCGGAGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAg                                                                           >  7:346208/1‑149 (MQ=255)
                                                                  aTGACTTCACCCACTTCCGGCGTTAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGat                                                                           <  1:430605/149‑2 (MQ=255)
                                                                    gACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTa                                                                         >  3:49590/1‑149 (MQ=255)
                                                                           cccACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGa                                                                   <  8:69892/148‑1 (MQ=255)
                                                                           cccACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGAt                                                                  >  1:292189/1‑149 (MQ=255)
                                                                            ccACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGAtt                                                                 <  7:297816/149‑1 (MQ=255)
                                                                            ccACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGAtt                                                                 >  6:72359/1‑149 (MQ=255)
                                                                             cACTTCCGGCGTGAACATCTGCGCGAAGGCCTGAACAATTGGTGCGGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGTGACGATATATATCACCGCCGCACCTAACGTAGCGCCGAGTACCGCGATTg                                                                <  8:346208/149‑1 (MQ=255)
                                                                             cACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTg                                                                <  3:70918/149‑1 (MQ=255)
                                                                             cACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTg                                                                >  6:389595/1‑149 (MQ=255)
                                                                                  ccGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGcac                                                            <  2:292189/148‑1 (MQ=255)
                                                                                  ccGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGca                                                             <  5:72359/147‑1 (MQ=255)
                                                                                      cGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTAc                                                       >  1:311876/1‑149 (MQ=255)
                                                                                       gTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGACACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACg                                                      >  3:3864/1‑149 (MQ=255)
                                                                                            cttCTGCGCGATGGCCACACCAAATGTTGCCGTGGAATTTGCCAGCTCCATTTTTGCCACAATCCCGGAAATAACGTTGGTGTATAGTATATAAATCAACGCCGCACCT‑ACGTACAGCCGAGTACCGCGATTGGCACATTACGTTCCgg                                                <  2:399827/147‑1 (MQ=255)
                                                                                             aTCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCgg                                                <  2:398677/149‑1 (MQ=255)
                                                                                                 gcgcGAAGGCCAGACAAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGttt                                            >  8:359944/1‑149 (MQ=255)
                                                                                                   gcgAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTc                                          >  1:21838/1‑149 (MQ=255)
                                                                                                   gcgAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTc                                          <  7:360780/149‑1 (MQ=255)
                                                                                                     gAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGa                                                                            >  3:195544/1‑113 (MQ=255)
                                                                                                     gAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGa                                                                            <  4:195544/113‑1 (MQ=255)
                                                                                                     gAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTCCa                                        <  1:404822/149‑1 (MQ=255)
                                                                                                         gCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTCCACTAc                                    >  6:384205/1‑149 (MQ=255)
                                                                                                           cAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAGTCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTCCACTACAt                                  >  4:404963/1‑149 (MQ=255)
                                                                                                           cAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTCCACTACAt                                  <  2:311876/149‑1 (MQ=255)
                                                                                                                cAAATGGTGCCGTTGAATTTGCCCGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCTCCGAGTACCGCGATTGGCACATGACGTTCCGGG‑TTTCCAATACATAAGTAt                            >  6:340527/1‑149 (MQ=255)
                                                                                                                cAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTCCACTACATCAGTa                             >  2:81144/1‑149 (MQ=255)
                                                                                                                    tGGTGCCGTTGAATTTGCCAGCTCCATATTTGACACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCgg                                                >  8:335456/1‑126 (MQ=255)
                                                                                                                    tGGTGCCGTTGAATTTGCCAGCTCCATATTTGACACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCgg                                                <  7:335456/126‑1 (MQ=255)
                                                                                                                      gTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTCCACTACATCAGTATTCGc                        <  7:359944/148‑1 (MQ=255)
                                                                                                                         ccGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTCCACTACATCAGTATTCGCacac                    >  3:341486/1‑149 (MQ=255)
                                                                                                                           gTTGAATTTGCCAGCTCCATATTTGGCGCAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTCCACTACATCAGTATTCGCACACg                   >  2:285184/1‑148 (MQ=255)
                                                                                                                           gTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTCCACTACATCAGTATTCGCACACGc                  <  5:437157/149‑1 (MQ=255)
                                                                                                                            ttGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTCCACTACATCAGTATTCGCACACGCa                 <  4:341486/149‑1 (MQ=255)
                                                                                                                                      cAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTCCACTACATCAGTATTCGCACACGCAGACTCCAGAc       <  1:81144/149‑1 (MQ=255)
                                                                                                                                      cAGCTCCATATTTGACACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTCCACTACATCAGTATTCGCACACGCAGACTCCAGAc       <  4:3864/149‑1 (MQ=255)
                                                                                                                                       aGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTc                                                   >  7:379988/1‑104 (MQ=255)
                                                                                                                                       aGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTc                                                   <  8:379988/104‑1 (MQ=255)
                                                                                                                                       aGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTCCACTACATCAGTATTCGCACACGCAGACTCCAGAc       >  7:219418/1‑148 (MQ=255)
                                                                                                                                       aGCTCCATATTTGGAACAATCCCGGCAATAACGTTGGTGGAGAAGATATAAATCACCGCCGCACCTAACGTACAGCCGAGTACAGCGATTGGCACATTACTTTCCGGGTTTTCCACTACATCAGTATTCGCACAAGCAGACTCAAGAc       >  6:320989/1‑148 (MQ=255)
                                                                                                                                         cTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGAACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTCCACTACATCAGTATTCGCACACGCAGACTCCAGAcc      >  7:49334/1‑147 (MQ=255)
                                                                                                                                           ccATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTCCACTACATCAGTATTCGCACACGCAGACTCCAGACCAAg   >  3:387947/1‑148 (MQ=255)
                                                                                                                                            cATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTCCACTACATCAGTATTCGCACACGCAGACTCCAGACCAAGa  >  6:376017/1‑148 (MQ=255)
                                                                                                                                            cATATTTGGCACAATCCCGGCAATCACGATGGTGGACACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGCTTATCCAATACATCAGTAt                            >  6:151540/1‑122 (MQ=255)
                                                                                                                                            cATATTTGGCACAATCCAGGCAATCACGTTGGTGGAGACGATATAAATCACCGCGGCACCTAACGTTCCGCCGAGTAGCGCGATTGGCACATTAGGTTCCTGGTTTTCCACTACATCAGTAt                            <  5:151540/122‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
ACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCATATTTGGCACAATCCCGGCAATCACGTTGGTGGAGACGATATAAATCACCGCCGCACCTAACGTACCGCCGAGTACCGCGATTGGCACATTACGTTCCGGGTTTTCCACTACATCAGTATTCGCACACGCAGACTCCAGACCAAGA  >  NC_000913/717370‑717657

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