Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 595,520 T→A 16.3% intergenic (‑77/‑80) cusR ← / → cusC response regulator in two‑component regulatory system with CusS/copper/silver efflux system, outer membrane component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913595,5200TA16.3% 330.5 / 52.1 140intergenic (‑77/‑80)cusR/cusCresponse regulator in two‑component regulatory system with CusS/copper/silver efflux system, outer membrane component
Reads supporting (aligned to +/- strand):  ref base T (63/53);  new base A (13/10);  total (76/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.58e-01

Only 100 of 143 total aligned reads displayed.

CTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCT  >  NC_000913/595422‑595619
                                                                                                  |                                                                                                   
cTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                    <  1:166603/100‑1 (MQ=255)
  tCGACTATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAGTCGCTTATTGGCAAAATGACAATTTTg                                                                                                  <  1:4064515/100‑1 (MQ=255)
    gACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTc                                                                                                >  1:3546452/1‑100 (MQ=255)
    gACAATCAACAGATTCATATTTCCTCCGCCTGTTGCTCGGGCAACTCTAGAGTAGCGGGAACCGATGGCAATCGCTTATTCCCCAACCGACAATTTTgcc                                                                                                >  2:1857223/1‑98 (MQ=255)
      cAATCAACAGTTTCATCTTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCa                                                                                               <  2:1484199/99‑1 (MQ=255)
      cAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAt                                                                                              >  2:4130914/1‑100 (MQ=255)
       aaTCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAtt                                                                                             >  2:1113343/1‑100 (MQ=255)
        aTCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGCTCAGATGGCAATCGCTTATTGGCCAAATGACAATTTTGTCAttt                                                                                            >  1:1737791/1‑100 (MQ=255)
        aTCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTAAttt                                                                                            >  1:851937/1‑100 (MQ=255)
         tCAACAGTTTCATAGTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttt                                                                                            <  2:1952881/99‑1 (MQ=255)
          cAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAttttt                                                                                          >  2:412346/1‑100 (MQ=255)
           aaCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTc                                                                                         >  1:1731204/1‑100 (MQ=255)
           aaCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgc                                                                                         <  2:1198712/100‑3 (MQ=255)
           aaCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgc                                                                                         <  2:2170128/100‑3 (MQ=255)
            aCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCt                                                                                        >  2:164944/1‑100 (MQ=255)
            aCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgcc                                                                                        <  2:36088/100‑4 (MQ=255)
            aCAGTTTCATATTGCCTCCGCATGTTGCCCGGGCAATTCTAGAATAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgcc                                                                                        <  2:3615929/100‑4 (MQ=255)
             cAGTTTCATATTTCCTCCGCATGTTGCCCTGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCt                                                                                        <  2:3021595/99‑1 (MQ=255)
             cAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTg                                                                                       <  2:3725432/100‑1 (MQ=255)
             cAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgcca                                                                                       <  2:3894754/100‑5 (MQ=39)
              aatttCATGTTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAtttt                                                                                           <  2:1990212/93‑1 (MQ=255)
              aGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGTCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTg                                                                                       <  1:3871027/99‑1 (MQ=255)
                 ttCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGCTCAGATGGCCATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCa                                                                                    >  2:2182120/1‑99 (MQ=255)
                  tcacatTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCAcn                                                                                  <  1:3189031/96‑2 (MQ=255)
                  tCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTCTTGGCAGAATGACAATTTTGTCATTTTTCTGTCAc                                                                                   >  2:2529205/1‑99 (MQ=255)
                  tCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCAcc                                                                                  <  1:2160864/100‑1 (MQ=255)
                   cATATTTCCTCCGCATGTTGCCCGGGCAATTCTCTGGTTGCGGGATTAGATGGCAATCGCTTATTGGCAAAATGACGATTTTGTCATTTTGCTGTCACCg                                                                                 <  2:1743168/100‑1 (MQ=255)
                   cATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGATTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgccaataagc                                                                                 <  2:1967579/100‑11 (MQ=37)
                   cATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCg                                                                                 <  1:1208556/100‑1 (MQ=255)
                        ttCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaaaa                                                                            >  1:1693512/1‑100 (MQ=255)
                        ttCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaaaa                                                                            <  2:3724272/100‑1 (MQ=255)
                        ttCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaaa                                                                             <  2:3206922/99‑1 (MQ=255)
                         ncctccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAAt                                                                           >  1:3047161/2‑100 (MQ=255)
                          cctgcGCATGTTGCCCGGGATATTCTAGAGTAGCGGATCAGCTGGCAATCGCTTAGTGGCAAAACGACAATTTTGTCATTTTTCTGTCACCGGAAAATCa                                                                         <  1:2351893/100‑1 (MQ=255)
                          cctccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATc                                                                          <  1:164944/100‑1 (MQ=255)
                              cGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCaga                                                                       >  2:2357077/1‑99 (MQ=255)
                                cATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGcc                                                                    <  1:2251335/100‑1 (MQ=255)
                                  tGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTg                                                                  >  1:3450862/1‑100 (MQ=255)
                                  tGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTg                                                                  <  1:334485/100‑1 (MQ=255)
                                   gTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTgg                                                                 >  2:1083058/1‑100 (MQ=255)
                                     tGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTTACCGGAAAATCAGAGCCTGGCg                                                               >  2:1411096/1‑100 (MQ=255)
                                     tGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCg                                                               <  2:1847847/100‑1 (MQ=255)
                                      gCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGa                                                              <  1:3853814/100‑1 (MQ=255)
                                            gCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAg                                                        <  1:224348/100‑1 (MQ=255)
                                             cAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGt                                                       <  1:2233893/100‑1 (MQ=255)
                                             cAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGt                                                       >  2:2863872/1‑100 (MQ=255)
                                                 tCTAGAGTAGCGGGATCAGATGGCTATCGCTTATTGGCAAAATGACAATTTTGTCAGTTGTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggc                                                   <  2:1014859/100‑1 (MQ=255)
                                                  cTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcg                                                  >  1:553515/1‑100 (MQ=255)
                                                    agagTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcggc                                                <  1:994175/100‑1 (MQ=255)
                                                     gagTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCa                                               <  1:4063640/100‑1 (MQ=255)
                                                       gTAGCTGGATCAGATGGCGATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCa                                               <  2:2493499/98‑1 (MQ=255)
                                                       gTAGCGGGAGGAGATGGCAATCGCCTATTGGCAAAACGACAATTTTGTCAAGTTTATGTCACCGGAAAATCAGAACCTCGCGTGTAAAGTTGGCGGCAt                                              >  2:73982/1‑99 (MQ=255)
                                                         aGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATaaa                                           >  1:3676803/1‑100 (MQ=255)
                                                            gggATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAAt                                         >  1:1903319/1‑99 (MQ=255)
                                                             ggATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCa                                       >  2:2804912/1‑100 (MQ=255)
                                                              gATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCAc                                      <  1:3064604/100‑1 (MQ=255)
                                                              gATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCACAAAATCAc                                      <  2:2855064/100‑1 (MQ=255)
                                                               aTCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCAcc                                     <  1:3746015/100‑1 (MQ=255)
                                                                tCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCa                                    >  1:3199579/1‑100 (MQ=255)
                                                                  aGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGa                                  >  2:3137516/1‑100 (MQ=255)
                                                                  aGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGa                                  <  2:1363135/100‑1 (MQ=255)
                                                                  aGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGa                                  <  1:1968158/100‑1 (MQ=255)
                                                                   gATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGaa                                 >  2:3908380/1‑100 (MQ=255)
                                                                   gATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGaa                                 >  2:2583943/1‑100 (MQ=255)
                                                                   gATGGCAATCGCTTATTAGCAAAATGACAATTTAGTCATTTTTCTGTCACGGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCaatag                                   >  2:275946/1‑94 (MQ=255)
                                                                    aTGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGaaa                                >  1:2580246/1‑100 (MQ=255)
                                                                        cAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTAt                            <  2:755083/100‑1 (MQ=255)
                                                                            cGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGc                        >  2:1487609/1‑100 (MQ=255)
                                                                            cGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGc                        >  1:3350248/1‑100 (MQ=255)
                                                                             gCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGcc                       <  2:2362790/100‑1 (MQ=255)
                                                                             gCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGcc                       >  2:1235029/1‑100 (MQ=255)
                                                                             gCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGcc                       >  1:350831/1‑100 (MQ=255)
                                                                               gtgacagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTa                     >  1:1091062/9‑100 (MQ=38)
                                                                                 gacagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATg                   >  2:375722/7‑100 (MQ=38)
                                                                                 gacagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATg                   >  1:2906664/7‑100 (MQ=38)
                                                                                  ttGGCAAAATGACAATTTTGGCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGt                  >  2:2570808/1‑100 (MQ=255)
                                                                                  acagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGATAATGGTGAGTAAAGTTGACGGCATAAAATCACCAGAAATTATGAGCCTATGt                  >  2:1863997/6‑100 (MQ=37)
                                                                                  acagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGt                  >  2:3316921/6‑100 (MQ=39)
                                                                                   tGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc                 <  1:1306993/100‑1 (MQ=255)
                                                                                   tGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc                 <  1:3943781/100‑1 (MQ=255)
                                                                                   cagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc                 >  2:4165289/5‑100 (MQ=39)
                                                                                   cagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc                 >  2:2127915/5‑100 (MQ=39)
                                                                                    ggCAAAATGACAATTTTGTCATTTTTCTGGCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAGATCAACAGAAATTAGGAGCCTATGgct                >  2:4115709/1‑100 (MQ=255)
                                                                                    agaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtct                >  2:3159077/4‑100 (MQ=255)
                                                                                     gagaaaTGACAAAATTGTCATTTTTCTGTCAACTAAAAATATGAGCCTGTCTAGTAAATTTGGCGGCATAAAATAACCAGAAATAATGACCCTAtttatc               >  2:2747440/4‑95 (MQ=255)
                                                                                     gaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtctc               >  1:2471749/3‑100 (MQ=255)
                                                                                     gaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtctc               >  1:2154262/3‑100 (MQ=255)
                                                                                     gaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtctc               >  1:4655/3‑100 (MQ=255)
                                                                                     gCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtctc               >  1:1466271/1‑100 (MQ=255)
                                                                                      cAAAATGACAATTTTGTCATTTTTCTCTCACCGTAAAAGCAGATCCTGGCGAGGAAATTTGTCGTCATAAAATCACCAGAAATTATGAGCCTATGTCTcc              >  2:345487/1‑100 (MQ=255)
                                                                                       aaaaTGATTATTTTGACGGTTTTCTGTAACCGTAAAATAAGAGCCTGGGCAGTAAAGTTGGCCGCATATAATCACCAGAACTTATGAGCCTATGTCTCCt             <  2:3886765/100‑1 (MQ=255)
                                                                                        aaaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCtt            >  1:2674152/1‑100 (MQ=255)
                                                                                         aaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTg           <  1:1508978/100‑1 (MQ=255)
                                                                                         aaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTg           <  1:3358759/100‑1 (MQ=255)
                                                                                            gACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTaa        <  1:3370103/100‑1 (MQ=255)
                                                                                            gACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTaa        <  2:3191412/100‑1 (MQ=255)
                                                                                              cAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAAc      >  2:3782881/1‑100 (MQ=255)
                                                                                              cAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAAc      >  2:1914080/1‑100 (MQ=255)
                                                                                                aTTTTGTCATTTTTCTGTCACCGAAAAATCAGACTCTTTCCAGTAAAGTTGGGTTCATAAAATCACCACAAATTATGAGCCTCTGTTTCTTTGTaa        <  2:1360180/96‑1 (MQ=255)
                                                                                                  tttGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCt  >  2:3288585/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCT  >  NC_000913/595422‑595619

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