BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,688,218	Δ1 bp	6.8%	coding (292/1107 nt)	PP_3253 →	carboxylate‑amine ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,688,218	0	T	.	6.8%	307.3 / 16.6	103	coding (292/1107 nt)	PP_3253	carboxylate‑amine ligase
Reads supporting (aligned to +/- strand): ref base T (51/44); new base . (2/5); total (53/50)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

Only 100 of 105 total aligned reads displayed.

TCCGCAGCCAGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCC  >  NC_002947/3688069‑3688365

tCCGCAGCCAGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGt                                                                                                                                                      >  1:815476/1‑150 (MQ=255)

tCCGCAGCCAGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGt                                                                                                                                                      <  2:815476/150‑1 (MQ=255)

   gCAGCCAGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGgc                                                                                                                                                   >  1:400894/1‑150 (MQ=255)

        cAGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGctgctg                                                                                                                                              >  1:1006251/1‑150 (MQ=255)

          gATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCtgctgc                                                                                                                                             <  1:287980/149‑1 (MQ=255)

          gATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCtgctgc                                                                                                                                             >  2:287980/1‑149 (MQ=255)

          gATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTACGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCAGATGCAGCGGCAGGCCAACCCTGGGCCTCGTGGCtgctgct                                                                                                                                            <  2:366151/150‑2 (MQ=255)

             ctaGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGaa                                                                                                                                         <  1:413457/148‑1 (MQ=255)

                gCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGcc                                                                                                                                      <  1:831349/150‑1 (MQ=255)

                         gCCGGTGTTCACCAACCTGTATGAAGCAAGTGAGTTTCTCCAGCGCAACCGTCAGCGGATGCCCGTGGCGTTCGTGACCTGGATAGCAAGGG‑‑CAAGGGCGGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGcaca                                                                                                                            >  1:2280901/1‑150 (MQ=255)

                          ccGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGTGCCGCG‑GGCTTCTGCAGCAGCCACGCGGcgcag                                                                                                                           <  1:47020/150‑2 (MQ=38)

                           cGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGTGCCGCG‑GGCTTCTGCAGCAGCCACGCGGcgcagg                                                                                                                          <  2:423209/150‑3 (MQ=38)

                             gtgtTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCGGAACAGGGGGTGGTGCTGTATGGGCCGGCAAGCCACGTCGGCGCCGCGTTGCTGCTGCAGAAGCC‑CGCGGCACAAggg                                                                                                                        <  2:2017651/150‑1 (MQ=255)

                             gtgtTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCACCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAggg                                                                                                                        <  2:1752566/150‑1 (MQ=255)

                                ttCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCAc                                                                                                                     >  2:1923652/1‑150 (MQ=255)

                                ttCACCAACCTGGATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCAc                                                                                                                     >  2:1399895/1‑150 (MQ=255)

                                 tCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCTTGTGCCGCG‑GGCTTCTGCAGCAGCCACGCGGCGCAggggtggc                                                                                                                    <  1:1989127/150‑5 (MQ=37)

                                  cACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACta                                                                                                                   >  1:1156385/1‑150 (MQ=255)

                                  cACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACta                                                                                                                   >  2:1449502/1‑150 (MQ=255)

                                   aCCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACtat                                                                                                                  >  2:1996253/1‑150 (MQ=255)

                                    ccAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACtata                                                                                                                 >  2:1934906/1‑150 (MQ=255)

                                      aaCCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGTGCACTATAag                                                                                                               <  1:510576/150‑1 (MQ=255)

                                      aaCCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAag                                                                                                               <  2:2154344/150‑1 (MQ=255)

                                        ccTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagca                                                                                                             >  2:1952886/1‑150 (MQ=255)

                                           gTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCt                                                                                                          >  1:1984977/1‑150 (MQ=255)

                                           gTATGAAGCACGTGACTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCt                                                                                                          <  2:1622853/150‑1 (MQ=255)

                                              tGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGtt                                                                                                       >  1:508089/1‑150 (MQ=255)

                                              tGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCATGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGtt                                                                                                       >  1:581037/1‑150 (MQ=255)

                                                  gCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTcgac                                                                                                   <  1:753953/150‑1 (MQ=255)

                                                    aCGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTcgacga                                                                                                 >  2:2010084/1‑150 (MQ=255)

                                                     cGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCgacgac                                                                                                >  2:1010794/1‑150 (MQ=255)

                                                      gTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACt                                                                                               >  1:2159235/1‑150 (MQ=255)

                                                          gTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCg                                                                                           >  1:77886/1‑150 (MQ=255)

                                                               tcCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCAtg                                                                                      >  2:938882/1‑150 (MQ=255)

                                                                 cAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑cg                                                                                                                                   <  1:1981819/103‑1 (MQ=255)

                                                                 cAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑cg                                                                                                                                   >  2:1981819/1‑103 (MQ=255)

                                                                  aGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGAGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTgg                                                                                   >  1:1705491/1‑150 (MQ=255)

                                                                   gcgcAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGcc                                                                                         <  2:1015541/143‑1 (MQ=255)

                                                                   gcgcAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGcc                                                                                         >  1:1015541/1‑143 (MQ=255)

                                                                      ccaCCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGACGCGTGGCTGCTGCAGAAGCC‑CCCGGCACAAGGGCACTATACGCTGCTGTTCGACGACTACCCCCATGTGACAAg                                                                               <  2:689744/148‑1 (MQ=255)

                                                                      cAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAg                                                                               <  2:1006251/150‑1 (MQ=255)

                                                                        aCCGTCAGCGGCTGCGCGTGGCATT‑G‑GACGAAGAGGGGATGGGGCCATACGCTGCCGAAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGc                                                                             <  2:42490/150‑1 (MQ=255)

                                                                         ccGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACAACTGCGCCGCTGGGTTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCg                                                                            <  1:2034998/150‑1 (MQ=255)

                                                                          cGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGTGCCGCG‑GGCTTCTGCAGCAGCCACGCgg                                                                                                                                <  2:1150516/97‑1 (MQ=25)

                                                                          cGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGTGCCGCG‑GGCTTCTGCAGCAGCCACGCgg                                                                                                                                >  1:1150516/1‑97 (MQ=25)

                                                                           gTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTa                                                                          <  1:1126495/150‑1 (MQ=255)

                                                                             cAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGc                                                                        <  1:1996253/150‑1 (MQ=255)

                                                                             cAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGc                                                                        <  1:1952886/150‑1 (MQ=255)

                                                                              aGCGGCTGCGCGTGGCATT‑G‑GCCGAAAAGGGGAGGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGcc                                                                       <  2:1705491/150‑1 (MQ=255)

                                                                                  gCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCt                                                                   <  1:1010794/150‑1 (MQ=255)

                                                                                      cgcgTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTTGCAAGGCGTAGCCTGCTCAAt                                                               <  1:1449502/150‑1 (MQ=255)

                                                                                        cgTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCt                                                                      >  2:607371/1‑141 (MQ=255)

                                                                                        cgTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCt                                                                      <  1:607371/141‑1 (MQ=255)

                                                                                         gTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGc                                                            <  2:2159235/150‑1 (MQ=255)

                                                                                             cATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGc                                                        <  2:1698541/150‑1 (MQ=255)

                                                                                              aTT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGcc‑c                                                                                                                                    <  2:1984592/73‑1 (MQ=255)

                                                                                              aTT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGcc‑c                                                                                                                                    >  1:1984592/1‑73 (MQ=255)

                                                                                                t‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCAtg                                                     <  1:854556/150‑1 (MQ=255)

                                                                                                  g‑gCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCAtgt                                                    >  1:447311/1‑150 (MQ=255)

                                                                                                  g‑gCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGAAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCAtgt                                                    >  2:1484353/1‑150 (MQ=255)

                                                                                                     ccGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTgg                                                                                                                                                    <  1:1690318/53‑1 (MQ=255)

                                                                                                     ccGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTgg                                                                                                                                                    >  2:1690318/1‑53 (MQ=255)

                                                                                                       gAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCG‑GGCTTCTGCAGCAGCCACGCgg                                                                                                                                >  1:109630/1‑70 (MQ=21)

                                                                                                       gAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCG‑GGCTTCTGCAGCAGCCACGCgg                                                                                                                                <  2:109630/70‑1 (MQ=21)

                                                                                                       gAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCACTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCAc                                                <  1:611346/150‑1 (MQ=255)

                                                                                                           agGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTg                                            <  1:1484353/150‑1 (MQ=255)

                                                                                                            ggggATGGGGCCATACGCTGCGGCAAGCCACACCTGCGCCGCGTGGCTGCTGAAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTGAATGGCCTGCATGTGCACGTTgg                                           >  1:974825/1‑150 (MQ=255)

                                                                                                             gggATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGc                                          <  2:400894/150‑1 (MQ=255)

                                                                                                             gggATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGc                                          <  2:1984977/150‑1 (MQ=255)

                                                                                                               gATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGt                                        >  2:1049615/1‑150 (MQ=255)

                                                                                                               gATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTCGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGt                                        >  2:893861/1‑150 (MQ=255)

                                                                                                                  ggggCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTgcc                                     >  1:1530646/1‑150 (MQ=255)

                                                                                                                    ggCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTCATCTTCTCAATGGCCTGCATGTGCACGTTGGCGTgccgc                                   >  1:122738/1‑150 (MQ=255)

                                                                                                                    ggCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTgccgc                                   >  2:1263710/1‑150 (MQ=255)

                                                                                                                    ggCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTgccgc                                   >  1:1234609/1‑150 (MQ=255)

                                                                                                                     gCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAATGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTgccgcc                                  >  1:142386/1‑150 (MQ=255)

                                                                                                                     gCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTgccgcc                                  >  1:1810435/1‑150 (MQ=255)

                                                                                                                       cATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTcgac                                                                                                   >  1:703697/1‑83 (MQ=255)

                                                                                                                       cATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTcgac                                                                                                   <  2:703697/83‑1 (MQ=255)

                                                                                                                        aTACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTcgac                                                                                                   <  2:2147009/82‑1 (MQ=255)

                                                                                                                        aTACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTcgac                                                                                                   >  1:2147009/1‑82 (MQ=255)

                                                                                                                          aCGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCtt                             <  1:1733877/150‑1 (MQ=255)

                                                                                                                          aCGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCtt                             <  1:1584381/150‑1 (MQ=255)

                                                                                                                          aCGCTGCGGCAAGCCACCCCTGCGCCGCATGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCtt                             <  2:581037/150‑1 (MQ=255)

                                                                                                                            gCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCTACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTtgt                           <  2:189505/150‑1 (MQ=255)

                                                                                                                            gCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTtgt                           <  1:877601/150‑1 (MQ=255)

                                                                                                                             cTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTtgtg                          >  1:2171137/1‑150 (MQ=255)

                                                                                                                              tGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGa                         >  1:2199779/1‑150 (MQ=255)

                                                                                                                                cGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTCCCGCCATGTGGCAAGGCGTCGAGGGCGCAGTAGTCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGAcc                       <  1:988956/150‑1 (MQ=255)

                                                                                                                                  gCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGc                     >  2:475389/1‑150 (MQ=255)

                                                                                                                                     aGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATg                  >  2:455662/1‑150 (MQ=255)

                                                                                                                                     aGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATg                  <  1:1934906/150‑1 (MQ=255)

                                                                                                                                      gCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGc                 <  1:2105379/150‑1 (MQ=255)

                                                                                                                                       ccACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCa                >  1:2101734/1‑150 (MQ=255)

                                                                                                                                          ccccTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACt             >  2:2065726/1‑150 (MQ=255)

                                                                                                                                            ccTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGa           <  1:1263710/150‑1 (MQ=255)

                                                                                                                                            ccTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGGAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGa           >  1:1313410/1‑150 (MQ=255)

                                                                                                                                             cTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGAt          >  2:568658/1‑150 (MQ=255)

                                                                                                                                                 gccgcGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAAc      >  1:907147/1‑150 (MQ=255)

                                                                                                                                                     cgTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGcc  <  2:86840/150‑1 (MQ=255)

TCCGCAGCCAGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCC  >  NC_002947/3688069‑3688365

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