Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,185,875 (G)7→8 intergenic (‑79/‑42) rcnR ← / → rcnA transcriptional repressor of rcnA/membrane protein conferring nickel and cobalt resistance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,185,8681.G84.9% 210.8 / 32.3 73intergenic (‑72/‑49)rcnR/rcnAtranscriptional repressor of rcnA/membrane protein conferring nickel and cobalt resistance
Reads supporting (aligned to +/- strand):  ref base . (1/10);  new base G (57/5);  total (58/15)
Fisher's exact test for biased strand distribution p-value = 4.93e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGTTTATCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGT  >  NC_000913/2185769‑2185959
                                                                                                     |                                                                                                
tCTGTTTATCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACt                                                                                                   >  2:1715807/1‑100 (MQ=255)
   gTTTATCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTg‑‑gg                                                                                            >  1:2516146/1‑100 (MQ=255)
     ttATCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTg‑‑gggg                                                                                         >  2:2662869/1‑100 (MQ=255)
     ttATCACAAATAGTATGAGACATGGCAACTCCTGGTTAAAAAAAAAAAAAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTg‑‑gg                                                                                            <  2:367570/98‑1 (MQ=255)
      tATCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCATTATTAATCTACTGGGGGGTAGTATCAGGTACTg‑‑ggggg                                                                                        <  2:2731175/100‑1 (MQ=255)
       aTCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTg‑‑gggggg                                                                                       >  1:2283239/1‑100 (MQ=255)
       aTCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTg‑‑gggggg                                                                                       <  2:2166615/100‑1 (MQ=255)
       aTCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATATACTGGGGGCTAGTACCAGCCACTgcaggg                                                                                            >  1:530743/1‑94 (MQ=255)
        tCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTg‑‑ggggggg                                                                                      >  1:2315086/1‑99 (MQ=255)
         cACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTg‑‑ggggggga                                                                                     >  2:2759057/1‑98 (MQ=255)
         cACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTg‑‑ggggggga                                                                                     >  2:2131589/1‑98 (MQ=255)
         cACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTg‑‑ggggggga                                                                                     >  1:115372/1‑98 (MQ=255)
          acgaaatGTATGAAACACGGCAACACCTGGTTAAAAAGAAAAAGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTg‑‑ggggggga                                                                                     <  2:1769302/93‑3 (MQ=255)
           cGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTg‑‑gggggggagt                                                                                   >  2:1320099/1‑96 (MQ=255)
            ggATTGTAGGAGACAGGTCAACAGCTGGTTAACAAGAAGATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTg‑‑ggggggga                                                                                     <  2:1807920/97‑3 (MQ=255)
             gATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAg                                                                                 >  1:1473291/1‑100 (MQ=255)
             gATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAg                                                                                 >  1:28246/1‑100 (MQ=255)
              aTTGTATGAGACATGGAAAAACCGGATTAACAAGAATATGAAAAATCATAGCACTATTAATCTGCTGGGGGGTAGTATCAGGTACTg‑‑gggggggagta                                                                                  <  2:2245346/98‑6 (MQ=255)
                 gTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATc                                                                             >  2:1869773/1‑100 (MQ=255)
                 gTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATc                                                                             >  1:1615675/1‑100 (MQ=255)
                  tATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCa                                                                            >  2:518533/1‑100 (MQ=255)
                  tATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCa                                                                            >  1:1778368/1‑100 (MQ=255)
                   gtGAGAAATGGCACCACAGAGTTAACAAGAATATGAATACTCATAGCTCTATTAATCAACTGGGTGGTAGTATCAGGTACTG‑‑GAGGGGGAGTAGAATCAg                                                                           <  2:1298890/99‑1 (MQ=255)
                   aTGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAg                                                                           >  1:1158949/1‑100 (MQ=255)
                    tGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGa                                                                          >  2:595076/1‑100 (MQ=255)
                     gagaCATGTGAACACATGGTTAACAAAAATATGAAAAATTATAGCAGTATGAAGAGACAGGGGCGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCa                                                                            <  2:347648/97‑1 (MQ=255)
                      agaCATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGAtt                                                                        >  1:1979973/1‑100 (MQ=255)
                      agaCATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGAtt                                                                        >  1:1978878/1‑100 (MQ=255)
                      agaCATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGAt                                                                         >  2:1232229/1‑99 (MQ=255)
                      agaCATGGCAACACCTGGTTAAAAAGAATATGAAAACTCATAGCACTATTAATCTACTGGGGGGTAGTAGCAGGGACTGG‑‑GGGGGGAGTAGAATcgggtt                                                                        >  1:1983225/1‑95 (MQ=255)
                          aTGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCg                                                                    >  1:1358585/1‑100 (MQ=255)
                           tgtcAACACCTGGGTAAAAAGAATATGAAAAATAATAACTGTATTAAGCTACTGGGGGGTAGTATCAAGGACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGa                                                                   <  2:2298658/97‑1 (MQ=255)
                           tGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGa                                                                   >  1:1465833/1‑100 (MQ=255)
                            ggCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGaa                                                                  >  2:1990265/1‑100 (MQ=255)
                            ggCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGaa                                                                  >  2:669994/1‑100 (MQ=255)
                              cgaccccTAGTTAACAAGGAGATGAAAGATGATAGCTCCATTAGTGTACTGGGGGGTAGTATCAGTTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAAtt                                                                <  2:1333814/95‑1 (MQ=255)
                                 cacCTGGTTAACACGAATATGAAAAATCATAACACTATTAATGTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAAt                                                             >  1:1950535/1‑100 (MQ=255)
                                  acCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATa                                                            >  1:1418421/1‑100 (MQ=255)
                                   ccTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATAc                                                           >  2:526809/1‑100 (MQ=255)
                                     tGGGTAACAAGAATATGAAAAAGGATAGAGGTATTGAGAGAGTGGGGGGTAGTATCAAGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAAt                                                             <  2:2332329/96‑1 (MQ=255)
                                      ggTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTaagatcgg                                                        >  2:1108099/1‑94 (MQ=255)
                                          aactaaaaTATGAAAAATGATAGCATTGATGATGTATTGGTGGGTATTAACAGGTACTG‑‑GGGGGGAAGTAGAATCAGATTGCCGAATTAATACTAAGAat                                                    <  2:22517/94‑1 (MQ=255)
                                           aCAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAatt                                                   >  1:2617986/1‑100 (MQ=255)
                                           aCAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAatt                                                   >  2:28375/1‑100 (MQ=255)
                                           aCAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAatt                                                   >  1:18519/1‑100 (MQ=255)
                                           aCAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAatt                                                   >  1:102851/1‑100 (MQ=255)
                                            cAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGCGTAGAATCAGATTGCCGAATTAATACTAAGAatta                                                  >  2:2321073/1‑100 (MQ=255)
                                            cAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAatta                                                  >  2:904093/1‑100 (MQ=255)
                                             aaGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAattat                                                 >  2:953980/1‑100 (MQ=255)
                                              aGAATATGAAAAATTATGGTATTATGAAGATACGGGGTGGTAGTATCAGGTAATG‑‑GGAGGGAAGTAGAATCAGATTGCCGAATTAATACTAAGAattat                                                 <  2:2315071/99‑1 (MQ=255)
                                               gAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAAttatta                                               >  2:1882794/1‑100 (MQ=255)
                                                aaTATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATtattat                                              >  2:575542/1‑100 (MQ=255)
                                                  gatGAAATATCATAGGTCTATGGATTTACTGGAGAGTATTATCAAGTACTGATAATGAAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCa                                            <  2:1164616/99‑1 (MQ=255)
                                                   atGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCAt                                           >  2:993379/1‑100 (MQ=255)
                                                    taaaaaaTCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATg                                          >  2:1385546/3‑100 (MQ=255)
                                                    tGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATg                                          >  2:21589/1‑100 (MQ=255)
                                                    tGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATg                                          >  1:2343633/1‑100 (MQ=255)
                                                    tGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATg                                          >  2:2658733/1‑100 (MQ=255)
                                                     gAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGa                                         >  2:1702298/1‑100 (MQ=255)
                                                       aaaaTCATAGGAGTATAAAGGAACGGGAGGGTAGTATAAAGTAATG‑‑GGGGGGAAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGa                                         <  2:1323015/98‑1 (MQ=255)
                                                        aaaTTATAGCACTATTAATTTACTGGGGGGTAGTATCAAATAATG‑‑GAAGAGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCa                                            <  2:1562097/94‑1 (MQ=255)
                                                         aaTCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGa                                     >  2:1503271/1‑100 (MQ=255)
                                                            cATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAAt                                   >  1:1928614/1‑99 (MQ=255)
                                                                gCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACa                              >  2:750592/1‑100 (MQ=255)
                                                                gCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACa                              >  2:1317860/1‑100 (MQ=255)
                                                                gCACTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACa                              >  1:2107469/1‑100 (MQ=255)
                                                                   cTATTAATCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAAct                           >  1:1665581/1‑100 (MQ=255)
                                                                       taatgtgcgggggggTAGTGTCAGGGAGTGGTGGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAActct                         <  1:2067278/90‑1 (MQ=255)
                                                                        aaTCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACtcttct                      >  2:1913672/1‑100 (MQ=255)
                                                                        aaTCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACtcttct                      >  1:1042151/1‑100 (MQ=255)
                                                                        aaTCTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACtcttct                      >  2:1914314/1‑100 (MQ=255)
                                                                           cTACTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCa                   >  2:1059362/1‑100 (MQ=255)
                                                                              cTGGGGGGTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCa                >  2:1374980/1‑100 (MQ=255)
                                                                                ggggggTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAg              >  2:861084/1‑100 (MQ=255)
                                                                                ggggggTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAg              >  2:1112798/1‑100 (MQ=255)
                                                                                atctggTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAg              <  2:820161/96‑1 (MQ=255)
                                                                                  ggggTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGa            >  2:533348/1‑100 (MQ=255)
                                                                                    ggTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAgagat          >  1:820162/1‑96 (MQ=255)
                                                                                    ggTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGaaa          >  1:136980/1‑100 (MQ=255)
                                                                                    ggTAGTATCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGaaa          >  1:136239/1‑100 (MQ=255)
                                                                                     gtagtaTCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAAc         >  1:2868534/1‑100 (MQ=255)
                                                                                       attaTCAGGTAGTG‑‑GGTAGGAAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACg        <  2:494064/97‑1 (MQ=255)
                                                                                         taTCAGGTAAGTG‑‑AGGTGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAAc         <  2:2516130/96‑1 (MQ=255)
                                                                                          aTCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTg    >  1:440699/1‑100 (MQ=255)
                                                                                          aTCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTg    >  1:441001/1‑100 (MQ=255)
                                                                                          aTCAGGTACTG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTgg   >  2:2580395/1‑100 (MQ=255)
                                                                                           tCAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTgg   >  1:2308411/1‑100 (MQ=255)
                                                                                            cAGGTACTGG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGt  >  2:1154040/1‑100 (MQ=255)
                                                                                                tACTG‑‑GTGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTgg   <  1:1112802/95‑1 (MQ=255)
                                                                                                     |                                                                                                
TCTGTTTATCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGTACTG‑‑GGGGGGAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGT  >  NC_000913/2185769‑2185959

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: