Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 597,290 (A)7→8 coding (160/333 nt) cusF → copper/silver export system periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913597,2831.A99.7% 1360.2 / ‑3.2 313coding (153/333 nt)cusFcopper/silver export system periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (148/164);  total (149/164)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 364 total aligned reads displayed.

TAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGCAGACGAAAATGAGTGA  >  NC_000913/597184‑597378
                                                                                                    |                                                                                               
tAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCa                                                                                                <  1:1021552/101‑1 (MQ=255)
tAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCa                                                                                                >  1:1022387/1‑101 (MQ=255)
 aaTGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCaa                                                                                               <  1:1021661/101‑1 (MQ=255)
  aTGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCaaa                                                                                              >  2:1022392/1‑101 (MQ=255)
  aTGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCaaa                                                                                              <  2:1021938/101‑1 (MQ=255)
   tGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCaaaa                                                                                             >  2:1022395/1‑101 (MQ=255)
   tGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCaaaa                                                                                             <  1:1021611/101‑1 (MQ=255)
      ccAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCaaaaaaa                                                                                          >  1:1022397/1‑101 (MQ=255)
        aGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCaaaaaaaat                                                                                        >  2:1022401/1‑99 (MQ=255)
        aGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCaaaaaaaat                                                                                        <  1:1021977/101‑3 (MQ=255)
         ggCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCaaaaaaaatc                                                                                       <  1:1021732/101‑4 (MQ=255)
           cTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCaaaaaaaatcac                                                                                     >  2:1022405/1‑96 (MQ=255)
            tAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCAcc                                                                                     >  2:1022410/1‑101 (MQ=255)
            tAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCAcc                                                                                     >  1:1022409/1‑101 (MQ=255)
            tAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCaaaaaaaatca                                                                                      <  2:1021866/99‑5 (MQ=255)
             aaCGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCAcca                                                                                    >  1:1022417/1‑101 (MQ=255)
             aaCGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCAcca                                                                                    <  2:1021982/101‑1 (MQ=255)
              aCGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCAccat                                                                                   >  2:1022418/1‑101 (MQ=255)
                gAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCAccatcc                                                                                 >  2:1022422/1‑101 (MQ=255)
                gAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCAccatcc                                                                                 >  1:1022424/1‑101 (MQ=255)
                 aaCATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCAccatcca                                                                                <  2:1021578/101‑1 (MQ=255)
                  aCATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCAccatccat                                                                               >  2:1022430/1‑101 (MQ=255)
                       atcatGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGa                                                                          >  1:1022438/1‑101 (MQ=255)
                       atcatGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGa                                                                          >  2:1022435/1‑101 (MQ=255)
                        tcatGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGAt                                                                         <  2:1022003/101‑1 (MQ=255)
                        tcatGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGAt                                                                         <  2:1021796/101‑1 (MQ=255)
                        tcatGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGAt                                                                         <  1:1021899/101‑1 (MQ=255)
                          atGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATcc                                                                       >  1:1022443/1‑101 (MQ=255)
                           tGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCg                                                                      <  2:1021960/101‑1 (MQ=255)
                             aaaCCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGAt                                                                         >  2:1022447/1‑96 (MQ=255)
                             aaaCCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGAt                                                                    >  1:1022446/1‑101 (MQ=255)
                             aaaCCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGAt                                                                    <  1:1021888/101‑1 (MQ=255)
                                ccAAGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATtgc                                                                 >  1:1022452/1‑101 (MQ=255)
                                 cATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCAccatc                                                                                  <  1:1021766/83‑1 (MQ=255)
                                 cATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATtgct                                                                <  2:1022162/101‑1 (MQ=255)
                                 cATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATtgct                                                                >  2:1022455/1‑101 (MQ=255)
                                  aTGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATtgctg                                                               <  1:1021366/101‑1 (MQ=255)
                                     aGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCg                                                            >  2:1022462/1‑101 (MQ=255)
                                     aGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCg                                                            >  1:1022463/1‑101 (MQ=255)
                                      gCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCTATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGt                                                           <  2:1021691/101‑1 (MQ=255)
                                      gCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGt                                                           <  1:1021708/101‑1 (MQ=255)
                                        gAAGCACGACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGa                                                         >  2:1022468/1‑101 (MQ=255)
                                         aaGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGaa                                                        <  1:1021810/101‑1 (MQ=255)
                                         aaGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGaa                                                        <  1:1022286/101‑1 (MQ=255)
                                           gCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACt                                                      <  2:1021780/101‑1 (MQ=255)
                                            cacaACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGacc                                                       <  1:1021748/99‑3 (MQ=255)
                                             acaACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTgg                                                    <  2:1021851/101‑1 (MQ=255)
                                             acaACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTgg                                                    >  1:1022473/1‑101 (MQ=255)
                                             acaACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTgg                                                    >  2:1022472/1‑101 (MQ=255)
                                              caACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAAc                                                       <  2:1021969/97‑1 (MQ=255)
                                               aaCCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGcc                                                  <  2:1021658/101‑1 (MQ=255)
                                               aaCCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGcc                                                  >  2:1022477/1‑101 (MQ=255)
                                               aaCCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGcc                                                  <  2:1021763/101‑1 (MQ=255)
                                                 ccACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCgg                                                <  1:1021840/101‑1 (MQ=255)
                                                   acaGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGgag                                              >  2:1022486/1‑101 (MQ=255)
                                                   acaGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGgag                                              >  1:1022483/1‑101 (MQ=255)
                                                     aGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGcc                                                  <  2:1021818/95‑1 (MQ=255)
                                                     aGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGAt                                            <  1:1021832/101‑1 (MQ=255)
                                                      ggTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATg                                           <  1:1021802/101‑1 (MQ=255)
                                                        ttattaGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCTATTGCTGCCGTGAACTGGCCGGAGATGAc                                         >  2:1022493/1‑101 (MQ=255)
                                                        ttattaGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGAc                                         >  1:1022492/1‑101 (MQ=255)
                                                        ttattaGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGAc                                         <  2:1021725/101‑1 (MQ=255)
                                                         tattaGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGAcc                                        >  1:1022496/1‑101 (MQ=255)
                                                           ttaGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCAt                                      <  2:1021792/101‑1 (MQ=255)
                                                            taGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATg                                     >  1:1022504/1‑101 (MQ=255)
                                                            taGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATg                                     >  2:1022502/1‑101 (MQ=255)
                                                             aGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATgc                                    <  2:1021918/101‑1 (MQ=255)
                                                                gcCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCCCCACCCATCACGATCCGATTGCTGCCGTGAACTGGCCgg                                                <  2:1022021/86‑1 (MQ=255)
                                                                gcCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATg                                           >  2:1022508/1‑91 (MQ=255)
                                                                  cACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCttt                               <  2:1021890/101‑1 (MQ=255)
                                                                  cACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCttt                               >  1:1022511/1‑101 (MQ=255)
                                                                   aCTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTa                              <  1:1021998/101‑1 (MQ=255)
                                                                     tGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTAcc                            <  1:1022063/101‑1 (MQ=255)
                                                                     tGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTAcc                            >  2:1022514/1‑101 (MQ=255)
                                                                      ggCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCa                           <  1:1021861/101‑1 (MQ=255)
                                                                       gCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCCCCACCCATCACGATCCGATTGCTGCCGTGAACTGGcc                                                  <  2:1021423/77‑1 (MQ=255)
                                                                        cGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATc                         >  1:1022518/1‑101 (MQ=255)
                                                                          tGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCAc                       >  2:1022525/1‑101 (MQ=255)
                                                                          tGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCAc                       >  1:1022524/1‑101 (MQ=255)
                                                                             tAAAGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCAcccc                    <  1:1021933/101‑1 (MQ=255)
                                                                               aaGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGc                  <  2:1022340/101‑1 (MQ=255)
                                                                               aaGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGc                  >  2:1022532/1‑101 (MQ=255)
                                                                               aaGGGTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGc                  >  1:1022531/1‑101 (MQ=255)
                                                                                 gggTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGCAg                <  1:1022039/101‑1 (MQ=255)
                                                                                   gTATCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGCAGAc              <  1:1022151/101‑1 (MQ=255)
                                                                                     aTCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGCAGACGa            >  2:1022539/1‑101 (MQ=255)
                                                                                     aTCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGCAGACGa            >  1:1022540/1‑101 (MQ=255)
                                                                                      tCGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGCAGACGaa           <  2:1021904/101‑1 (MQ=255)
                                                                                       cGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGCAGACGaaa          >  2:1022543/1‑101 (MQ=255)
                                                                                       cGATCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGCAGACGaaa          <  2:1021874/101‑1 (MQ=255)
                                                                                         aTCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGCAGACGAAAAt        >  2:1022551/1‑101 (MQ=255)
                                                                                         aTCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGCAGACGAAAAt        >  1:1022548/1‑101 (MQ=255)
                                                                                          tCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGt                                                           <  2:1021698/49‑1 (MQ=39)
                                                                                          tCTGGAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGCAGACGAAAATg       >  1:1022555/1‑101 (MQ=255)
                                                                                             ggAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGCAGACGAAAATGAGt    >  1:1022559/1‑101 (MQ=255)
                                                                                             ggAAAGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGCAGACGAAAATGAGt    <  1:1021971/101‑1 (MQ=255)
                                                                                             ggAAAGCAAAAAAAACCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGCAGACGAAAATGAGt    <  2:1022089/101‑1 (MQ=255)
                                                                                               aaaGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGCAGACGAAAATGAGTGa  <  2:1021912/101‑1 (MQ=255)
                                                                                               aaaGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGCAGACGAAAATGAGTGa  <  1:1022110/101‑1 (MQ=255)
                                                                                               aaaGCAAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGCAGACGAAAATGAGTGa  <  2:1022237/101‑1 (MQ=255)
                                                                                                    |                                                                                               
TAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCTTTACCATCACCCCGCAGACGAAAATGAGTGA  >  NC_000913/597184‑597378

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: