BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,688,230	A→C	6.7%	K102Q (AAG→CAG)	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,230	0	A	C	6.7%	304.8 / 11.9	105	K102Q (AAG→CAG)	PP_3253	carboxylate‑amine ligase
Reads supporting (aligned to +/- strand): ref base A (52/46); new base C (2/5); total (54/51)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 108 total aligned reads displayed.

CGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTT  >  NC_002947/3688082‑3688378

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

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

    cTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGcc‑c                                                                                                                                                 >  2:2034998/1‑150 (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)

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

                     cACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACta                                                                                                                                >  1:1156385/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                                                                                                            >  2:1126495/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:1952886/150‑1 (MQ=255)

                                                                cAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGc                                                                                     <  1:1996253/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                                                                                   <  1:607371/141‑1 (MQ=255)

                                                                           cgTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCt                                                                                   >  2:607371/1‑141 (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                                                                                                                                                 >  1:1984592/1‑73 (MQ=255)

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

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

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

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

                                                                                          gAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCACTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCAc                                                             <  1:611346/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                                                >  1:1234609/1‑150 (MQ=255)

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

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

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

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

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

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

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

                                                                                                             aCGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCtt                                          <  1:1733877/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                               <  1:1934906/150‑1 (MQ=255)

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

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

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

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

                                                                                                                            aCCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAAc                           <  1:475389/150‑1 (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)

                                                                                                                                            gctgctGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCt           <  1:872623/150‑1 (MQ=255)

                                                                                                                                              tgctgcAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGAcc                                    <  2:544064/123‑1 (MQ=255)

                                                                                                                                              tgctgcAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGAcc                                    >  1:544064/1‑123 (MQ=255)

                                                                                                                                                  gcAGAAGCC‑CGCGGCAAAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTg     <  1:1983824/150‑1 (MQ=255)

                                                                                                                                                    aGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCGTGCAACTGATCAACCGCCTTCTGCCCTg     <  1:1200588/148‑1 (MQ=255)

                                                                                                                                                    aGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCGTGCAACTGATCAACCGCCTTCTGCCCTg     >  2:1200588/1‑148 (MQ=255)

                                                                                                                                                     gAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGtt  >  2:607032/1‑150 (MQ=255)

                                                                                                                                                     gAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGtt  >  1:428175/1‑150 (MQ=255)

                                                                                                                                                      aaGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCAc                                                             >  2:425085/1‑90 (MQ=255)

                                                                                                                                                      aaGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCAc                                                             <  1:425085/90‑1 (MQ=255)

CGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATT‑G‑GCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTT  >  NC_002947/3688082‑3688378

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: `‑`