Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 595,520 T→A 16.7% 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.7% 209.9 / 33.9 89intergenic (‑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 (41/33);  new base A (11/4);  total (52/37)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

CTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCT  >  NC_000913/595422‑595619
                                                                                                  |                                                                                                   
cTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                    >  2:1225544/1‑100 (MQ=255)
cTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                    >  1:2078029/1‑100 (MQ=255)
cTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                    <  1:52226/100‑1 (MQ=255)
   cGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGt                                                                                                 >  1:935689/1‑100 (MQ=255)
    gACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTc                                                                                                >  2:2546410/1‑100 (MQ=255)
      cAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAt                                                                                              >  1:2261459/1‑100 (MQ=255)
       aaTCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAtt                                                                                             <  2:104480/100‑1 (MQ=255)
        aTCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAttt                                                                                            >  2:2483394/1‑100 (MQ=255)
          cAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAttttt                                                                                          >  2:139418/1‑100 (MQ=255)
          cAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAttttt                                                                                          >  2:352505/1‑100 (MQ=255)
           aaCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTc                                                                                         >  2:1584557/1‑100 (MQ=255)
           aaCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgc                                                                                         <  2:1143693/100‑3 (MQ=255)
            ccAGTTTCATATTTCCTCCGGATGTTGCCCGAACAGATCTAGAGTAGCAGGATCAGATTGCAATCGCGTAGTGGCCAAATAACTACTTTTTCATtcttct                                                                                        >  1:821644/2‑100 (MQ=255)
            aCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCt                                                                                        <  2:2099983/100‑1 (MQ=255)
            aCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCt                                                                                        <  1:1122894/100‑1 (MQ=255)
            aCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgcc                                                                                        <  1:1128027/100‑4 (MQ=255)
            aCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgc                                                                                         <  1:1251558/99‑3 (MQ=255)
            aCAGTCTCATAATCCCTCCGCATGTTGCCAGGGCAATTCTAGAGTAGCGGCATCAGATGGCAAGCGATCATTCGAAAAATGACAACTTTATCCttt                                                                                            >  1:821482/1‑96 (MQ=255)
               tttccATATTCCCTCCGCATGCCGCCCGGGCAATTCTACAGTATCGGTATCAGATTGCAATCGCTTATTGGCTATATGACAATTTTGTCATTTTTCTGTc                                                                                     <  1:2550634/96‑1 (MQ=255)
               gTTTCATATTTCCTCCGCATGTTGCCCGGGTAATTCTAGGGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTg                                                                                       <  2:2353408/98‑1 (MQ=255)
                 ttCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTATCGGGATAAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCAc                                                                                   >  1:1978018/1‑100 (MQ=255)
                 ttCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCAc                                                                                   >  2:842951/1‑100 (MQ=255)
                 ttCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCAc                                                                                   >  2:1725001/1‑100 (MQ=255)
                 ttCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGCCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCAc                                                                                   >  1:2279749/1‑100 (MQ=255)
                   cATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGATTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCAcc                                                                                  <  1:2478655/99‑1 (MQ=255)
                   cATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCg                                                                                 >  2:1121435/1‑100 (MQ=255)
                   cATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCg                                                                                 <  1:928540/100‑1 (MQ=255)
                   cATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCg                                                                                 >  1:1824291/1‑100 (MQ=255)
                      atTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaa                                                                              >  1:745844/1‑100 (MQ=255)
                      atTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaa                                                                              <  1:537719/100‑1 (MQ=255)
                        ttCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaaaa                                                                            >  1:1914733/1‑100 (MQ=255)
                         tcctccGCATGTTGCCCGGGCACTTCTAGACTTTCTCGATAAGATAGCAATCGCTTACTGGCAAAATCACAGTTTTGACATTTATCTGTCAGCGGaaaa                                                                            >  1:389755/1‑99 (MQ=255)
                         tcctccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAAt                                                                           <  1:1282693/100‑1 (MQ=255)
                         tcctccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAAt                                                                           >  1:1374068/1‑100 (MQ=255)
                                  tGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTg                                                                  <  2:1416739/100‑1 (MQ=255)
                                     tGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCg                                                               >  2:2430607/1‑100 (MQ=255)
                                       cccGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAg                                                             <  2:573953/100‑1 (MQ=255)
                                            gCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGACAATCAGAGCCTGGCGAGTAAAg                                                        <  1:1501403/100‑1 (MQ=255)
                                            gCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTaaa                                                         >  1:2096707/1‑99 (MQ=255)
                                            gCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAg                                                        <  2:1675614/100‑1 (MQ=255)
                                               aTTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTg                                                     >  2:954734/1‑100 (MQ=255)
                                                ttCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTgg                                                    >  1:1683057/1‑100 (MQ=255)
                                                 tCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggc                                                   >  2:1473002/1‑100 (MQ=255)
                                                  cTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcg                                                  <  2:2534631/100‑1 (MQ=255)
                                                  cTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcg                                                  >  2:494592/1‑100 (MQ=255)
                                                   tAGAGTAGCGGTATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGGCATTTTTCTGGCACCGGAAAATCAGAGCCTGCCGAGTAAAGTTggcgg                                                 <  1:591399/100‑1 (MQ=255)
                                                   tAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcgg                                                 <  1:1879153/100‑1 (MQ=255)
                                                   tAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcgg                                                 <  1:2279177/100‑1 (MQ=255)
                                                     gagTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCa                                               >  1:1678638/1‑100 (MQ=255)
                                                         aGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATaaa                                           >  2:1037633/1‑100 (MQ=255)
                                                         aGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATaaa                                           >  2:1182366/1‑100 (MQ=255)
                                                          gCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATaaaa                                          >  1:2403754/1‑100 (MQ=255)
                                                           cGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAAt                                         <  2:2070390/100‑1 (MQ=255)
                                                           cGGGATCAGATGGCAATCGCTTATGGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAAt                                         <  2:1314843/100‑1 (MQ=255)
                                                            gggATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATc                                        >  2:1979556/1‑100 (MQ=255)
                                                             ggATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCa                                       <  1:352503/100‑1 (MQ=255)
                                                              gATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCAc                                      <  1:1182362/100‑1 (MQ=255)
                                                               aTCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCAcc                                     >  1:1046639/1‑100 (MQ=255)
                                                                tCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCa                                    <  2:491546/100‑1 (MQ=255)
                                                                 cAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAg                                   <  2:829052/100‑1 (MQ=255)
                                                                  aGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGa                                  >  2:85705/1‑100 (MQ=255)
                                                                   gATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGaa                                 <  1:842948/100‑1 (MQ=255)
                                                                      ggCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAAtt                              <  1:270256/100‑1 (MQ=255)
                                                                           tCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAg                         <  2:487403/100‑1 (MQ=255)
                                                                             gCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGcc                       >  1:2403893/1‑100 (MQ=255)
                                                                               ttATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTa                     <  1:1331103/100‑1 (MQ=255)
                                                                                 aTTGGCAAAGTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATg                   <  1:2594090/100‑1 (MQ=255)
                                                                                 aTTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATg                   >  2:1405665/1‑100 (MQ=255)
                                                                                 aTTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATg                   <  1:2503048/100‑1 (MQ=255)
                                                                                  ttGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGt                  >  1:1016003/1‑100 (MQ=255)
                                                                                  acagaaaaaTGACAAAATTGCCATTTTTATGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATAATGAGCCTATGt                  >  1:1974680/6‑100 (MQ=255)
                                                                                   tGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc                 >  2:2358175/1‑100 (MQ=255)
                                                                                   nagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc                 >  2:719682/5‑100 (MQ=255)
                                                                                   cagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc                 >  2:2337622/5‑100 (MQ=39)
                                                                                   cagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc                 >  1:2033484/5‑100 (MQ=39)
                                                                                   cagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc                 >  2:43383/5‑100 (MQ=39)
                                                                                   cagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGt                  >  1:206404/5‑99 (MQ=39)
                                                                                    ggCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtct                >  2:634689/1‑100 (MQ=255)
                                                                                    ggCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATAACCAGAACTTATGGGCCTATGtct                >  1:2556729/1‑100 (MQ=255)
                                                                                     gaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtctc               >  2:2221968/3‑100 (MQ=255)
                                                                                     gaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtctc               >  2:2448389/3‑100 (MQ=255)
                                                                                     gaaaaaTGACAAAATTGTCATTTTTCTGTCAACGGAAAATCAGAGCCTGGCGAGTAACGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtctc               >  1:2122228/3‑100 (MQ=255)
                                                                                      aaaaaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTcc              >  1:378816/2‑100 (MQ=255)
                                                                                       aaaaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCt             <  2:2021961/100‑1 (MQ=255)
                                                                                       aaaaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCt             <  1:361132/100‑1 (MQ=255)
                                                                                       aaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCt             >  2:464010/1‑100 (MQ=255)
                                                                                        aaaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCtt            <  2:1683066/100‑1 (MQ=255)
                                                                                          aTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGt          >  2:1816737/1‑100 (MQ=255)
                                                                                           tGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTa         <  2:769192/100‑1 (MQ=255)
                                                                                              cAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAAc      <  2:2207938/100‑1 (MQ=255)
                                                                                               aaTTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACt     >  1:2144760/1‑100 (MQ=255)
                                                                                                aTTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACtt    >  1:412754/1‑100 (MQ=255)
                                                                                                 ttttGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTc   <  1:2430592/100‑1 (MQ=255)
                                                                                                  tttGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTACAGTTGGCGGCATAAAATCACCAGAAATTATGAGCTTATGTCTCCTTATAAACTTCt  >  1:619415/1‑100 (MQ=255)
                                                                                                  tttGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCt  <  2:2382785/100‑1 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: