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	99.1%	664.3 / ‑2.1	211	coding (77/1191 nt)	b0967	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base . (2/0); new base C (100/109); total (102/109)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 233 total aligned reads displayed.

GCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑CCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTT  >  NC_000913/1029024‑1029203

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

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

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

    aTAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑cccc                                                                                     <  2:5983550/93‑1 (MQ=255)

    aTAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑cccc                                                                                     >  1:1636117/1‑93 (MQ=255)

        gcgcCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑cccccgga                                                                                 >  2:3374604/1‑89 (MQ=255)

           ccgcGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCGATGCGGGCAACGGCCCCCGGaaaa                                                                               <  1:3847324/93‑1 (MQ=255)

             gcgTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGa                                                                             >  2:2226840/1‑93 (MQ=255)

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

             gcgTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGa                                                                             >  1:5103028/1‑93 (MQ=255)

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

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

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

                 gCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCa                                                                         >  2:3790194/1‑93 (MQ=255)

                  cTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCAc                                                                        >  1:66638/1‑93 (MQ=255)

                   tAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACg                                                                       <  2:1201665/93‑1 (MQ=255)

                   tAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACg                                                                       >  1:2071/1‑93 (MQ=255)

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

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

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

                      ccAGTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGAt                                                                    <  1:5692944/93‑1 (MQ=255)

                       cATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATg                                                                   <  2:3138160/93‑1 (MQ=255)

                        aTTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATgg                                                                  <  1:511211/93‑1 (MQ=255)

                         ttttCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGc                                                                 <  2:1993407/93‑1 (MQ=255)

                            tCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGcgac                                                              >  2:2848950/1‑93 (MQ=255)

                             cccTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGcgacg                                                             >  1:3240335/1‑93 (MQ=255)

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

                               cTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGaa                                                           <  2:1842246/93‑1 (MQ=255)

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

                                 gatgatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAga                                                         >  2:5596102/1‑92 (MQ=255)

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

                                     atCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATg                                                     <  2:675351/93‑1 (MQ=255)

                                     atCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATg                                                     <  1:4729494/93‑1 (MQ=255)

                                     atCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATg                                                     <  1:6165278/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                         acaaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGAttt                                                 <  2:3542193/93‑1 (MQ=255)

                                          caaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGAtttt                                                >  2:5729001/1‑93 (MQ=255)

                                          caaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGAtttt                                                <  2:4294524/93‑1 (MQ=255)

                                          caaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGAtttt                                                <  1:1698074/93‑1 (MQ=255)

                                           aaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGGAATGAttttt                                               >  2:422503/1‑93 (MQ=255)

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

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

                                              atCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTcgc                                            >  2:5764402/1‑93 (MQ=255)

                                               tCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTcgcg                                           >  1:1529900/1‑93 (MQ=255)

                                               tCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTcgcg                                           >  1:600345/1‑93 (MQ=255)

                                                cGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCgcgc                                          >  1:5586308/1‑93 (MQ=255)

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

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

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

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

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

                                                     gTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCttt                                     >  1:1484543/1‑93 (MQ=255)

                                                      tttCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTATTGATTTTTCGCGCCCTTTg                                    <  1:2544537/93‑1 (MQ=255)

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

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

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

                                                        tCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACTCACGGATGGCGACGAAGCAATGATTATTTGCGCCATTTTGc                                  >  1:3004241/1‑93 (MQ=255)

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

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

                                                         cACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCt                                 >  2:482004/1‑93 (MQ=255)

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

                                                             gAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAaca                             >  1:2230633/1‑93 (MQ=255)

                                                             gAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAaca                             <  1:242495/93‑1 (MQ=255)

                                                             gAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAaca                             <  2:3046310/93‑1 (MQ=255)

                                                              aGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAacac                            >  1:1380169/1‑93 (MQ=255)

                                                                gctggctTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTa                          <  2:3959805/93‑1 (MQ=255)

                                                                  tggctTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTaaa                        <  1:339127/93‑1 (MQ=255)

                                                                  tggctTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTaaa                        >  1:2146920/1‑93 (MQ=255)

                                                                  tggctTTACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTaaa                        <  2:4746934/93‑1 (MQ=255)

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

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

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

                                                                       ttACCTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTac                   >  1:1258699/1‑93 (MQ=255)

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

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

                                                                         tccTTCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTacac                 <  2:4848171/92‑1 (MQ=255)

                                                                            ttCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCa              >  1:2392810/1‑93 (MQ=255)

                                                                            ttCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCa              <  2:5090721/93‑1 (MQ=255)

                                                                            ttCCATGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCa              <  2:2513800/93‑1 (MQ=255)

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

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

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

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

                                                                                aTGCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAAt          <  2:2816337/93‑1 (MQ=255)

                                                                                aTGCGGCCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAAt          <  2:2244307/93‑1 (MQ=255)

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

                                                                                  gCGGGCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAAttt        <  2:1486622/93‑1 (MQ=255)

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

                                                                                     ggCAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACt     >  1:4769262/1‑93 (MQ=255)

                                                                                       cAACGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTCCttt   <  1:4068490/93‑1 (MQ=255)

                                                                                        aaCGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACtttt  <  2:5212012/93‑1 (MQ=255)

                                                                                        aaCGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACtttt  >  1:4276790/1‑93 (MQ=255)

                                                                                        aaCGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACtttt  <  1:5550637/93‑1 (MQ=255)

                                                                                        aaCGGCCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACtttt  <  2:3286014/93‑1 (MQ=255)

GCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑CCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTT  >  NC_000913/1029024‑1029203

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