Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 595,519 T→A 18.9% 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,5190TA18.9% 339.2 / 72.1 150intergenic (‑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 (55/65);  new base A (18/10);  total (74/76)
Fisher's exact test for biased strand distribution p-value = 9.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

Only 100 of 159 total aligned reads displayed.

ATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAAC  >  NC_000913/595420‑595615
                                                                                                    |                                                                                                
  cTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                >  2:2319983/1‑100 (MQ=255)
  cTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                >  1:626385/1‑100 (MQ=255)
  cTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                >  1:1240755/1‑100 (MQ=255)
  cTTCGACAATAAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                >  2:269281/1‑100 (MQ=255)
     cggcAATCAACCGGTTCAGATTTCCTTCGCAGGTTGGCCGGGGCATTTTTGGGGAGGGGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                <  1:2462466/95‑1 (MQ=255)
     cGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGt                                                                                             <  2:2506341/100‑1 (MQ=255)
      gACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGACCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTg                                                                                              >  1:1696464/1‑98 (MQ=255)
       aCAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCa                                                                                           <  1:530753/100‑1 (MQ=255)
       aCAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCa                                                                                           <  2:3391321/100‑1 (MQ=255)
         aaTCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAtt                                                                                         >  2:2127949/1‑100 (MQ=255)
          aTCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGTCAAAATGACAATTTTGTCAtt                                                                                         >  1:942304/1‑99 (MQ=255)
          aTCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAttt                                                                                        >  1:1564236/1‑100 (MQ=255)
          aTCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAttt                                                                                        >  1:2758308/1‑100 (MQ=255)
          aTCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttt                                                                                        >  2:3741710/1‑100 (MQ=255)
           tCAACAGTTTCATATTTCCTGCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAtttt                                                                                       <  1:3468105/100‑1 (MQ=255)
           tCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAtttt                                                                                       <  1:2757229/100‑1 (MQ=255)
             aaCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTc                                                                                     >  2:729237/1‑100 (MQ=255)
             aaCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTc                                                                                     >  2:3786380/1‑100 (MQ=255)
              actgTTTTATATTTCCTCCGCAGGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgcc                                                                                    <  1:44054/97‑4 (MQ=255)
              aCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCt                                                                                    <  1:664925/100‑1 (MQ=255)
              aCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCt                                                                                    <  2:2866081/100‑1 (MQ=255)
              aCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCt                                                                                    >  2:3212446/1‑100 (MQ=255)
              aCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCt                                                                                    >  2:1359486/1‑100 (MQ=255)
               cAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTg                                                                                   <  1:58966/100‑1 (MQ=255)
               cAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgcca                                                                                   <  2:2508590/100‑5 (MQ=39)
               cAGTTTCAAATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGTATCAGATGGCAATCGCTTATTGGCAAAATGTCAAAATTGTCAttttgcca                                                                                   <  1:3544261/100‑5 (MQ=255)
                aGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGt                                                                                  <  2:487953/100‑1 (MQ=255)
                   ttCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCAc                                                                               >  2:1196973/1‑100 (MQ=255)
                   ttCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCAc                                                                               <  2:1899042/100‑1 (MQ=255)
                   ttCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCAc                                                                               <  1:1454436/100‑1 (MQ=255)
                    tCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCAcc                                                                              <  2:3745808/100‑1 (MQ=255)
                    tCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgccaataag                                                                              <  1:2073831/100‑10 (MQ=38)
                      atatTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCgg                                                                            >  2:60441/1‑100 (MQ=255)
                      atatTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCgg                                                                            >  1:888651/1‑100 (MQ=255)
                      atatTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTAATTTTTCTGTCACCg                                                                             >  1:3267961/1‑99 (MQ=255)
                      atatTTCCAGCGCATATTGCCCGGGCAATTGCATATTACGGGGATGAGATGGCAATAGCGTAGTGGCAAAATGACAATTTTg                                                                                              <  1:3693412/82‑1 (MQ=255)
                         attCATCTGCACGTTGCCCGGGCAATTCAAAAATAGCGGGATCAGATGGCAGTCAATTATTAGAAAAATGACAATGTTGTCGTTTTACTgaaac                                                                               >  1:2309718/2‑90 (MQ=255)
                          ttCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaaaa                                                                        <  1:3085135/100‑1 (MQ=255)
                             ctccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCa                                                                     <  2:2768733/100‑1 (MQ=255)
                              tccGCATGTTGTCCGGGCAATTCTAGAGTAGTTTTATCAGATGGCAATCGCTTATTGGAAAAATGACAATTTTGTCATTTTTCTTTCACCGGAAAATCa                                                                     <  1:1229646/99‑1 (MQ=255)
                               ccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCaga                                                                   >  2:789645/1‑100 (MQ=255)
                                cGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCagag                                                                  <  2:75110/100‑1 (MQ=255)
                                cGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCagag                                                                  >  2:2070230/1‑100 (MQ=255)
                                 gCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGc                                                                 >  2:3258531/1‑100 (MQ=255)
                                      ttGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGc                                                            <  1:2406796/100‑1 (MQ=255)
                                           ccGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGt                                                        >  1:2071600/1‑100 (MQ=255)
                                           ccGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGt                                                        <  2:3480724/100‑1 (MQ=255)
                                           ccGGGCAATTCTAAAGTAGCGGGATCAGATGGCCATCACTTATTGGAAAAATGAAAATATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAg                                                         >  1:3860316/1‑99 (MQ=255)
                                            cGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTa                                                       >  2:3439059/1‑100 (MQ=255)
                                             gggCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGCAAAATCAGAGCCTGGCGAGTaa                                                      <  1:1908708/100‑1 (MQ=255)
                                               gCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAg                                                    <  1:3786348/100‑1 (MQ=255)
                                                   ttCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGt                                                   <  1:1418898/97‑1 (MQ=255)
                                                    tCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggc                                               <  1:1007915/100‑1 (MQ=255)
                                                     cTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcg                                              >  2:3879471/1‑100 (MQ=255)
                                                       agagTAGCGGTATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcggc                                            >  2:1959146/1‑100 (MQ=255)
                                                       agagTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcggc                                            >  1:3798858/1‑100 (MQ=255)
                                                           tAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATaa                                        <  2:3300784/100‑1 (MQ=255)
                                                             gCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATaaaa                                      <  1:3342737/100‑1 (MQ=255)
                                                             gCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATaaaa                                      <  2:1466480/100‑1 (MQ=255)
                                                               gggATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGAAAAGTTGGCGGCATAAAATc                                    >  2:3442338/1‑100 (MQ=255)
                                                                ggATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCa                                   <  2:3669421/100‑1 (MQ=255)
                                                                   tCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGGCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCa                                <  1:1583141/100‑1 (MQ=255)
                                                                    cAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAg                               >  1:1163469/1‑100 (MQ=255)
                                                                    cAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAg                               >  1:2738790/1‑100 (MQ=255)
                                                                     agctGGCACTCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGCCGAGTAAAGTTTACGGCCTAAAATCACCa                                >  1:752612/4‑98 (MQ=255)
                                                                      gATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGCGTAAAGTTGGCGGCATAAAATCACCAGaa                             >  1:2607587/1‑100 (MQ=255)
                                                                      gATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGaa                             >  1:1042382/1‑100 (MQ=255)
                                                                      gATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGaa                             >  1:2876750/1‑100 (MQ=255)
                                                                      gATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGaa                             >  1:2011732/1‑100 (MQ=255)
                                                                        tGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAAt                           <  2:3056676/100‑1 (MQ=255)
                                                                          gCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTa                         <  2:2841220/100‑1 (MQ=255)
                                                                           cAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTAt                        <  2:1908607/100‑1 (MQ=255)
                                                                              tCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATg                       <  1:1158156/98‑1 (MQ=255)
                                                                              tCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAg                     <  1:3260613/100‑1 (MQ=255)
                                                                                gCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGcc                   <  2:1342873/100‑1 (MQ=255)
                                                                                 cTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCt                  >  1:881540/1‑100 (MQ=255)
                                                                                    nacagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATg               >  1:3702549/7‑100 (MQ=38)
                                                                                    aTTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATg               <  2:948134/100‑1 (MQ=255)
                                                                                    aTTGGCAAAATGACAATTTTGTCATTTTTCTGACACCGGTAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATg               <  1:2727082/100‑1 (MQ=255)
                                                                                     ttGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGt              >  1:678727/1‑100 (MQ=255)
                                                                                      cagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc             >  2:263094/5‑100 (MQ=39)
                                                                                      cagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc             >  2:3659892/5‑100 (MQ=39)
                                                                                      cagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTCGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc             >  1:859085/5‑100 (MQ=39)
                                                                                       ggCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc             >  1:2292796/1‑99 (MQ=255)
                                                                                       agaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtct            >  2:2635736/4‑100 (MQ=255)
                                                                                        gaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAACTCAGAGCCTGGCGAGTAAAGTTAGCGGCATAAAAACACCAGAAATAATAAGCCTATGtct            >  1:446620/3‑99 (MQ=255)
                                                                                        gaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtctc           >  2:1706348/3‑100 (MQ=255)
                                                                                        gaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtctc           >  2:2673006/3‑100 (MQ=255)
                                                                                        gaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGAATAAAATCCCCAGAAATTATGAGCCTATGtct            >  1:1308789/3‑99 (MQ=255)
                                                                                         naaaaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTcc          >  1:745246/2‑100 (MQ=255)
                                                                                         cAAAATGACAATTTTGTCATTTTTCTTTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtctc           <  1:2952597/99‑1 (MQ=255)
                                                                                         cAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTcc          >  2:5341/1‑100 (MQ=255)
                                                                                         cAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTcc          <  1:1784161/100‑1 (MQ=255)
                                                                                          aaaaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCt         <  2:888659/100‑1 (MQ=255)
                                                                                          aaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCt         >  2:2247499/1‑100 (MQ=255)
                                                                                          aaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCt         >  2:2364344/1‑100 (MQ=255)
                                                                                            aaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTg       <  1:1658827/100‑1 (MQ=255)
                                                                                             aTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTg       <  2:466128/99‑1 (MQ=255)
                                                                                               gACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTaa    >  1:1821868/1‑100 (MQ=255)
                                                                                                 cAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAAc  >  1:94189/1‑100 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: