Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 595,519 T→A 17.1% intergenic (‑76/‑81) 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,5190TA17.1% 323.5 / 56.7 140intergenic (‑76/‑81)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 (62/54);  new base A (14/10);  total (76/64)
Fisher's exact test for biased strand distribution p-value = 8.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-01

Only 100 of 144 total aligned reads displayed.

ATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAAC  >  NC_000913/595420‑595615
                                                                                                   |                                                                                                
aTCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAt                                                                                                  >  1:3843199/1‑100 (MQ=255)
aTCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAt                                                                                                  >  1:3843419/1‑100 (MQ=255)
  cTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                <  1:166603/100‑1 (MQ=255)
    tCGACTATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAGTCGCTTATTGGCAAAATGACAATTTTg                                                                                              <  1:4064515/100‑1 (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)
           tCAACAGTTTCATAGTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttt                                                                                        <  2:1952881/99‑1 (MQ=255)
            cAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAttttt                                                                                      >  2:412346/1‑100 (MQ=255)
             aaCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTc                                                                                     >  1:1731204/1‑100 (MQ=255)
             aaCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgc                                                                                     <  2:2170128/100‑3 (MQ=255)
             aaCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgc                                                                                     <  2:1198712/100‑3 (MQ=255)
             aaCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgc                                                                                     <  1:3479928/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)
                    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)
                     cATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCg                                                                             >  2:3231698/1‑100 (MQ=255)
                      atatTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTACTGTCACCgg                                                                            >  2:2749863/1‑100 (MQ=255)
                          ttCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaaaa                                                                        <  2:3724272/100‑1 (MQ=255)
                          ttCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaaaa                                                                        >  1:1693512/1‑100 (MQ=255)
                           ncctccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAAt                                                                       >  1:3047161/2‑100 (MQ=255)
                            cctgcGCATGTTGCCCGGGATATTCTAGAGTAGCGGATCAGCTGGCAATCGCTTAGTGGCAAAACGACAATTTTGTCATTTTTCTGTCACCGGAAAATCa                                                                     <  1:2351893/100‑1 (MQ=255)
                            cctccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAAt                                                                       >  2:3439277/1‑99 (MQ=255)
                            cctccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATc                                                                      <  1:164944/100‑1 (MQ=255)
                                cGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCaga                                                                   >  2:2357077/1‑99 (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)
                                            gggCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTaa                                                      <  1:325470/100‑1 (MQ=255)
                                              gCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAg                                                    <  1:224348/100‑1 (MQ=255)
                                               cAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGt                                                   >  2:2863872/1‑100 (MQ=255)
                                               cAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGt                                                   <  1:2233893/100‑1 (MQ=255)
                                                   tCTAGAGTAGCGGGATCAGATGGCTATCGCTTATTGGCAAAATGACAATTTTGTCAGTTGTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggc                                               <  2:1014859/100‑1 (MQ=255)
                                                    cTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcg                                              >  1:553515/1‑100 (MQ=255)
                                                     tAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcgg                                             >  1:2293994/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)
                                                              gggATCAGATGGCCATCGCTTATTGGCAAAATGACAAATTTGTTATTTTTCTGTCACCGGAAAATCATAGCCTTGCGAGTAAAGTTATCGTCATAAAATc                                    >  2:2234960/1‑100 (MQ=255)
                                                              gggATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAAt                                     >  1:1903319/1‑99 (MQ=255)
                                                                gATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCACAAAATCAc                                  <  2:2855064/100‑1 (MQ=255)
                                                                 aTCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCAcc                                 <  1:3746015/100‑1 (MQ=255)
                                                                    aGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGa                              <  2:1363135/100‑1 (MQ=255)
                                                                    aGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGa                              >  2:3137516/1‑100 (MQ=255)
                                                                    aGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGa                              <  1:1968158/100‑1 (MQ=255)
                                                                     gATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGaa                             >  2:2583943/1‑100 (MQ=255)
                                                                     gATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGaa                             >  2:3908380/1‑100 (MQ=255)
                                                                     gATGGCAATCGCTTATTAGCAAAATGACAATTTAGTCATTTTTCTGTCACGGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCaatag                               >  2:275946/1‑94 (MQ=255)
                                                                      aTGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGaaa                            >  1:2580246/1‑100 (MQ=255)
                                                                        ggCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAAtt                          <  2:2835682/100‑1 (MQ=255)
                                                                          cAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTAt                        <  2:755083/100‑1 (MQ=255)
                                                                              cGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGc                    >  1:3350248/1‑100 (MQ=255)
                                                                              cGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGc                    >  2:3190043/1‑100 (MQ=255)
                                                                              cGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGc                    >  2:1487609/1‑100 (MQ=255)
                                                                               gCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGcc                   <  2:2362790/100‑1 (MQ=255)
                                                                               gCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGcc                   >  1:350831/1‑100 (MQ=255)
                                                                               gCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGcc                   >  2:1235029/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:2127915/5‑100 (MQ=39)
                                                                                     cagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc             >  2:4165289/5‑100 (MQ=39)
                                                                                      ggCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtct            >  1:800085/1‑100 (MQ=255)
                                                                                      ggCAAAATGACAATTTTGTCATTTTTCTGGCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAGATCAACAGAAATTAGGAGCCTATGgct            >  2:4115709/1‑100 (MQ=255)
                                                                                       gaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtctc           >  1:4655/3‑100 (MQ=255)
                                                                                       gaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtctc           >  1:2471749/3‑100 (MQ=255)
                                                                                       gaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtctc           >  1:2154262/3‑100 (MQ=255)
                                                                                       gCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtctc           >  1:1466271/1‑100 (MQ=255)
                                                                                         aaaaTGATTATTTTGACGGTTTTCTGTAACCGTAAAATAAGAGCCTGGGCAGTAAAGTTGGCCGCATATAATCACCAGAACTTATGAGCCTATGTCTCCt         <  2:3886765/100‑1 (MQ=255)
                                                                                         aaaaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCt         <  1:3120494/100‑1 (MQ=255)
                                                                                          aaaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCtt        >  1:2674152/1‑100 (MQ=255)
                                                                                           aaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTg       <  1:3358759/100‑1 (MQ=255)
                                                                                           aaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTg       <  1:1508978/100‑1 (MQ=255)
                                                                                              gACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTaa    <  2:3191412/100‑1 (MQ=255)
                                                                                              gACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTaa    <  1:3370103/100‑1 (MQ=255)
                                                                                                cAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAAc  >  2:3782881/1‑100 (MQ=255)
                                                                                                cAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAAc  >  1:3612681/1‑100 (MQ=255)
                                                                                                cAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAAc  >  2:1914080/1‑100 (MQ=255)
                                                                                                  aTTTTGTCATTTTTCTGTCACCGAAAAATCAGACTCTTTCCAGTAAAGTTGGGTTCATAAAATCACCACAAATTATGAGCCTCTGTTTCTTTGTaa    <  2:1360180/96‑1 (MQ=255)
                                                                                                   |                                                                                                
ATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAAC  >  NC_000913/595420‑595615

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: