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

TCTGTTTATCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACT‑G‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGT  >  NC_000913/2185769‑2185959
                                                                                                     |                                                                                                
tCTGTTTATCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACt                                                                                                   >  2:1715807/1‑100 (MQ=255)
   gTTTATCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACT‑g‑‑g‑g                                                                                            >  1:2516146/1‑100 (MQ=255)
     ttATCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACT‑g‑‑g‑g‑gg                                                                                         >  2:2662869/1‑100 (MQ=255)
     ttATCACAAATAGTATGAGACATGGCAACTCCTGGTTAAAAAAAAAAAAAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACT‑g‑‑g‑g                                                                                            <  2:367570/98‑1 (MQ=255)
      tATCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCATTATTAATCTACTGGGGGGTAGTATCAGGT‑ACT‑g‑‑g‑g‑ggg                                                                                        <  2:2731175/100‑1 (MQ=255)
       aTCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACT‑g‑‑g‑g‑gggg                                                                                       >  1:2283239/1‑100 (MQ=255)
       aTCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACT‑g‑‑g‑g‑gggg                                                                                       <  2:2166615/100‑1 (MQ=255)
       aTCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATATACTGGGGGCTAGTACCAGCC‑ACT‑gcaggg                                                                                            >  1:530743/1‑94 (MQ=255)
        tCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACT‑g‑‑g‑g‑ggggg                                                                                      >  1:2315086/1‑99 (MQ=255)
         cACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACT‑g‑‑g‑g‑ggggga                                                                                     >  2:2759057/1‑98 (MQ=255)
         cACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACT‑g‑‑g‑g‑ggggga                                                                                     >  2:2131589/1‑98 (MQ=255)
         cACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACT‑g‑‑g‑g‑ggggga                                                                                     >  1:115372/1‑98 (MQ=255)
          acgaaatGTATGAAACACGGCAACACCTGGTTAAAAAGAAAAAGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACT‑g‑‑g‑g‑ggggga                                                                                     <  2:1769302/93‑3 (MQ=255)
           cGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACT‑g‑‑g‑g‑gggggagt                                                                                   >  2:1320099/1‑96 (MQ=255)
            ggATTGTAGGAGACAGGTCAACAGCTGGTTAACAAGAAGATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACT‑g‑‑g‑g‑ggggga                                                                                     <  2:1807920/97‑3 (MQ=255)
             gATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAg                                                                                 >  1:1473291/1‑100 (MQ=255)
             gATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAg                                                                                 >  1:28246/1‑100 (MQ=255)
              aTTGTATGAGACATGGAAAAACCGGATTAACAAGAATATGAAAAATCATAGCACTATTAATCTGCTGGGGGGTAGTATCAGGT‑ACT‑g‑‑g‑g‑gggggagta                                                                                  <  2:2245346/98‑6 (MQ=255)
                 gTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATc                                                                             >  2:1869773/1‑100 (MQ=255)
                 gTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATc                                                                             >  1:1615675/1‑100 (MQ=255)
                  tATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCa                                                                            >  2:518533/1‑100 (MQ=255)
                  tATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCa                                                                            >  1:1778368/1‑100 (MQ=255)
                   gtGAGAAATGGCACCACAGAGTTAACAAGAATATGAATACTCATAGCTCTATTAATCAACTGGGTGGTAGTATCAGGT‑ACT‑G‑‑GAG‑GGGG‑AGTAGAATCAg                                                                           <  2:1298890/99‑1 (MQ=255)
                   aTGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAg                                                                           >  1:1158949/1‑100 (MQ=255)
                    tGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGa                                                                          >  2:595076/1‑100 (MQ=255)
                     gagaCATGTGAACACATGGTTAACAAAAATATGAAAAATTATAGCAGTATGAAGAGACAGGGGCGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCa                                                                            <  2:347648/97‑1 (MQ=255)
                      agaCATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGAtt                                                                        >  1:1979973/1‑100 (MQ=255)
                      agaCATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGAtt                                                                        >  1:1978878/1‑100 (MQ=255)
                      agaCATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGAt                                                                         >  2:1232229/1‑99 (MQ=255)
                      agaCATGGCAACACCTGGTTAAAAAGAATATGAAAACTCATAGCACTATTAATCTACTGGGGGGTAGTAGCAGGG‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATcgggtt                                                                        >  1:1983225/1‑95 (MQ=255)
                          aTGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCg                                                                    >  1:1358585/1‑100 (MQ=255)
                           tgtcAACACCTGGGTAAAAAGAATATGAAAAATAATAACTGTATTAAGCTACTGGGGGGTAGTATCAAGG‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGa                                                                   <  2:2298658/97‑1 (MQ=255)
                           tGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGa                                                                   >  1:1465833/1‑100 (MQ=255)
                            ggCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGaa                                                                  >  2:1990265/1‑100 (MQ=255)
                            ggCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGaa                                                                  >  2:669994/1‑100 (MQ=255)
                              cgaccccTAGTTAACAAGGAGATGAAAGATGATAGCTCCATTAGTGTACTGGGGGGTAGTATCAGTT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAAtt                                                                <  2:1333814/95‑1 (MQ=255)
                                 cacCTGGTTAACACGAATATGAAAAATCATAACACTATTAATGTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAAt                                                             >  1:1950535/1‑100 (MQ=255)
                                  acCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATa                                                            >  1:1418421/1‑100 (MQ=255)
                                   ccTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATAc                                                           >  2:526809/1‑100 (MQ=255)
                                     tGGGTAACAAGAATATGAAAAAGGATAGAGGTATTGAGAGAGTGGGGGGTAGTATCAAGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAAt                                                             <  2:2332329/96‑1 (MQ=255)
                                      ggTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTaagatcgg                                                        >  2:1108099/1‑94 (MQ=255)
                                          aactaaaaTATGAAAAATGATAGCATTGATGATGTATTGGTGGGTATTAACAGGT‑ACT‑G‑‑G‑G‑GGGGAAGTAGAATCAGATTGCCGAATTAATACTAAGAat                                                    <  2:22517/94‑1 (MQ=255)
                                           aCAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAatt                                                   >  1:2617986/1‑100 (MQ=255)
                                           aCAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAatt                                                   >  2:28375/1‑100 (MQ=255)
                                           aCAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAatt                                                   >  1:18519/1‑100 (MQ=255)
                                           aCAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAatt                                                   >  1:102851/1‑100 (MQ=255)
                                            cAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑CGTAGAATCAGATTGCCGAATTAATACTAAGAatta                                                  >  2:2321073/1‑100 (MQ=255)
                                            cAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAatta                                                  >  2:904093/1‑100 (MQ=255)
                                             aaGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAattat                                                 >  2:953980/1‑100 (MQ=255)
                                              aGAATATGAAAAATTATGGTATTATGAAGATACGGGGTGGTAGTATCAGGT‑AAT‑G‑‑G‑GAGGGA‑AGTAGAATCAGATTGCCGAATTAATACTAAGAattat                                                 <  2:2315071/99‑1 (MQ=255)
                                               gAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAAttatta                                               >  2:1882794/1‑100 (MQ=255)
                                                aaTATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATtattat                                              >  2:575542/1‑100 (MQ=255)
                                                  gatGAAATATCATAGGTCTATGGATTTACTGGAGAGTATTATCAAGT‑ACT‑GA‑T‑A‑ATGA‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCa                                            <  2:1164616/99‑1 (MQ=255)
                                                   atGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCAt                                           >  2:993379/1‑100 (MQ=255)
                                                    taaaaaaTCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATg                                          >  2:1385546/3‑100 (MQ=255)
                                                    tGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATg                                          >  2:21589/1‑100 (MQ=255)
                                                    tGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATg                                          >  1:2343633/1‑100 (MQ=255)
                                                    tGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATg                                          >  2:2658733/1‑100 (MQ=255)
                                                     gAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGa                                         >  2:1702298/1‑100 (MQ=255)
                                                       aaaaTCATAGGAGTATAAAGGAACGGGAGGGTAGTATAAAGT‑AAT‑G‑‑G‑G‑GGGGAAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGa                                         <  2:1323015/98‑1 (MQ=255)
                                                        aaaTTATAGCACTATTAATTTACTGGGGGGTAGTATCAAAT‑AAT‑G‑‑GAA‑GAGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCa                                            <  2:1562097/94‑1 (MQ=255)
                                                         aaTCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGa                                     >  2:1503271/1‑100 (MQ=255)
                                                            cATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAAt                                   >  1:1928614/1‑99 (MQ=255)
                                                                gCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACa                              >  2:750592/1‑100 (MQ=255)
                                                                gCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACa                              >  2:1317860/1‑100 (MQ=255)
                                                                gCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACa                              >  1:2107469/1‑100 (MQ=255)
                                                                   cTATTAATCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAAct                           >  1:1665581/1‑100 (MQ=255)
                                                                       taatgtgcgggggggTAGTGTCAGGG‑AGT‑GGTG‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAActct                         <  1:2067278/90‑1 (MQ=255)
                                                                        aaTCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACtcttct                      >  2:1913672/1‑100 (MQ=255)
                                                                        aaTCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACtcttct                      >  1:1042151/1‑100 (MQ=255)
                                                                        aaTCTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACtcttct                      >  2:1914314/1‑100 (MQ=255)
                                                                           cTACTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCa                   >  2:1059362/1‑100 (MQ=255)
                                                                              cTGGGGGGTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCa                >  2:1374980/1‑100 (MQ=255)
                                                                                ggggggTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAg              >  2:861084/1‑100 (MQ=255)
                                                                                ggggggTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAg              >  2:1112798/1‑100 (MQ=255)
                                                                                atctggTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAg              <  2:820161/96‑1 (MQ=255)
                                                                                  ggggTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGa            >  2:533348/1‑100 (MQ=255)
                                                                                    ggTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAgagat          >  1:820162/1‑96 (MQ=255)
                                                                                    ggTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGaaa          >  1:136980/1‑100 (MQ=255)
                                                                                    ggTAGTATCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGaaa          >  1:136239/1‑100 (MQ=255)
                                                                                     gtagtaTCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAAc         >  1:2868534/1‑100 (MQ=255)
                                                                                       attaTCAGGT‑AGT‑G‑‑G‑GTAGGA‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACg        <  2:494064/97‑1 (MQ=255)
                                                                                         taTCAGGTAAGT‑G‑‑A‑G‑GTGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAAc         <  2:2516130/96‑1 (MQ=255)
                                                                                          aTCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTg    >  1:440699/1‑100 (MQ=255)
                                                                                          aTCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTg    >  1:441001/1‑100 (MQ=255)
                                                                                          aTCAGGT‑ACT‑G‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTgg   >  2:2580395/1‑100 (MQ=255)
                                                                                           tCAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTgg   >  1:2308411/1‑100 (MQ=255)
                                                                                            cAGGT‑ACTGG‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGt  >  2:1154040/1‑100 (MQ=255)
                                                                                                t‑ACT‑G‑‑GTG‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTgg   <  1:1112802/95‑1 (MQ=255)
                                                                                                     |                                                                                                
TCTGTTTATCACGGATTGTATGAGACATGGCAACACCTGGTTAACAAGAATATGAAAAATCATAGCACTATTAATCTACTGGGGGGTAGTATCAGGT‑ACT‑G‑‑G‑G‑GGGG‑AGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGT  >  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: ‑