BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,848,839	C→A	5.2%	L80M (CTG→ATG)	cysM →	cysteine synthase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,848,839	0	C	A	5.2%	495.4 / 8.6	151	L80M (CTG→ATG)	cysM	cysteine synthase B
Reads supporting (aligned to +/- strand): ref base C (77/66); new base A (4/4); total (81/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 158 total aligned reads displayed.

CCTCCTGCTCAAGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTG  >  NC_002947/1848691‑1848988

cctcctGCTCAAGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCt                                                                                                                                                         <  1:1026035/150‑1 (MQ=255)

  tcctGCTCAAGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCTAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTgg                                                                                                                                                       >  2:2120849/1‑150 (MQ=255)

  tcctGCTCAAGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTgg                                                                                                                                                       <  1:1378005/150‑1 (MQ=255)

     tGCTCAAGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGa                                                                                                                                                    >  1:2751525/1‑150 (MQ=255)

     tGCTCAAGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCCATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGa                                                                                                                                                    >  2:2252642/1‑150 (MQ=255)

          aaGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATggcg                                                                                                                                               <  1:612040/150‑1 (MQ=255)

           aGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCTTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCAGGGCGATgcccg                                                                                                                                              <  1:1869783/150‑5 (MQ=255)

           aGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATggcgg                                                                                                                                              <  2:1215351/150‑1 (MQ=255)

           aGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATggcgg                                                                                                                                              <  2:1375135/150‑1 (MQ=255)

            gCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATggcggc                                                                                                                                             >  1:1157475/1‑150 (MQ=255)

                 aaGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGa                                                                                                                                        <  2:2751525/150‑1 (MQ=255)

                      aaCAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑Ag                                                                                                                                  >  2:1852552/1‑150 (MQ=255)

                        cAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATg                                                                                                                                                  >  2:703781/1‑133 (MQ=255)

                         aaTCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGc                                                                                                                               >  1:2694227/1‑150 (MQ=255)

                         aaTCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGCCATCGCCAGGGCGATgcccgt                                                                                                                                             >  2:1510579/1‑132 (MQ=255)

                         aaTCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGCCATCGCCAGGGCGATgcccgt                                                                                                                                             <  1:1510579/137‑6 (MQ=255)

                          aTCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACGCGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑Agg                                                                                                                                 <  1:1827566/147‑1 (MQ=255)

                          aTCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACGCGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑Agg                                                                                                                                 >  2:1827566/1‑147 (MQ=255)

                           tCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTa                                                                                                                             >  1:2261591/1‑150 (MQ=255)

                            ccTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTAc                                                                                                                            <  2:3020340/150‑1 (MQ=255)

                             cTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGCCATCGCCAGGGCGATgcccgt                                                                                                                                             >  2:2176632/1‑128 (MQ=255)

                             cTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGCCATCGCCAGGGCGATgcccgt                                                                                                                                             <  1:2176632/133‑6 (MQ=255)

                              tGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGGCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGAAGAGGCGACCACAGGGcagaa                                                                                                                           >  2:172983/1‑146 (MQ=255)

                              tGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACaa                                                                                                                          <  2:14506/150‑1 (MQ=255)

                                     tCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACATGATGATc                                                                                                                   <  1:2719553/150‑1 (MQ=255)

                                     tCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATc                                                                                                                   <  2:2694227/150‑1 (MQ=255)

                                     tCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATc                                                                                                                   >  1:742324/1‑150 (MQ=255)

                                      cGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAga                                                                                                                        >  1:381833/1‑144 (MQ=255)

                                      cGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAga                                                                                                                        <  2:381833/144‑1 (MQ=255)

                                        gTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTg                                                                                                                >  1:794820/1‑150 (MQ=255)

                                         tGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATg                                                                                                                                                  >  1:1523304/1‑116 (MQ=255)

                                         tGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATg                                                                                                                                                  <  2:1523304/116‑1 (MQ=255)

                                          gAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCTTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAgat                                                                                                                       <  2:3001314/141‑1 (MQ=255)

                                          gAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCTTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAgat                                                                                                                       >  1:3001314/1‑141 (MQ=255)

                                          gAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATg                                                                                                                                                  >  1:2440101/1‑115 (MQ=255)

                                          gAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATg                                                                                                                                                  <  2:2440101/115‑1 (MQ=255)

                                             aGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCAt                                                                                                           >  2:1091415/1‑150 (MQ=255)

                                             aGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGACAGGCGATACCCTGATAGCAACTACTGCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCAt                                                                                                           >  1:1145162/1‑150 (MQ=255)

                                              gACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATAGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATg                                                                                                          <  1:2858625/150‑1 (MQ=255)

                                               aCCGCCCGGCCCTGTCGATGATTGCTCGCGCCGAATTGCGGGGACAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGc                                                                                                         <  1:662075/150‑1 (MQ=255)

                                                   cccGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGa                                                                                                     >  2:2569432/1‑150 (MQ=255)

                                                    ccGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAgatg                                                                                                                      <  2:203711/132‑1 (MQ=255)

                                                    ccGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAgatg                                                                                                                      >  1:203711/1‑132 (MQ=255)

                                                       gCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAAc                                                                                                 >  2:53638/1‑150 (MQ=255)

                                                           tGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCa                                                                                             <  2:742324/150‑1 (MQ=255)

                                                               gatgatCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCAccgc                                                                                         >  1:2802829/1‑150 (MQ=255)

                                                                atgatCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCAccgcc                                                                                        >  1:78756/1‑150 (MQ=255)

                                                                  gatCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGa                                                                                      <  2:266252/150‑1 (MQ=255)

                                                                   atCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCACTGGCGATGGCGGCGGCGATCA‑Agg                                                                                                                                 <  2:1174052/106‑1 (MQ=255)

                                                                        tcgcGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑Ca‑a                                                                              >  1:479214/1‑150 (MQ=255)

                                                                        tcgcGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑Ca‑a                                                                              >  1:1018877/1‑150 (MQ=255)

                                                                        tcgcGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑Ca‑a                                                                              <  2:1157475/150‑1 (MQ=255)

                                                                           cgcCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATggcgg                                                                                                                                              <  1:2092936/86‑1 (MQ=255)

                                                                            gcCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGgcc                                                                          >  2:362548/1‑150 (MQ=255)

                                                                                      ggggCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑Ag                                                                                                                                  >  1:1315682/1‑86 (MQ=255)

                                                                                       gggCCAGATCAAGCAAGGCGCTCCACTGATGGAAGCCACCACCTGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGc                                                               <  2:1568648/150‑1 (MQ=255)

                                                                                          ccAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGCCATCGCCAGGGCGATgccc                                                                                                                                               <  2:1941703/70‑4 (MQ=37)

                                                                                            aGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATg                                                          <  1:2629845/150‑1 (MQ=255)

                                                                                             gATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATgg                                                         >  2:846580/1‑150 (MQ=255)

                                                                                             gATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCACGAACTCTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATgg                                                         >  1:2337928/1‑150 (MQ=255)

                                                                                              aTCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGgc                                                        >  1:702998/1‑150 (MQ=255)

                                                                                               tCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑Ag                                                                                                                                  <  1:1734374/77‑1 (MQ=255)

                                                                                                  aGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCCATAGGCCGTCATGGCGGCCttgcgctc                                                      <  2:78756/150‑8 (MQ=255)

                                                                                                     cAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGc                                                 >  1:175152/1‑150 (MQ=255)

                                                                                                     aaGGAGAAACCCTAATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGc                                                 <  2:173976/149‑1 (MQ=255)

                                                                                                        gCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACt                                                                                                >  1:158506/1‑102 (MQ=255)

                                                                                                        gCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACt                                                                                                <  2:158506/102‑1 (MQ=255)

                                                                                                        gCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGa                                              >  2:803910/1‑150 (MQ=255)

                                                                                                              cccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTgg                                        <  2:1018877/150‑1 (MQ=255)

                                                                                                               ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGgc                                                        <  2:1313375/133‑1 (MQ=255)

                                                                                                               ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGgc                                                        >  1:1313375/1‑133 (MQ=255)

                                                                                                               ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGc                                                 <  2:916988/140‑1 (MQ=255)

                                                                                                               ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGc                                                 >  1:916988/1‑140 (MQ=255)

                                                                                                               ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGt                                       <  1:2451983/150‑1 (MQ=255)

                                                                                                               ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGt                                       <  1:1190947/150‑1 (MQ=255)

                                                                                                               ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGt                                       <  1:2120849/150‑1 (MQ=255)

                                                                                                               ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGt                                       >  1:2150070/1‑150 (MQ=255)

                                                                                                               ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATAGAGCCGCGCTGAtggtggg                                       >  2:2477417/1‑144 (MQ=255)

                                                                                                               ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACCGCCTATGGCGCCGAGCTGATCCTGGt                                       >  2:1073726/1‑150 (MQ=255)

                                                                                                               ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGCCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGCCCTATGGCGCCGACCTGCTCCTGGt                                       >  1:2161885/1‑150 (MQ=255)

                                                                                                                 tGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGa                                     >  1:944350/1‑150 (MQ=255)

                                                                                                                   aTCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAAc                                                                                                 >  1:7208/1‑90 (MQ=255)

                                                                                                                    tCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCt                                          <  2:582752/142‑1 (MQ=255)

                                                                                                                      gAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAag                                >  1:1932661/1‑150 (MQ=255)

                                                                                                                       aaGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAagg                               <  1:803910/150‑1 (MQ=255)

                                                                                                                          ccACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAaggagg                            >  2:1461092/1‑150 (MQ=255)

                                                                                                                            aCCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAaggaggag                          >  1:354781/1‑150 (MQ=255)

                                                                                                                              cTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAaggag                             >  2:2899129/1‑145 (MQ=255)

                                                                                                                              cTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAggg                        >  1:2334012/1‑150 (MQ=255)

                                                                                                                               tCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAgggg                       <  1:1852552/150‑1 (MQ=255)

                                                                                                                                  ggCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATg                    <  2:2807984/150‑1 (MQ=255)

                                                                                                                                     aaCACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGaa                 <  1:1461092/150‑1 (MQ=255)

                                                                                                                                       cacGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAgg               >  1:3071184/1‑150 (MQ=255)

                                                                                                                                        acGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAAGAGGAGGGGATGGAAGgc              >  1:2153016/1‑150 (MQ=255)

                                                                                                                                          gggCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGgcgc            >  2:1769258/1‑150 (MQ=255)

                                                                                                                                               tCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAg                <  2:2355735/141‑1 (MQ=255)

                                                                                                                                               tCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAg                >  1:2355735/1‑141 (MQ=255)

                                                                                                                                                  cccTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAaggag                             <  2:2504264/125‑1 (MQ=255)

                                                                                                                                                  cccTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCAACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAaggag                             >  1:2504264/1‑125 (MQ=255)

                                                                                                                                                    cTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGACCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTg  <  1:2911202/150‑1 (MQ=255)

CCTCCTGCTCAAGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCA‑AGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCG‑CA‑AGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTG  >  NC_002947/1848691‑1848988

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