BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,325,759	A→G	G135G (GGT→GGC)	b2225 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,325,759	0	A	G	98.6%	442.4 / ‑2.5	140	G135G (GGT→GGC)	b2225	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (77/61); total (79/61)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

Only 100 of 144 total aligned reads displayed.

GCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTACCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGCCTTCAGTGAAAATATA  >  NC_000913/2325669‑2325848

gCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGcc                                                                                         >  2:3947933/1‑93 (MQ=255)

 ccATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGccc                                                                                        >  1:292234/1‑93 (MQ=255)

 ccATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGccc                                                                                        >  1:3328670/1‑93 (MQ=255)

   aTCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCgg                                                                                      >  1:1176834/1‑93 (MQ=255)

   aTCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCgg                                                                                      >  2:1131007/1‑93 (MQ=255)

      gggCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGt                                                                                   >  2:1582137/1‑93 (MQ=255)

      gggCGTAACAACGTGAAGGTCGAGCTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGt                                                                                   <  1:3717632/93‑1 (MQ=255)

       ggCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGtt                                                                                  >  2:2432115/1‑93 (MQ=255)

       ggCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGtt                                                                                  >  1:4311131/1‑93 (MQ=255)

          ggaacaccGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAg                                                                               <  1:3827355/91‑1 (MQ=255)

          gTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAg                                                                               >  2:1332922/1‑93 (MQ=255)

          gTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAg                                                                               >  1:4523009/1‑93 (MQ=255)

          gTAACAACGTGAAGGTCGAGGGCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCCCGAATCGTTCCTGTGCCCGGCGTTGAg                                                                               >  2:3398118/1‑93 (MQ=255)

             acaacGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTaa                                                                            <  2:2701796/93‑1 (MQ=255)

              caacGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTaaa                                                                           >  2:4621441/1‑93 (MQ=255)

               aacGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTaaaa                                                                          >  1:1701566/1‑93 (MQ=255)

               aacGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTaaaa                                                                          <  1:905575/93‑1 (MQ=255)

                 cGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACt                                                                        >  2:505006/1‑93 (MQ=255)

                   tGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGc                                                                      >  1:455468/1‑93 (MQ=255)

                      aggtcgaggtcgGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCAtt                                                                   >  1:635623/1‑93 (MQ=255)

                      aggtcgaggtcgGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCAtt                                                                   <  2:726627/93‑1 (MQ=255)

                      agctcgagctcgGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCAtt                                                                   <  1:6261/90‑1 (MQ=255)

                        gtcgaggtcgGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCAtttt                                                                 >  1:587744/1‑93 (MQ=255)

                          cgaggtcgGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCAtttttt                                                               >  1:1963190/1‑93 (MQ=255)

                           gaggtcgGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTg                                                              <  1:1436114/93‑1 (MQ=255)

                            aggtcgTTATTGTCCGTACCCAACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGt                                                             <  2:2111834/93‑1 (MQ=255)

                            aggtcgGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGt                                                             <  2:2837362/93‑1 (MQ=255)

                             ggtcgGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTc                                                            >  1:2808655/1‑93 (MQ=255)

                             gctcgGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTc                                                            <  2:3640669/91‑1 (MQ=255)

                              gtcgGTTTTGCCCGTAACCCAGGAGAAACCCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCg                                                           <  2:1770822/93‑1 (MQ=255)

                              gtcgGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCg                                                           <  2:2934168/93‑1 (MQ=255)

                              gtcgGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCg                                                           <  1:4195950/93‑1 (MQ=255)

                              gtcgGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCg                                                           <  1:4184414/93‑1 (MQ=255)

                               tcgGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCgg                                                          >  1:3479918/1‑93 (MQ=255)

                                 ggTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCt                                                        >  2:1728182/1‑93 (MQ=255)

                                  gTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCtt                                                       <  2:311148/93‑1 (MQ=255)

                                  gTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCtt                                                       <  1:239867/93‑1 (MQ=255)

                                    aTTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTg                                                     <  1:2499415/93‑1 (MQ=255)

                                    aTTGCCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTg                                                     <  1:4260552/93‑1 (MQ=255)

                                     ttGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGc                                                    >  1:3292232/1‑93 (MQ=255)

                                      tGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGcc                                                   >  1:2192227/1‑93 (MQ=255)

                                      tGTCCGTACCCCACGAGAGACCCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGcc                                                   <  2:4513423/93‑1 (MQ=255)

                                       gTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGccc                                                  <  2:4010992/93‑1 (MQ=255)

                                          cGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATc                                               <  1:3846693/93‑1 (MQ=255)

                                            tATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCg                                             <  2:2763844/93‑1 (MQ=255)

                                              tCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCggg                                           >  1:860267/1‑93 (MQ=255)

                                              tCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCggg                                           <  1:2817906/93‑1 (MQ=255)

                                               cccACGAGAGAGCCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGc                                          <  1:4012399/93‑1 (MQ=255)

                                               cccACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGc                                          <  2:2079044/93‑1 (MQ=255)

                                                ccACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCCTCCGGGCt                                         >  1:1016608/1‑93 (MQ=255)

                                                ccACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCt                                         >  1:1477467/1‑93 (MQ=255)

                                                 cACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCtg                                        >  1:4351995/1‑93 (MQ=255)

                                                 cACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAAATGCATTTTTTGTCGGCTTTGCCCATCCGGGCtg                                        >  1:1327180/1‑93 (MQ=255)

                                                    gagagaACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCtgatg                                     >  2:3886884/1‑93 (MQ=255)

                                                     ggagaACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCtgatga                                    <  2:2888755/92‑1 (MQ=255)

                                                      gagaACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGAc                                   >  1:639546/1‑93 (MQ=255)

                                                       agaACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGAcc                                  >  1:1673829/1‑93 (MQ=255)

                                                       agaACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGAcc                                  <  2:2648286/93‑1 (MQ=255)

                                                       agaACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGAcc                                  >  1:137877/1‑93 (MQ=255)

                                                       agaACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGAcc                                  >  1:2522850/1‑93 (MQ=255)

                                                        gaACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCt                                 >  1:3914778/1‑93 (MQ=255)

                                                        gaACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCt                                 >  1:4495278/1‑93 (MQ=255)

                                                        gaACCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCt                                 >  2:498058/1‑93 (MQ=255)

                                                        gaACCAGCCGTAAACGTGCACGAATCGTTCCTGTACCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCt                                 >  2:237998/1‑93 (MQ=255)

                                                         aaCCAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTg                                >  2:1427628/1‑93 (MQ=255)

                                                            cacgcGTAAACGTGCACGAATCGTTCCTGTGCCCGGCATTGAGTAAAACTGCATTATTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCAc                             <  2:598775/89‑1 (MQ=255)

                                                            cAGCCGTAAACGTGCACTAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCAc                             <  2:1967745/93‑1 (MQ=255)

                                                            cAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCAc                             <  2:1355316/93‑1 (MQ=255)

                                                            cAGCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCAc                             <  1:3735706/93‑1 (MQ=255)

                                                              gCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGc                           >  2:2349591/1‑93 (MQ=255)

                                                              gCCGTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGc                           <  2:2528603/93‑1 (MQ=255)

                                                                 gTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGt                        >  2:4355942/1‑93 (MQ=255)

                                                                 gTAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGt                        <  2:3026621/93‑1 (MQ=255)

                                                                  tAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGtt                       >  1:391536/1‑93 (MQ=255)

                                                                  tAAACGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGtt                       <  1:259739/93‑1 (MQ=255)

                                                                    aaCGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTtgc                     >  1:1201895/1‑93 (MQ=255)

                                                                    aaCGTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTtgc                     <  1:1264477/93‑1 (MQ=255)

                                                                    aaCGTGCACGAATCGTTCCTATGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTtgc                     >  1:306941/1‑93 (MQ=255)

                                                                       gTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTtgctgc                  >  2:3902613/1‑93 (MQ=255)

                                                                       gTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTtgctgc                  >  1:202602/1‑93 (MQ=255)

                                                                       gTGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACATGCACGCTGTtgctgc                  >  1:3826565/1‑93 (MQ=255)

                                                                        tGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGcc                 >  1:3045695/1‑93 (MQ=255)

                                                                        tGCACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGcc                 <  2:3636720/93‑1 (MQ=255)

                                                                          cACGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGCCtt               <  2:3311769/93‑1 (MQ=255)

                                                                            cGAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGCCTTCa             <  1:4001358/93‑1 (MQ=255)

                                                                             gAATCGTTCCTGTGCCCGTCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGCCTTCAg            <  2:4633189/93‑1 (MQ=255)

                                                                             gAATCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGCCTTCAg            <  1:2992200/93‑1 (MQ=255)

                                                                                tCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGCCTTCAGTGa         <  2:2882173/93‑1 (MQ=255)

                                                                                tCGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGCCTTCAGTGa         <  1:4555318/93‑1 (MQ=255)

                                                                                 cGTTCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGCCTTCAGTGaa        <  2:1157671/93‑1 (MQ=255)

                                                                                   ttCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGCCTTCAGTGaaaa      >  2:4561865/1‑93 (MQ=255)

                                                                                   ttCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGCCTTCAGTGaaaa      <  1:2660161/93‑1 (MQ=255)

                                                                                   gtCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCGTTGCCCATCCGGGCTGAGGACCTGCACGCTGTTGCTGCCTTCAGTGaaaa      <  1:2660160/92‑1 (MQ=255)

                                                                                    tCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGCCTTCAGTGAAAat     <  2:2575631/93‑1 (MQ=255)

                                                                                    tCCTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGCCTTCAGTGAAAat     <  2:649467/93‑1 (MQ=255)

                                                                                     ccTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGCCTTCAGTGAAAata    <  2:3781628/93‑1 (MQ=255)

                                                                                     ccTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGCCTTCAGTGAAAata    >  1:1160693/1‑93 (MQ=255)

                                                                                      cTGTGCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGCCTTCAGTGAAAatat   >  1:143256/1‑93 (MQ=255)

                                                                                       tgtgCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGCCTTCAGTGAAAAtata  >  1:543886/1‑93 (MQ=255)

                                                                                       tgtgCCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGCCTTCAGTGAAAAtata  >  2:596963/1‑93 (MQ=255)

GCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAGCCGTAAACGTGCACGAATCGTTCCTGTACCCGGCGTTGAGTAAAACTGCATTTTTTGTCGGCTTTGCCCATCCGGGCTGATGACCTGCACGCTGTTGCTGCCTTCAGTGAAAATATA  >  NC_000913/2325669‑2325848

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