BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,141,550:1	+A	5.0%	coding (424/972 nt)	PP_2757 →	sugar‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,141,550	1	.	A	5.0%	406.7 / 13.9	100	coding (424/972 nt)	PP_2757	sugar‑binding protein
Reads supporting (aligned to +/- strand): ref base . (50/45); new base A (2/3); total (52/48)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 104 total aligned reads displayed.

CGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGCCGATCATCCGCAATTGCAGATCGTCCAGCC  >  NC_002947/3141404‑3141698

cGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGCCAACG‑CGCGGCCGGCcgcgacgt                                                                                                                                                   <  2:2076902/150‑5 (MQ=255)

    aaTTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGc                                                                                                                                               <  1:1389033/150‑1 (MQ=255)

         aGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAAc                                                                                                                                          <  1:832148/150‑1 (MQ=255)

           aCCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCAGCTGTTCGGCCAACG‑CGCGGCCGGCCGCGACGTTGTCCGag                                                                                                                                           <  2:448764/150‑2 (MQ=25)

            cccTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACCACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGc                                                                                                                                       <  1:1598935/150‑1 (MQ=255)

             ccTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCt                                                                                                                                      >  2:1380887/1‑150 (MQ=255)

                    aCCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTgcgc                                                                                                                               <  1:75103/150‑1 (MQ=255)

                    aCCCGGAGCGTGATCGACCCCTGGCTCAAGCGCAGCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTgcgc                                                                                                                               <  2:360153/150‑1 (MQ=255)

                      ccTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCAGCG‑CGTTGACCTAACAGCTGGTGCGCGa                                                                                                                             >  1:1001059/1‑150 (MQ=255)

                       cTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGac                                                                                                                            >  2:2326238/1‑150 (MQ=255)

                              tGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGTCGAACAGCTGGTGCGCGACACcggcg                                                                                                                     <  2:1793180/150‑1 (MQ=255)

                                 tCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTgc                                                                                                                                 >  2:1925456/1‑135 (MQ=255)

                                 tCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTgc                                                                                                                                 <  1:1925456/135‑1 (MQ=255)

                                  cGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACG‑CGCGGCCGGCcgcgac                                                                                                                                                     >  1:840885/1‑112 (MQ=255)

                                   gACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAg                                                                                                                <  1:2326238/150‑1 (MQ=255)

                                    aCCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACG‑CGCGGCCGGCCGCGACGTTgtc                                                                                                                                               >  1:2443458/1‑116 (MQ=255)

                                    aCCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACG‑CGCGGCCGGCCGCGACGTTgtc                                                                                                                                               <  2:2443458/118‑3 (MQ=255)

                                             tCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATcc                                                                                                      >  1:676383/1‑150 (MQ=255)

                                             tCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACAACGGCGGCAAGGGCCGCATcc                                                                                                      >  2:425111/1‑150 (MQ=255)

                                             tCAAGCGCATCCAGGCAGCAGGTATCCCACTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCTCG‑CGTTGGCCGAACAGCTAGTGCGCGACACCGGCGGCAAGGGCCGCATcc                                                                                                      >  2:1251248/1‑150 (MQ=255)

                                               aaGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTg                                                                                                    <  1:1097681/150‑1 (MQ=255)

                                               aaGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTg                                                                                                    <  2:676383/150‑1 (MQ=255)

                                                   ccATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGgtgt                                                                                                <  2:2128471/149‑1 (MQ=255)

                                                      tCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCa                                                                                             >  2:1784706/1‑150 (MQ=255)

                                                       ccAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCaa                                                                                            <  1:1784706/150‑1 (MQ=255)

                                                        cAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAAc                                                                                           <  1:1380887/150‑1 (MQ=255)

                                                         aGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACg                                                                                          >  1:1336265/1‑150 (MQ=255)

                                                                   tATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACgg                                                                                >  2:2359422/1‑150 (MQ=255)

                                                                    aTCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTg                                                                                                                                     <  2:1395584/96‑1 (MQ=255)

                                                                    aTCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTg                                                                                                                                     >  1:1395584/1‑96 (MQ=255)

                                                                       ccGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGgt                                                                                                  <  1:1845433/128‑1 (MQ=255)

                                                                       ccGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGgt                                                                                                  >  2:1845433/1‑128 (MQ=255)

                                                                       ccGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTg                                                                            >  1:766377/1‑150 (MQ=255)

                                                                        cGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGTCAAGAGCCGCATCCTGGTGGGCAACGGCGTCTGCGGtgtgg                                                                           >  2:2360706/1‑149 (MQ=255)

                                                                             ttCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGCCACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTc                                                                      <  1:95380/150‑1 (MQ=255)

                                                                              tCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCTGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCt                                                                     <  1:2249453/150‑1 (MQ=255)

                                                                                  cATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGc                                                                 >  1:1870129/1‑150 (MQ=255)

                                                                                   aTTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCg                                                                >  1:727022/1‑150 (MQ=255)

                                                                                   aTTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGACGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCCGTGCg                                                                >  2:1288526/1‑150 (MQ=255)

                                                                                    ttGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGa                                                               >  2:1472020/1‑150 (MQ=255)

                                                                                      gATGTGCCATCGCAGTACAGCCTGAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTACTGCGGGGTGCCGGTCTGTGCGATc                                                             >  2:1193442/1‑150 (MQ=255)

                                                                                        tgtgCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTgg                                                                                                   >  1:575225/1‑110 (MQ=255)

                                                                                        tgtgCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTgg                                                                                                   <  2:575225/110‑1 (MQ=255)

                                                                                        tgtgCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCg                                                           >  1:498561/1‑150 (MQ=255)

                                                                                           gCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACgggg                                                                              >  2:806052/1‑128 (MQ=255)

                                                                                           gCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACgggg                                                                              <  1:806052/128‑1 (MQ=255)

                                                                                             cATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGc                                                                                        <  2:1956424/116‑1 (MQ=255)

                                                                                             cATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGc                                                                                        >  1:1956424/1‑116 (MQ=255)

                                                                                              aTCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGGCACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTacga                                                     >  2:1658955/1‑150 (MQ=255)

                                                                                              aTCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTg                                                                                                    >  2:2394107/1‑103 (MQ=255)

                                                                                              aTCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTg                                                                                                    <  1:2394107/103‑1 (MQ=255)

                                                                                              aTCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTa                                                                                   <  2:1968543/120‑1 (MQ=255)

                                                                                              aTCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTa                                                                                   >  1:1968543/1‑120 (MQ=255)

                                                                                              aTCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTacga                                                     <  1:1422856/150‑1 (MQ=255)

                                                                                              aTCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTacga                                                     >  1:1286560/1‑150 (MQ=255)

                                                                                              aTCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTacga                                                     >  1:2393766/1‑150 (MQ=255)

                                                                                              aTCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTacga                                                     <  2:1336265/150‑1 (MQ=255)

                                                                                               tCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTacgac                                                    >  2:1985321/1‑150 (MQ=255)

                                                                                                 gCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTg                                                                                                                                     <  1:490112/67‑1 (MQ=255)

                                                                                                 gCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTg                                                                                                                                     >  2:490112/1‑67 (MQ=255)

                                                                                                 gCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTAcgacga                                                  <  1:2359422/150‑1 (MQ=255)

                                                                                                 gCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTAcgacga                                                  <  2:727022/150‑1 (MQ=255)

                                                                                                 gCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTAcgacga                                                  <  1:1678421/150‑1 (MQ=255)

                                                                                                   aGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGt                                                <  2:430338/150‑1 (MQ=255)

                                                                                                     tACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTg                                              >  1:239679/1‑150 (MQ=255)

                                                                                                       cAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTgaa                                            <  2:239679/150‑1 (MQ=255)

                                                                                                        aGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTgaag                                           >  2:754705/1‑150 (MQ=255)

                                                                                                        aGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTgaag                                           >  1:1612275/1‑150 (MQ=255)

                                                                                                          ccTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTgaagaa                                         >  1:2105572/1‑150 (MQ=255)

                                                                                                          ccTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTgaagaa                                         <  1:1319567/150‑1 (MQ=255)

                                                                                                          ccTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCAACCCGGTCAGTCGAGttcgagag                                         >  2:1319567/1‑144 (MQ=255)

                                                                                                          ccTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCGGTGCGATCCGCTGCGACGAGTTgaagaa                                         >  1:1303447/1‑150 (MQ=255)

                                                                                                           cTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGaagaag                                        <  2:498561/150‑1 (MQ=255)

                                                                                                            tCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGt                                                >  2:2023756/1‑141 (MQ=255)

                                                                                                            tCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGt                                                <  1:2023756/141‑1 (MQ=255)

                                                                                                            tCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAgg                                       >  2:36144/1‑150 (MQ=255)

                                                                                                                caacaCCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACt                                   <  2:1286560/150‑1 (MQ=255)

                                                                                                                 aacaCCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTg                                  >  2:2441362/1‑150 (MQ=255)

                                                                                                                   caccacC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCGACGACGAGTTGAAGAAGGTACTGGc                                <  2:1077040/150‑1 (MQ=255)

                                                                                                                    accacc‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAg                                                 <  2:2032775/132‑1 (MQ=255)

                                                                                                                    accacc‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAg                                                 >  1:2032775/1‑132 (MQ=255)

                                                                                                                     ccacc‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGCCg                              <  2:2393766/150‑1 (MQ=255)

                                                                                                                      cacc‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCtgtg                                                                  <  2:608241/113‑1 (MQ=255)

                                                                                                                      cacc‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCtgtg                                                                  >  1:608241/1‑113 (MQ=255)

                                                                                                                      cacc‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGCCgc                             >  2:4646/1‑149 (MQ=255)

                                                                                                                       acc‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACg                                                                                          >  1:190495/1‑88 (MQ=255)

                                                                                                                       acc‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACg                                                                                          <  2:190495/88‑1 (MQ=255)

                                                                                                                       acc‑‑‑GCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGCCGat                            <  2:1001059/150‑1 (MQ=255)

                                                                                                                         c‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACAGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGCCGatca                          >  2:695125/1‑150 (MQ=255)

                                                                                                                               ggACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGTGAAGTTTTAGAAGGTACTGGCCGATCATcc                       >  1:1181437/1‑150 (MQ=255)

                                                                                                                                gACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGCCGATCATCCg                      >  1:345454/1‑150 (MQ=255)

                                                                                                                                gACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGCCGATCATACg                      >  1:479258/1‑150 (MQ=255)

                                                                                                                                 aCAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGCCGATCATCCGc                     >  2:1483021/1‑150 (MQ=255)

                                                                                                                                 aCAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATACGCTACGACGAGTTGAAGAAGGTACTGGCCGATCATCCGc                     >  2:1483582/1‑150 (MQ=255)

                                                                                                                                    aacgcgcgGCCGGCCGCGACGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGCCGa                             >  1:669468/5‑140 (MQ=255)

                                                                                                                                    aacgcgcgGCCGGCCGCGACGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGCCGa                             <  2:669468/136‑1 (MQ=255)

                                                                                                                                        cgcgGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGCCGATCATCCGCAATTGCa              <  1:1666006/150‑1 (MQ=255)

                                                                                                                                          cggccggcCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGCCGATCATCCGCAATTGCAGa            <  2:1870129/150‑1 (MQ=255)

                                                                                                                                            gccggccgCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGCCGATCATCCGCAATTGCAGATc          <  1:1658955/150‑1 (MQ=255)

                                                                                                                                                    cg‑cgTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGCCGATCATCCGCAATTGCAGATCGTCCAGcc  <  2:1612275/150‑1 (MQ=255)

CGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACC‑‑‑TCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGCCGATCATCCGCAATTGCAGATCGTCCAGCC  >  NC_002947/3141404‑3141698

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