BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,555	(G)_5→6	coding (55/618 nt)	b3646 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,550	1	.	G	98.2%	1238.8 / 11.0	399	coding (50/618 nt)	b3646	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (2/5); new base G (176/216); total (178/221)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 451 total aligned reads displayed.

CACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑GGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGT  >  NC_000913/3816458‑3816637

 acaAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑g                                                                                         >  2:7751962/1‑93 (MQ=255)

  caAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑gg                                                                                        <  2:1303944/93‑1 (MQ=255)

  aaaaGCGCTACCCTGCCAGACAGATATTTAGGGAGAAAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACAGCCTAAGCCATGAC‑gg                                                                                        <  2:4602435/92‑1 (MQ=255)

   aaaTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggg                                                                                       >  2:702567/1‑93 (MQ=255)

    aaTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggg‑g                                                                                     >  1:4082618/1‑93 (MQ=255)

    aaTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggg‑g                                                                                     >  2:6417148/1‑93 (MQ=255)

     aTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggg‑gg                                                                                    >  2:5933193/1‑93 (MQ=255)

     aTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggg‑gg                                                                                    <  1:4954094/93‑1 (MQ=255)

      tCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggg‑ggg                                                                                   <  1:7358966/93‑2 (MQ=255)

        gCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggg‑gggcg                                                                                 >  2:42975/1‑90 (MQ=255)

         cTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggg‑gggcgc                                                                                >  2:228716/1‑89 (MQ=255)

         cTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggg‑gggcgc                                                                                <  1:5443994/93‑5 (MQ=255)

         cTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggg‑gggcgc                                                                                <  1:4323729/93‑5 (MQ=255)

            cccTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTgg                                                                              <  2:7090704/93‑1 (MQ=255)

              cTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTggcg                                                                            >  2:575595/1‑93 (MQ=255)

                 ccAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGcc                                                                         >  2:5255953/1‑93 (MQ=255)

                   agacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCgg                                                                       >  2:5064486/1‑93 (MQ=255)

                   agacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCgg                                                                       <  2:7404692/93‑1 (MQ=255)

                    gacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGga                                                                      <  1:6881796/93‑1 (MQ=255)

                      cagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGgacg                                                                    >  1:6543730/1‑93 (MQ=255)

                        gaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGgacgac                                                                  >  1:4479204/1‑93 (MQ=255)

                        gaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGgacgac                                                                  <  1:5702180/93‑1 (MQ=255)

                        gaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGgacgac                                                                  <  1:886321/93‑1 (MQ=255)

                          tttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACgg                                                                >  2:5173097/1‑93 (MQ=255)

                          tttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACgg                                                                <  2:5427493/93‑1 (MQ=255)

                           ttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGgc                                                               >  1:7818501/1‑93 (MQ=255)

                           ttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGgc                                                               <  1:1081506/93‑1 (MQ=255)

                            tttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGgcg                                                              <  2:1245396/93‑1 (MQ=255)

                            tttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGgcg                                                              >  1:1611646/1‑93 (MQ=255)

                              tAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTCCTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCa                                                            <  2:6094112/93‑1 (MQ=255)

                              tAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCa                                                            <  1:261825/93‑1 (MQ=255)

                               aGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCag                                                           >  2:4891396/1‑92 (MQ=255)

                               aGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCAt                                                           <  2:1762216/93‑1 (MQ=255)

                                gggAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATg                                                          <  1:3867931/93‑1 (MQ=255)

                                gggAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATg                                                          <  1:3983770/93‑1 (MQ=255)

                                  gagagaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGa                                                        >  2:1092282/1‑93 (MQ=255)

                                  gagagaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGa                                                        >  1:3043705/1‑93 (MQ=255)

                                   agagaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGAt                                                       <  2:675492/93‑1 (MQ=255)

                                   agagaACCATGCTGTTACACATTGTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGAt                                                       <  1:6830600/93‑1 (MQ=255)

                                      gagccATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑GGGTGGCGCTGGCGGCCGGACGACGGCGCATGGATACa                                                    <  1:1254668/90‑1 (MQ=255)

                                      gaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACa                                                    <  1:3632757/93‑1 (MQ=255)

                                      gaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACa                                                    >  2:5887273/1‑93 (MQ=255)

                                      gaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACa                                                    >  1:7623000/1‑93 (MQ=255)

                                        aCCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACAtt                                                  <  2:6222851/93‑1 (MQ=255)

                                        aCCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACAtt                                                  >  1:5291430/1‑93 (MQ=255)

                                        aCCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACAtt                                                  <  2:2788113/93‑1 (MQ=255)

                                         ccATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACAttt                                                 >  1:6387793/1‑93 (MQ=255)

                                         ccAGGCTGTTCCACATTTTGTATTTGGTTGGCTTTACTGCCGAGGCAATGACGGGG‑GGCGCTGGCCGCCGGACGACGTCGCATGGATACAttt                                                 <  2:7638500/93‑1 (MQ=255)

                                          cATGCTGTTACACATTATGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTg                                                <  1:6198706/93‑1 (MQ=255)

                                            tGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGc                                              <  2:2392693/93‑1 (MQ=255)

                                            tGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGc                                              <  1:3437125/93‑1 (MQ=255)

                                               tGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTa                                           >  2:3511909/1‑93 (MQ=255)

                                               tGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTa                                           >  2:7470827/1‑93 (MQ=255)

                                                  tACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAatt                                        >  1:2260248/1‑93 (MQ=255)

                                                  tACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAatt                                        <  1:6753483/93‑1 (MQ=255)

                                                   acaaATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCGGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAatta                                       <  1:101228/93‑1 (MQ=255)

                                                    cacaTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattat                                      <  2:4401872/93‑1 (MQ=255)

                                                     acaTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGGCGGCGCAGGGATACATTTTGCGTTattatt                                     >  2:3679258/1‑93 (MQ=255)

                                                     acaTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattatt                                     >  1:1906099/1‑93 (MQ=255)

                                                      caTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattatta                                    <  1:2382505/93‑1 (MQ=255)

                                                      caTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattatta                                    >  1:7483713/1‑93 (MQ=255)

                                                       aTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattat                                   >  1:3243835/1‑93 (MQ=255)

                                                       aTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattat                                   <  2:3357780/93‑1 (MQ=255)

                                                       aTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattat                                   >  2:3662675/1‑93 (MQ=255)

                                                        ttttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattatt                                  <  1:2542961/93‑1 (MQ=255)

                                                        ttttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattatt                                  <  1:2818709/93‑1 (MQ=255)

                                                         tttGTATTTGGTTGGCATTACTGCCGAGGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTg                                 <  1:3217985/93‑1 (MQ=255)

                                                          ttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGc                                >  2:6526317/1‑93 (MQ=255)

                                                           tGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCg                               >  2:1147242/1‑93 (MQ=255)

                                                            gTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGa                              <  1:1457281/93‑1 (MQ=255)

                                                             tATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGAc                             <  1:4596880/93‑1 (MQ=255)

                                                             tATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGAc                             <  2:4100396/93‑1 (MQ=255)

                                                             tATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGAc                             <  2:3559273/93‑1 (MQ=255)

                                                             tATTTGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGAc                             >  2:2128611/1‑93 (MQ=255)

                                                               tttGGTTGGCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTCATTATTATTGCGACCg                           >  1:5630322/1‑93 (MQ=255)

                                                                 tggttggCATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGcc                         <  1:135799/93‑1 (MQ=255)

                                                                        cATTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCaa                  <  2:2568962/93‑1 (MQ=255)

                                                                         aTTACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAt                 >  2:2467485/1‑93 (MQ=255)

                                                                          ttACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAtt                <  2:2968753/93‑1 (MQ=255)

                                                                          ttACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAtt                >  2:2728737/1‑93 (MQ=255)

                                                                          ttACTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAtt                >  1:7699774/1‑93 (MQ=255)

                                                                            aCTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTgg              >  1:3148984/1‑93 (MQ=255)

                                                                             cTGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGc             <  1:3466818/93‑1 (MQ=255)

                                                                              tGCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCCCCACCGCAATTGGCg            >  2:3249421/1‑93 (MQ=255)

                                                                               gCCGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCgg           <  2:5619974/93‑1 (MQ=255)

                                                                                ccGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCggg          >  1:6129638/1‑92 (MQ=255)

                                                                                ccGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGa          <  1:2049017/93‑1 (MQ=255)

                                                                                 cGAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAg         <  1:921108/93‑1 (MQ=255)

                                                                                 cGAAGCAATGACGGGG‑GGAGCTGGCCGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAg         <  1:4838713/93‑1 (MQ=255)

                                                                                  gAAGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAgg        >  1:2192976/1‑93 (MQ=255)

                                                                                  gAAGCCATGACGGGG‑GGCGCGGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATGATTGCGACCGCCACCGCAATTGGCGGAgg        <  2:3871166/93‑1 (MQ=255)

                                                                                  gAAGCCATGAC‑GGGAGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAgg        <  2:1012333/93‑1 (MQ=255)

                                                                                   aaGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAggg       <  2:1641571/93‑1 (MQ=255)

                                                                                    aGCCATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGt      <  1:702288/93‑1 (MQ=255)

                                                                                      ccATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCa    >  2:7290717/1‑93 (MQ=255)

                                                                                       cATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAg   <  1:6315365/93‑1 (MQ=255)

                                                                                       cATGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAg   <  2:7896659/93‑1 (MQ=255)

                                                                                        aTGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGt  >  2:1378293/1‑93 (MQ=255)

                                                                                        aTGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGt  <  1:5969720/93‑1 (MQ=255)

                                                                                        aTGACGGGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGt  <  1:3710847/93‑1 (MQ=255)

CACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑GGG‑GGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGT  >  NC_000913/3816458‑3816637

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