BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,003,385	Δ1 bp	coding (186/681 nt)	ycbF →	putative chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,003,384	0	A	.	89.2%	311.2 / 35.2	96	coding (185/681 nt)	ycbF	putative chaperone
Reads supporting (aligned to +/- strand): ref base A (2/1); major base . (53/30); minor base G (7/3); total (62/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01

Only 100 of 117 total aligned reads displayed.

AGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAG  >  NC_000913/1003292‑1003471

agaATCTATATTTTTTTTACTGACTAATCCCAGTAAGAAATCCTGTTTTATAAAAAGAAAAACCACCCCCCCAACGCGTTGGGCCGGGGGGGa                                                                                         <  2:4390824/93‑1 (MQ=255)

agaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTACCAGTAAATTCAACCGCCCAACGCGTTGGGCCGGGGGGGa                                                                                         <  2:4071837/93‑1 (MQ=255)

agaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGa                                                                                         <  2:840264/93‑1 (MQ=255)

 gaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGag                                                                                        >  2:3904353/1‑92 (MQ=255)

 gaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGag                                                                                        >  1:2210184/1‑92 (MQ=255)

  aaTCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGagc                                                                                       >  2:624973/1‑91 (MQ=255)

   aTCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGagcg                                                                                      >  2:2227870/1‑90 (MQ=255)

   aTCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGagcg                                                                                      >  2:2227869/1‑90 (MQ=255)

    tCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGagcgt                                                                                     <  2:521342/93‑5 (MQ=255)

    tCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGagcgt                                                                                     >  1:1896593/1‑89 (MQ=255)

    tCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACACGTTGGGCCGGGGGGGagcgt                                                                                     >  1:1838751/1‑89 (MQ=255)

      gATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCg                                                                                  >  1:3890356/1‑93 (MQ=255)

       atatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGa                                                                                 >  2:14812/1‑93 (MQ=255)

       atatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGa                                                                                 >  1:1234074/1‑93 (MQ=255)

        tatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGAc                                                                                >  1:527285/1‑93 (MQ=255)

         atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAGCGTCGAc                                                                                >  1:2108807/1‑93 (MQ=255)

         atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACa                                                                               >  1:641135/1‑93 (MQ=255)

              ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGcc                                                                          <  1:4549799/93‑1 (MQ=255)

              ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGcc                                                                          <  1:1291063/93‑1 (MQ=255)

               tttACTGACTAATCCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCCCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCa                                                                         <  2:4633711/93‑1 (MQ=255)

               tttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCa                                                                         >  2:4502115/1‑93 (MQ=255)

               tttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTAGACAGTGCCa                                                                         <  1:743861/93‑1 (MQ=255)

                ttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAg                                                                        >  1:835147/1‑93 (MQ=255)

                ttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAg                                                                        >  1:549328/1‑93 (MQ=255)

                  aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCa                                                                      >  2:2015146/1‑93 (MQ=255)

                    tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCAca                                                                    >  2:2227924/1‑92 (MQ=255)

                     gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCACCCCCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCa                                                                   <  2:3065339/93‑1 (MQ=255)

                     gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCa                                                                   <  1:843046/93‑1 (MQ=255)

                      aCTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑GGCGTCGACAGTGCCAGCACCAt                                                                  >  1:4163862/1‑93 (MQ=255)

                      aCTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑CGCGTCGACAGTGCCAGCACCAt                                                                  >  1:3132801/1‑93 (MQ=255)

                       cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCAtt                                                                 >  1:908123/1‑93 (MQ=255)

                       cTAATACCAGTAAGAAATCCTGGCTTATAAACAGAAAAACCACCCCCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCAtt                                                                 <  2:3592322/93‑1 (MQ=255)

                         aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTAc                                                               <  1:3865866/93‑1 (MQ=255)

                         aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTAc                                                               <  2:3173458/93‑1 (MQ=255)

                          aTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCAACATTACt                                                              >  2:70453/1‑93 (MQ=255)

                           tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTg                                                             >  1:1971907/1‑93 (MQ=255)

                           tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTg                                                             >  2:8934/1‑93 (MQ=255)

                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTgg                                                            >  1:296070/1‑93 (MQ=255)

                            aCCAGCCAGGATTCCTGGTTAATAAACAGTAAAATCACCCCCCCATCGCG‑TGGGCCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTgg                                                            <  2:4523394/93‑1 (MQ=255)

                              cagtcagtAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGTGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGcc                                                          <  1:508147/93‑1 (MQ=255)

                              cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑CGCGTCGACAATGCCAGCACCATTACTGGcc                                                          >  2:4039171/1‑93 (MQ=255)

                                gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGc                                                        >  2:1854977/1‑93 (MQ=255)

                                gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGc                                                        >  1:1122424/1‑93 (MQ=255)

                                gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGc                                                        >  1:2353035/1‑93 (MQ=255)

                                  cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTc                                                      >  1:2898893/1‑93 (MQ=255)

                                   aGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                                     >  1:1076171/1‑93 (MQ=255)

                                   aGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                                     >  2:4405130/1‑93 (MQ=255)

                                    ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg                                                    >  2:862782/1‑93 (MQ=255)

                                    ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg                                                    <  2:4695555/93‑1 (MQ=255)

                                    ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCAAAGGCGTTGGGCCGGGGGGG‑GGGGTCGCCGGGGCCGGCACCATTACTGGCCGTTccc                                                    >  1:2087940/1‑92 (MQ=255)

                                     gAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGCGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg                                                    >  2:3331432/1‑93 (MQ=255)

                                      aaTCCTGGCTTATGAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCAATACTGGCCGCTCCGcc                                                  >  2:3610630/1‑93 (MQ=255)

                                       aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCATCACCATTACTGCCCGCTCCGCCa                                                 >  2:3072309/1‑93 (MQ=255)

                                        tCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa                                                 >  2:3588699/1‑93 (MQ=255)

                                          cTGTCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt                                              <  1:3159131/93‑1 (MQ=255)

                                          cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt                                              >  2:4168921/1‑93 (MQ=255)

                                           tGGCTTATAAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTa                                             <  1:2977030/93‑1 (MQ=255)

                                              tttattaACAGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtc                                          <  1:3638441/92‑1 (MQ=255)

                                              cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtc                                          >  1:1519189/1‑93 (MQ=255)

                                                tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtctc                                        >  2:4625280/1‑93 (MQ=255)

                                                tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtctc                                        >  2:1711093/1‑93 (MQ=255)

                                                  ttaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCt                                      >  2:3059475/1‑93 (MQ=255)

                                                   taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGACGCTCCGCCACTTATTCTCCTg                                     >  2:2131675/1‑93 (MQ=255)

                                                    aaCAGAAAAACCACCCGCCCACGGGGTGGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGa                                    <  2:308128/93‑1 (MQ=255)

                                                     acggTAAAATCAACCCCACAGCGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGaa                                   <  1:462325/90‑1 (MQ=255)

                                                     aaaGTAAAAACAACCCCCAAGGGGGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGaa                                   <  1:76831/91‑1 (MQ=255)

                                                      aaaaaaaaCCACCCCCCCACCGTGTGGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAg                                  <  2:3626799/89‑1 (MQ=255)

                                                        gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGGCGCTCCGCCACTTATTCTCCTGAAGc                                 >  2:3096139/1‑93 (MQ=255)

                                                        gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCCCTTATTCTCCTGAAGcc                                >  2:3707723/1‑93 (MQ=255)

                                                        gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                                >  2:1637215/1‑93 (MQ=255)

                                                        gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                                >  2:4497602/1‑93 (MQ=255)

                                                        gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                                >  1:130557/1‑93 (MQ=255)

                                                        gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                                >  2:119135/1‑93 (MQ=255)

                                                        aaaaaaCCACCCCCCCACGGCGTGGTGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                                <  1:3488402/91‑1 (MQ=255)

                                                         tAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGccc                               <  1:3271702/93‑1 (MQ=255)

                                                          aaaaTCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCg                              >  2:829792/1‑93 (MQ=255)

                                                           aaaTCAACCGCCAAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCgg                             <  2:3443634/93‑1 (MQ=255)

                                                              tCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTAc                          >  1:1692735/1‑93 (MQ=255)

                                                               cAACCGCCCACCGCGTTGGGCCGGGGGGGGTGCGTCGACAGTGCCAGGACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAACCCGGtcg                          >  2:4053562/1‑91 (MQ=255)

                                                                 aCCGCCCAACGCGTTGGGACGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGAc                       <  1:3686032/93‑1 (MQ=255)

                                                                    gCCCAACGCGTTGGGCCGGGGGGG‑GGCGTCGACAGTGCCAGCACCATTACTGGCCGGTCCGCACCTTATTATCTTGAATACCGGAACGACggg                    >  2:2275469/1‑93 (MQ=255)

                                                                    gCCCAACGCGTTGGGCCGGGGGGG‑GGCGGCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                    >  2:2275470/1‑93 (MQ=255)

                                                                    gCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCAGCCACTTATTTTCTTGAAGCCCGGTACGACAgg                    >  2:2275575/1‑93 (MQ=255)

                                                                      ccAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCa                  >  1:2751830/1‑93 (MQ=255)

                                                                        aaCGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACg                >  1:2279741/1‑93 (MQ=255)

                                                                        aaCGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACg                >  2:581819/1‑93 (MQ=255)

                                                                         aCGGGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGt               <  1:2123003/93‑1 (MQ=255)

                                                                         aCGCGTTGGGGCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACg                >  2:3993671/1‑93 (MQ=255)

                                                                          cgcgTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGtt              >  1:54507/1‑93 (MQ=255)

                                                                            gtgggggCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgc            <  1:254838/89‑1 (MQ=255)

                                                                             gTTGGGCCGGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgc            >  1:4641607/1‑93 (MQ=255)

                                                                             gTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcg           >  2:1876477/1‑93 (MQ=255)

                                                                              tggggCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcgc          <  1:2499126/91‑1 (MQ=255)

                                                                                  gCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGc      <  2:185730/93‑1 (MQ=255)

                                                                                    gggggggggAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCt     <  2:2303491/92‑1 (MQ=255)

                                                                                    cGGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCt     <  2:2303490/93‑1 (MQ=255)

                                                                                     cggggggGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTg    <  1:4055321/92‑1 (MQ=255)

                                                                                     cggggggGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTg    >  1:887658/2‑93 (MQ=255)

                                                                                       ggggggAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTgag  >  2:2217566/1‑93 (MQ=255)

                                                                                       ggggggAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTgag  >  1:873798/1‑93 (MQ=255)

AGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAG  >  NC_000913/1003292‑1003471

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