BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,029,120	+C	coding (73/1191 nt)	b0967 ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,029,116	1	.	C	98.6%	659.5 / 0.7	210	coding (77/1191 nt)	b0967	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base . (2/1); new base C (94/113); total (96/114)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 237 total aligned reads displayed.

CGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑CCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTT  >  NC_000913/1029025‑1029203

cGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑c                                                                                        <  1:5006943/93‑1 (MQ=255)

 gAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑cc                                                                                       >  2:6524486/1‑93 (MQ=255)

  aaTAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑ccc                                                                                      <  2:2469545/93‑1 (MQ=255)

  aaTAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑ccc                                                                                      >  2:6256937/1‑93 (MQ=255)

    tAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑ccccc                                                                                    >  2:466141/1‑92 (MQ=255)

    tAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑ccccc                                                                                    >  2:4108834/1‑92 (MQ=255)

      aGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑cccccgg                                                                                  >  2:1025430/1‑90 (MQ=255)

         gccgcGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑CCCCCGgaaa                                                                               >  2:2550381/1‑93 (MQ=255)

            gcgTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGa                                                                             <  2:3717398/93‑1 (MQ=255)

              gTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGAcc                                                                           >  2:2918990/1‑93 (MQ=255)

              gTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGAcc                                                                           <  2:3588786/93‑1 (MQ=255)

              gTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGAcc                                                                           >  2:1810614/1‑93 (MQ=255)

               tGCTAACCATTTTTCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGAccc                                                                          >  1:5796666/1‑93 (MQ=255)

               tGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGAccc                                                                          <  1:6268853/93‑1 (MQ=255)

               tGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGAccc                                                                          >  1:4798979/1‑93 (MQ=255)

                gCTAACCATTTTCCCTGATGATCAACAATATCGATTGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCa                                                                         >  1:5293739/1‑93 (MQ=255)

                gCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCa                                                                         <  2:4199301/93‑1 (MQ=255)

                gCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCa                                                                         <  2:2352620/93‑1 (MQ=255)

                gCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCa                                                                         <  1:3583691/93‑1 (MQ=255)

                 cTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCAc                                                                        <  2:6421336/93‑1 (MQ=255)

                    aCCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGa                                                                     >  1:2993926/1‑93 (MQ=255)

                    aCCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGa                                                                     <  1:4711050/93‑1 (MQ=255)

                     ccATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGAt                                                                    <  2:4841386/93‑1 (MQ=255)

                     ccATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGAt                                                                    >  1:6717906/1‑93 (MQ=255)

                     ccATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGAt                                                                    <  1:272518/93‑1 (MQ=255)

                       aTTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATgg                                                                  >  1:5486278/1‑93 (MQ=255)

                       aTTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATgg                                                                  >  2:3713183/1‑93 (MQ=255)

                        ttttCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGc                                                                 >  2:3943418/1‑93 (MQ=255)

                             ccTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGcgacga                                                            <  2:2659844/93‑1 (MQ=255)

                             ccTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGcgacga                                                            <  1:996802/93‑1 (MQ=255)

                              cTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGaa                                                           >  2:1691782/1‑93 (MQ=255)

                              cTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGaa                                                           >  2:2007140/1‑93 (MQ=255)

                                tatgatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGt                                                         <  2:2515815/92‑1 (MQ=255)

                                 atgatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGACGGCGACGAAGTa                                                        >  2:3253072/1‑93 (MQ=255)

                                 atgatAAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTTCATGCGGGCAACGGCCCCCGGCAAAGACCCACGGATGGCGACGAAGTa                                                        >  1:1752411/1‑93 (MQ=255)

                                  tgatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTaa                                                       <  2:6100346/93‑1 (MQ=255)

                                  tgatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTaa                                                       >  2:3424793/1‑93 (MQ=255)

                                   gatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAAt                                                      <  1:1013709/93‑1 (MQ=255)

                                   gatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAAt                                                      <  2:1159551/93‑1 (MQ=255)

                                     tCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGa                                                    >  1:838323/1‑93 (MQ=255)

                                     tCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGa                                                    >  1:5875350/1‑93 (MQ=255)

                                      caacaaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGAt                                                   <  2:2211159/93‑1 (MQ=255)

                                       aacaaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGAtt                                                  >  2:2257188/1‑93 (MQ=255)

                                          aaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGAttttt                                               >  2:5479319/1‑93 (MQ=255)

                                           atatCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTc                                              >  1:1253588/1‑93 (MQ=255)

                                           atatCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTc                                              <  1:4930933/93‑1 (MQ=255)

                                           atatCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTc                                              >  1:6588523/1‑93 (MQ=255)

                                              tCGATGGTTTCGCCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTcgcg                                           <  1:1523196/93‑1 (MQ=255)

                                                gATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGcc                                         >  2:1748919/1‑93 (MQ=255)

                                                 aTGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGccc                                        <  2:6467736/93‑1 (MQ=255)

                                                 aTGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGccc                                        >  2:5644137/1‑93 (MQ=255)

                                                  tGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCt                                       <  2:1263810/93‑1 (MQ=255)

                                                  tGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCt                                       >  2:3487054/1‑93 (MQ=255)

                                                  tGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCt                                       <  1:6624745/93‑1 (MQ=255)

                                                  tGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCt                                       <  1:5416848/93‑1 (MQ=255)

                                                   ggTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCtt                                      <  2:5786330/93‑1 (MQ=255)

                                                   ggTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCtt                                      <  2:4270408/93‑1 (MQ=255)

                                                    gTTTCACCGAGTCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCttt                                     <  1:2510177/93‑1 (MQ=255)

                                                    gTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCttt                                     <  1:6132893/93‑1 (MQ=255)

                                                    gTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCttt                                     <  2:5547990/93‑1 (MQ=255)

                                                    gTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCttt                                     <  2:4389394/93‑1 (MQ=255)

                                                     tttCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTg                                    >  1:2264199/1‑93 (MQ=255)

                                                     tttCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTg                                    <  2:141328/93‑1 (MQ=255)

                                                      ttCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTgg                                   >  1:1120253/1‑93 (MQ=255)

                                                      ttCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTgg                                   <  2:1375759/93‑1 (MQ=255)

                                                       tCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGc                                  <  2:826227/93‑1 (MQ=255)

                                                       tCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGc                                  <  1:4028909/93‑1 (MQ=255)

                                                        cACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCt                                 <  1:335319/93‑1 (MQ=255)

                                                          ccGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTaa                               >  1:6317235/1‑93 (MQ=255)

                                                           cGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAac                              <  2:5035873/93‑1 (MQ=255)

                                                              ggctggctTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACt                           >  2:334915/1‑93 (MQ=255)

                                                                ctggctTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTaa                         <  1:1831941/93‑1 (MQ=255)

                                                                ctggctTTAACTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTaa                         <  2:5449521/93‑1 (MQ=255)

                                                                  ggctTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAAc                       <  2:919349/93‑1 (MQ=255)

                                                                  ggctTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAAc                       <  2:5237567/93‑1 (MQ=255)

                                                                  ggctTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAAc                       <  1:755932/93‑1 (MQ=255)

                                                                   gctTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACg                      >  2:2868768/1‑93 (MQ=255)

                                                                   gctTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACg                      >  2:3158979/1‑93 (MQ=255)

                                                                    ctTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGt                     >  2:3232517/1‑93 (MQ=255)

                                                                    ctTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGt                     <  2:1436852/93‑1 (MQ=255)

                                                                    ctTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGt                     <  1:3480842/93‑1 (MQ=255)

                                                                     tttACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTa                    <  2:3029556/93‑1 (MQ=255)

                                                                      ttACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTac                   <  1:5654226/93‑1 (MQ=255)

                                                                       tACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTaca                  <  1:6038832/93‑1 (MQ=255)

                                                                       tACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTaca                  >  2:103/1‑93 (MQ=255)

                                                                        aCCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTacac                 <  1:6406281/93‑1 (MQ=255)

                                                                        aCCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTacac                 <  2:5996149/93‑1 (MQ=255)

                                                                            tCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTTCACTCAt             <  1:302412/93‑1 (MQ=255)

                                                                             ccATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATa            <  2:6295395/93‑1 (MQ=255)

                                                                              cATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATaa           >  2:4454175/1‑93 (MQ=255)

                                                                              cATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATaa           <  2:6632675/93‑1 (MQ=255)

                                                                               aTGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAAt          >  1:4270334/1‑93 (MQ=255)

                                                                                tGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAAtt         <  2:5807896/93‑1 (MQ=255)

                                                                                tGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAAtt         >  1:2798584/1‑93 (MQ=255)

                                                                                  cGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTa       <  1:3630107/93‑1 (MQ=255)

                                                                                   gggCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTAc      <  2:2365833/93‑1 (MQ=255)

                                                                                    ggCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACt     <  2:4355602/93‑1 (MQ=255)

                                                                                      cAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATCATTTACttt   >  1:1431341/1‑93 (MQ=255)

                                                                                      cAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACttt   <  1:3261244/93‑1 (MQ=255)

                                                                                      cAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACttt   <  1:1910024/93‑1 (MQ=255)

CGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑CCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTT  >  NC_000913/1029025‑1029203

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 20 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`