BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,989,898	A→G	Q378R (CAG→CGG)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,989,898	0	A	G	99.8%	1450.8 / NA	449	Q378R (CAG→CGG)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (197/251); total (197/252)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 466 total aligned reads displayed.

CGCGCAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCG  >  NC_000913/3989807‑3989987

  cgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGt                                                                                        <  1:3361582/93‑1 (MQ=255)

     aGTGTTGATCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCg                                                                                     <  1:2990695/93‑1 (MQ=255)

     aGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCg                                                                                     >  2:8561258/1‑93 (MQ=255)

     aGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCg                                                                                     >  2:3075249/1‑93 (MQ=255)

      gtgtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATTCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGt                                                                                    <  2:2182391/93‑1 (MQ=255)

      gtgtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGt                                                                                    <  2:3378798/93‑1 (MQ=255)

       tgtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGAGTg                                                                                   >  2:7586702/1‑93 (MQ=255)

        gtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGa                                                                                  <  2:7869975/93‑1 (MQ=255)

        gtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGa                                                                                  <  2:1587607/93‑1 (MQ=255)

        gtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGa                                                                                  <  1:3074094/93‑1 (MQ=255)

        gtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGa                                                                                  <  1:2879319/93‑1 (MQ=255)

         ttGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAg                                                                                 <  2:10018706/93‑1 (MQ=255)

         ttGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAg                                                                                 >  1:8437336/1‑93 (MQ=255)

          tGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAgg                                                                                <  1:6754671/93‑1 (MQ=255)

            aGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGgcg                                                                              >  1:5680527/1‑93 (MQ=255)

              ccAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGgcgcc                                                                            >  1:7355520/1‑92 (MQ=255)

              ccAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGca                                                                            >  1:4369100/1‑93 (MQ=255)

              ccAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGca                                                                            >  2:4021322/1‑93 (MQ=255)

              ccAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGca                                                                            >  2:6910736/1‑93 (MQ=255)

               cAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGcac                                                                           <  2:7479257/93‑1 (MQ=255)

                aGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGcaca                                                                          >  1:8682518/1‑93 (MQ=255)

                aGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGcaca                                                                          >  2:8868116/1‑93 (MQ=255)

                  ttAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAAc                                                                        >  1:7258053/1‑93 (MQ=255)

                    aGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGa                                                                      <  1:6891370/93‑1 (MQ=255)

                      tGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGt                                                                    <  1:2687607/93‑1 (MQ=255)

                      tGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGt                                                                    <  2:2355695/93‑1 (MQ=255)

                      tGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGt                                                                    >  2:4268171/1‑93 (MQ=255)

                         gCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGc                                                                 <  1:220870/93‑1 (MQ=255)

                             gTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGa                                                             >  2:675133/1‑93 (MQ=255)

                             gTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGa                                                             <  2:9326995/93‑1 (MQ=255)

                             gTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGa                                                             <  2:4783747/93‑1 (MQ=255)

                              tGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAc                                                            >  1:9908345/1‑93 (MQ=255)

                              tGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAc                                                            <  2:4482414/93‑1 (MQ=255)

                               ggggTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAcg                                                           <  1:9683704/93‑1 (MQ=255)

                               ggggTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAcg                                                           <  1:6634680/93‑1 (MQ=255)

                                gggTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAcgc                                                          <  2:2502055/93‑1 (MQ=255)

                                 ggTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAcgcg                                                         <  2:8413728/93‑1 (MQ=255)

                                  gTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCga                                                        <  1:7709697/93‑1 (MQ=255)

                                  gTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCga                                                        <  1:1879808/93‑1 (MQ=255)

                                  gTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCga                                                        <  2:3537372/93‑1 (MQ=255)

                                   ttGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgat                                                       >  2:1112833/1‑93 (MQ=255)

                                     gggACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatga                                                     <  2:267660/93‑1 (MQ=255)

                                     gggACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatga                                                     >  1:1298035/1‑93 (MQ=255)

                                     gggACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatga                                                     <  2:6210715/93‑1 (MQ=255)

                                     gggACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatga                                                     <  2:9697057/93‑1 (MQ=255)

                                        aCGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGc                                                  <  2:3198453/93‑1 (MQ=255)

                                         cGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCa                                                 <  1:10002552/93‑1 (MQ=255)

                                          gAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCag                                                >  1:6389997/1‑93 (MQ=255)

                                           aaGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCaga                                               <  2:6549062/93‑1 (MQ=255)

                                             gCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAATATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGc                                             <  1:8865029/93‑1 (MQ=255)

                                             gCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGc                                             <  2:496767/93‑1 (MQ=255)

                                             gCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGc                                             >  2:2779003/1‑93 (MQ=255)

                                             gCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGc                                             <  1:3152542/93‑1 (MQ=255)

                                             gCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGc                                             <  2:1306218/93‑1 (MQ=255)

                                             gCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGc                                             <  2:5363977/93‑1 (MQ=255)

                                              cTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCt                                            >  1:9525059/1‑93 (MQ=255)

                                               tCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGGAGAGCTa                                           >  2:9581804/1‑93 (MQ=255)

                                               tCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGGGCGTGAGGCGCACAACGAGGTGCTCGACGCGCTGATGCCGCGCTa                                           >  1:4771321/1‑93 (MQ=255)

                                                     ggCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGGGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTaa                                     >  1:1241701/1‑93 (MQ=255)

                                                     ggCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTaa                                     >  1:613980/1‑93 (MQ=255)

                                                     ggCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTaa                                     <  2:1856943/93‑1 (MQ=255)

                                                       cAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATc                                   <  2:8762401/93‑1 (MQ=255)

                                                       cAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATc                                   <  1:2490788/93‑1 (MQ=255)

                                                        aaTGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCt                                  >  1:3826968/1‑93 (MQ=255)

                                                          tGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGa                                >  2:2048155/1‑93 (MQ=255)

                                                           gCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGAt                               >  2:3623631/1‑93 (MQ=255)

                                                            cTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATc                              >  1:9029419/1‑93 (MQ=255)

                                                            cTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATc                              >  2:8263503/1‑93 (MQ=255)

                                                             tCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATcc                             >  2:5774104/1‑93 (MQ=255)

                                                              cGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCg                            >  2:3824884/1‑93 (MQ=255)

                                                                aTAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCt                          >  1:7547768/1‑93 (MQ=255)

                                                                aTAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCt                          >  2:8028700/1‑93 (MQ=255)

                                                                aTAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCt                          <  1:5393798/93‑1 (MQ=255)

                                                                  aaCCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATTCAGAGCTACCGTAATCTGATCCGCttt                        <  1:8174887/93‑1 (MQ=255)

                                                                    ccGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTgc                      <  2:9469233/93‑1 (MQ=255)

                                                                     cgcgATAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTgcg                     <  2:4575804/93‑1 (MQ=255)

                                                                       cgcTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGcgcg                   >  2:6374819/1‑93 (MQ=255)

                                                                        gcTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGt                  >  1:176833/1‑93 (MQ=255)

                                                                        gcTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGt                  >  1:4261588/1‑93 (MQ=255)

                                                                        gcTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGt                  <  2:295136/93‑1 (MQ=255)

                                                                         cTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTc                 >  1:8564234/1‑93 (MQ=255)

                                                                           aaCTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGc               >  2:1441019/1‑93 (MQ=255)

                                                                           aaCTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGc               >  2:913062/1‑93 (MQ=255)

                                                                             cTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCaa             <  1:9311851/93‑1 (MQ=255)

                                                                              tGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAAt            <  2:7311608/93‑1 (MQ=255)

                                                                              tGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAAt            <  2:5563840/93‑1 (MQ=255)

                                                                                  aGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAAcc        >  1:5083646/1‑93 (MQ=255)

                                                                                  aGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAAcc        <  2:707529/93‑1 (MQ=255)

                                                                                  aGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAAcc        <  1:626986/93‑1 (MQ=255)

                                                                                  aGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAAcc        >  1:8332910/1‑93 (MQ=255)

                                                                                   gATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCt       <  1:4495445/93‑1 (MQ=255)

                                                                                   gATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCt       >  1:3591185/1‑93 (MQ=255)

                                                                                    aTTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCtt      <  1:5496676/93‑1 (MQ=255)

                                                                                    aTTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCtt      <  1:6993464/93‑1 (MQ=255)

                                                                                     ttGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTa     <  2:980628/93‑1 (MQ=255)

                                                                                     ttGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTa     <  2:6004602/93‑1 (MQ=255)

                                                                                      tGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAg    <  1:3413172/93‑1 (MQ=255)

                                                                                      tGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAg    >  1:2054616/1‑93 (MQ=255)

                                                                                       gATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGc   <  2:3931685/93‑1 (MQ=255)

                                                                                        aTCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCg  <  1:4932852/93‑1 (MQ=255)

CGCGCAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCG  >  NC_000913/3989807‑3989987

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 29 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`