BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	474,195	Δ1 bp	5.2%	coding (793/1983 nt)	dnaG ←	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	474,195	0	C	.	5.2%	648.4 / 27.6	194	coding (793/1983 nt)	dnaG	DNA primase
Reads supporting (aligned to +/- strand): ref base C (92/91); new base . (2/8); total (94/100)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 201 total aligned reads displayed.

CTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCGATGATGCGC  >  NC_002947/474046‑474343

 tGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGGCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTcgacgacg                                                                                                                                                      >  1:75759/1‑150 (MQ=255)

       ccGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGAtc                                                                                                                                                <  2:1705069/150‑1 (MQ=255)

            aaaaCAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑tcgtc                                                                                                                                          <  2:452765/150‑1 (MQ=255)

            aaaaCAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑tcgtc                                                                                                                                          <  2:1872942/150‑1 (MQ=255)

             aaaCAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑tcgtcg                                                                                                                                         >  1:5067700/1‑150 (MQ=255)

              aacagaacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGA‑GATCATCGTCGTCGa                                                                                                                                        <  2:1121021/150‑1 (MQ=255)

              aacagaacacCACGCTGGGCACCACGCGCAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGCTGCCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGa                                                                                                                                        <  2:183463/150‑1 (MQ=255)

               acagaacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGA‑GATCATCGTCGTCGAg                                                                                                                                       <  2:836745/150‑1 (MQ=255)

               acagaacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGA‑GATCATCGTCGTCGAg                                                                                                                                       <  1:1384053/150‑1 (MQ=255)

                 agaacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGt                                                                                                                                     >  1:4240910/1‑150 (MQ=255)

                     cagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCg                                                                                                                                 >  1:4411949/1‑150 (MQ=255)

                      cgCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCgg                                                                                                                                <  2:3344235/149‑1 (MQ=255)

                      agCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCTCTCGACGA‑GATCATCGTCGTCGAgggctaca                                                                                                                                <  1:3743794/150‑7 (MQ=255)

                      agCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCgg                                                                                                                                <  2:3253168/150‑1 (MQ=255)

                      agCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCgg                                                                                                                                <  1:2765051/150‑1 (MQ=255)

                       gCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGt                                                                                                                               >  2:3184644/1‑150 (MQ=255)

                       gCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGACCTCGACGA‑GATCATCGTCGTCGAgggctacat                                                                                                                               <  1:829178/150‑8 (MQ=255)

                        cACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGtg                                                                                                                              >  2:4884159/1‑149 (MQ=255)

                              gggCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCtt                                                                                                                        <  1:3226285/150‑1 (MQ=255)

                              gggCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCtt                                                                                                                        <  2:750461/150‑1 (MQ=255)

                               ggCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTAGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTg                                                                                                                       <  2:385273/150‑1 (MQ=255)

                                     acGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGcc                                                                                                                 >  1:2546915/1‑150 (MQ=255)

                                       gcgGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTc                                                                                                               <  1:5357352/150‑1 (MQ=255)

                                       gcgGAACAGGCGCTGGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTc                                                                                                               <  2:2669925/150‑1 (MQ=255)

                                          gAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTa                                                                                                            <  1:2415548/150‑1 (MQ=255)

                                             cAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAg                                                                                                         <  1:5022037/150‑1 (MQ=255)

                                                gcgcTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGccc                                                                                                      >  2:4606428/1‑150 (MQ=255)

                                                gcgcTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGccc                                                                                                      >  1:4056541/1‑150 (MQ=255)

                                                  gcTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGt                                                                                                    >  1:3175250/1‑150 (MQ=255)

                                                   cTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTa                                                                                                   <  1:4959077/150‑1 (MQ=255)

                                                   cTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTa                                                                                                   <  1:4471600/150‑1 (MQ=255)

                                                   cTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTa                                                                                                   <  1:1185449/150‑1 (MQ=255)

                                                     tGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACa                                                                                                 <  2:4240984/150‑1 (MQ=255)

                                                      gAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAg                                                                                                <  1:3879477/150‑1 (MQ=255)

                                                           gCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCt                                                                                           >  1:4834928/1‑150 (MQ=255)

                                                             tctTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGc                                                                                         >  2:2300367/1‑150 (MQ=255)

                                                               ttCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGccc                                                                                       <  2:4838000/150‑1 (MQ=255)

                                                               ttCGCTGGGGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGccc                                                                                       <  2:2905666/150‑1 (MQ=255)

                                                                 cGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCt                                                                                     <  2:1473735/150‑1 (MQ=255)

                                                                    tggtggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGGCAGGGCAATGAGGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCGTCGTACAGCCCGTACAGCTCCTGCCCCTtgt                                                                                  >  2:5159410/1‑150 (MQ=255)

                                                                    tggtggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTtgt                                                                                  >  2:2590119/1‑150 (MQ=255)

                                                                    tggtggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTtgt                                                                                  >  1:4770553/1‑150 (MQ=255)

                                                                    tggtggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTtgt                                                                                  <  1:4766993/150‑1 (MQ=255)

                                                                      gtggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTgg                                                                                >  1:2474074/1‑150 (MQ=255)

                                                                             tGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCg                                                                         <  2:2913775/150‑1 (MQ=255)

                                                                              gCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCgg                                                                        <  2:1688216/150‑1 (MQ=255)

                                                                               ccAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCggg                                                                       <  1:533171/150‑1 (MQ=255)

                                                                               ccAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCggg                                                                       <  1:3196928/150‑1 (MQ=255)

                                                                               ccAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCggg                                                                       <  1:4175611/150‑1 (MQ=255)

                                                                               ccAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCggg                                                                       <  1:3043311/150‑1 (MQ=255)

                                                                               ccAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCggg                                                                       <  1:2680571/150‑1 (MQ=255)

                                                                                cAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCgggg                                                                      >  1:1812678/1‑150 (MQ=255)

                                                                                 aaGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGt                                                                     <  1:1316729/150‑1 (MQ=255)

                                                                                 aaGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGt                                                                     >  2:832652/1‑150 (MQ=255)

                                                                                 aaGGGTGGCCACGGCATAGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGt                                                                     >  2:522614/1‑150 (MQ=255)

                                                                                    ggTGGCCACGGCATTTCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCGTGCGCGCCTCGTACATCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTc                                                                  >  2:1108777/1‑150 (MQ=255)

                                                                                     gTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTcc                                                                 >  1:3541115/1‑150 (MQ=255)

                                                                                        gCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCggg                                                              <  2:358373/150‑1 (MQ=255)

                                                                                         ccACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCgggg                                                             >  2:405588/1‑150 (MQ=255)

                                                                                         ccACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCgggg                                                             >  2:2180449/1‑150 (MQ=255)

                                                                                          cACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGa                                                            <  1:2300367/150‑1 (MQ=255)

                                                                                           aCGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCg                                                                >  1:3548362/1‑145 (MQ=255)

                                                                                           aCGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAg                                                           >  2:4690279/1‑150 (MQ=255)

                                                                                           aCGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAg                                                           >  1:3463080/1‑150 (MQ=255)

                                                                                           aCGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAg                                                           >  1:2383333/1‑150 (MQ=255)

                                                                                              gaaTTGCGCAGGCATTGCTGTGCCAGGGCAATGACGTCAATGTATCCATAGACTAAGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCATCGTACAGCCCGTAAAGCTCCTGCCCCTTGTGGAAAACCGGCGTTTCCGGTGAGTTc                                                        <  2:3039403/148‑1 (MQ=255)

                                                                                               cATTGCGCAGGCCTTGCTGGGTCAGGGCAATGACGTCCATGTAGCCCTCGACGCCGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCa                                                       <  2:1740853/150‑1 (MQ=255)

                                                                                                     gcAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACt                                                 >  2:149572/1‑150 (MQ=255)

                                                                                                      cAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACtt                                                >  1:176466/1‑150 (MQ=255)

                                                                                                      cAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACtt                                                <  2:1000064/150‑1 (MQ=255)

                                                                                                         gCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTggg                                             >  2:1303247/1‑150 (MQ=255)

                                                                                                           cTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCt                                           >  1:1472239/1‑150 (MQ=255)

                                                                                                               cTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTgtc                                       <  2:4334994/150‑1 (MQ=255)

                                                                                                                tGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTgtcg                                      <  2:4442346/150‑1 (MQ=255)

                                                                                                                 gggCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTgtcgt                                     >  2:1928217/1‑150 (MQ=255)

                                                                                                                        ggCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAg                              >  1:1558995/1‑150 (MQ=255)

                                                                                                                        ggCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAg                              <  2:3515383/150‑1 (MQ=255)

                                                                                                                        ggCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAg                              >  2:1758316/1‑150 (MQ=255)

                                                                                                                              gACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCg                        >  2:2827288/1‑150 (MQ=255)

                                                                                                                              gACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCg                        <  2:4995867/150‑1 (MQ=255)

                                                                                                                                 gTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGcc                     >  2:3005643/1‑150 (MQ=255)

                                                                                                                                 gTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGcc                     >  1:2941820/1‑150 (MQ=255)

                                                                                                                                   ccATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCAc                   >  1:4576078/1‑150 (MQ=255)

                                                                                                                                      tGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGa                <  2:3543515/150‑1 (MQ=255)

                                                                                                                                       gTAGCCCTCGACGACGGTGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTtgtg                                                                                 >  1:4952314/1‑84 (MQ=255)

                                                                                                                                       gTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGaa               >  1:3931059/1‑150 (MQ=255)

                                                                                                                                        tAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAg              <  2:2513613/150‑1 (MQ=255)

                                                                                                                                          gCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGc            <  2:4939582/150‑1 (MQ=255)

                                                                                                                                            ccTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGaa               >  2:4424877/1‑145 (MQ=255)

                                                                                                                                            ccTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCga          >  1:1986472/1‑150 (MQ=255)

                                                                                                                                            ccTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCga          >  2:592491/1‑150 (MQ=255)

                                                                                                                                            ccTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCga          <  1:5164112/150‑1 (MQ=255)

                                                                                                                                            ccTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCga          >  2:2895117/1‑150 (MQ=255)

                                                                                                                                            ccTCGACGA‑GATCATCGTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCgg                                                                        >  1:1597964/1‑88 (MQ=38)

                                                                                                                                            ccTCGACGA‑GATCATCGTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCgg                                                                        <  2:1597964/88‑1 (MQ=38)

                                                                                                                                             cTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCgat         <  1:2827288/150‑1 (MQ=255)

                                                                                                                                              tCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCgatg        >  1:646035/1‑150 (MQ=255)

                                                                                                                                              tCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCgatg        >  2:3658484/1‑150 (MQ=255)

                                                                                                                                              tCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCgatg        >  2:2396597/1‑150 (MQ=255)

                                                                                                                                                gacgacgaTGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCgatgat      <  1:4746530/150‑1 (MQ=255)

CTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCGATGATGCGC  >  NC_002947/474046‑474343

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