BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,848,918	C→T	9.1%	A106V (GCC→GTC)	cysM →	cysteine synthase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,848,918	0	C	T	9.1%	293.8 / 15.9	110	A106V (GCC→GTC)	cysM	cysteine synthase B
Reads supporting (aligned to +/- strand): ref base C (41/58); new base T (4/6); total (46/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

Only 100 of 119 total aligned reads displayed.

TTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAACCCGATCGCCCACTACA  >  NC_002947/1848773‑1849067

ttGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATg                                                                                                                                                   >  2:1934208/1‑150 (MQ=255)

     gggCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGc                                                                                                                                              <  2:1589229/150‑1 (MQ=255)

     gggCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGc                                                                                                                                              >  2:398167/1‑150 (MQ=255)

         cAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGTCATTTCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGTCATGGCGGCCttg                                                                                                                                          <  1:25228/150‑3 (MQ=255)

         cAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTAt                                                                                                                                          <  1:1164956/150‑1 (MQ=255)

          aGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATg                                                                                                                                         <  2:1270071/150‑1 (MQ=255)

           gATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAAAACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATgg                                                                                                                                        >  2:654247/1‑150 (MQ=255)

            aTCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGAc                                                                                                                                                 >  1:861170/1‑140 (MQ=255)

            aTCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGgc                                                                                                                                       <  1:1909846/150‑1 (MQ=255)

             tCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGTCATGGCGGCCtt                                                                                                                                           >  2:2423408/1‑144 (MQ=255)

             tCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGgcg                                                                                                                                      >  1:980222/1‑150 (MQ=255)

                 gCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGa                                                                                                                                  <  2:2253039/150‑1 (MQ=255)

                   cAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGc                                                                                                                                <  2:907386/150‑1 (MQ=255)

                        gATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATc                                                                                                                           >  2:2389597/1‑150 (MQ=255)

                         aTACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATcc                                                                                                                          >  2:1835977/1‑150 (MQ=255)

                           aCCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGCCGGCCTATGGCGCCGAGCTGATCCTg                                                                                                                        <  1:322823/150‑1 (MQ=255)

                           aCCCTGATCGAAGCCACCTCCGGCAACACGGGCAGCGCGCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACACCTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTg                                                                                                                        <  2:568999/150‑1 (MQ=255)

                             ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTg                                                                                                                        <  2:159340/148‑1 (MQ=255)

                             ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTg                                                                                                                        >  1:159340/1‑148 (MQ=255)

                             ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGt                                                                                                                      <  1:398167/150‑1 (MQ=255)

                             ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGt                                                                                                                      >  1:2431954/1‑150 (MQ=255)

                             ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGt                                                                                                                      >  1:1514513/1‑150 (MQ=255)

                             ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGt                                                                                                                      <  2:1620089/150‑1 (MQ=255)

                              cTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTg                                                                                                                     >  2:1013566/1‑150 (MQ=255)

                                   cGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCaa                                                                                                                <  1:1835977/150‑1 (MQ=255)

                                        ccACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAaggagg                                                                                                           <  1:1564693/150‑1 (MQ=255)

                                          aCCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAaggaggag                                                                                                         >  2:2227996/1‑150 (MQ=255)

                                           ccTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAaggaggagg                                                                                                        <  2:710502/150‑1 (MQ=255)

                                             tCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAgggg                                                                                                      <  1:2389597/150‑1 (MQ=255)

                                              ccGGCAACGCGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGa                                                                                                     <  1:60376/150‑1 (MQ=255)

                                              ccGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGa                                                                                                     <  1:847256/150‑1 (MQ=255)

                                                 gCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATgg                                                                                                  >  1:1100811/1‑150 (MQ=255)

                                                   aaCACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGaa                                                                                                <  2:1100811/150‑1 (MQ=255)

                                                   aaCACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGaa                                                                                                >  2:2135749/1‑150 (MQ=255)

                                                            aTCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCAATAACGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGt                                                                                       >  2:1883511/1‑150 (MQ=255)

                                                                cccTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGAcc                                                                                   >  1:2342268/1‑150 (MQ=255)

                                                                cccTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGAcc                                                                                   >  1:1480647/1‑150 (MQ=255)

                                                                  cTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCCAGGAGGAGGGGATGGAAGGCGCCCGTGACCTg                                                                                 <  1:1013566/150‑1 (MQ=255)

                                                                          ggcggcggaGATCAATGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCCAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGaa                                                                         <  1:1333551/150‑1 (MQ=255)

                                                                             ggcggcgATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCt                                                                      >  1:1843885/1‑150 (MQ=255)

                                                                              gcggcgATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGc                                                                                                                                <  1:1386006/91‑1 (MQ=255)

                                                                              gcggcgATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGc                                                                                                                                >  2:1386006/1‑91 (MQ=255)

                                                                              gcggcgATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTg                                                                     >  1:108371/1‑150 (MQ=255)

                                                                               cggcgATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGc                                                                    <  1:9817/150‑1 (MQ=255)

                                                                                ggcgATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCa                                                                   <  1:654247/150‑1 (MQ=255)

                                                                                 gcgATCAAGGGCTACAAGATGATCCTGATCATGCCCGACCACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAg                                                                  <  1:1203415/150‑1 (MQ=255)

                                                                                  cgATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAgg                                                                 <  1:1111608/150‑1 (MQ=255)

                                                                                     tCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGgcgc                                                                                           <  2:1504249/121‑1 (MQ=255)

                                                                                        aGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGgcgc                                                                                                                                     <  2:962322/76‑1 (MQ=255)

                                                                                        aGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGgcgc                                                                                                                                     >  1:962322/1‑76 (MQ=255)

                                                                                        aGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAg                                                           <  2:2232371/150‑1 (MQ=255)

                                                                                             tACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCTAAGGCCGt                                                      <  1:2436526/150‑1 (MQ=255)

                                                                                               cAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTgg                                                    >  1:1558157/1‑150 (MQ=255)

                                                                                               cAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTgg                                                    >  2:925120/1‑150 (MQ=255)

                                                                                                aaGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGc                                                   >  1:1370981/1‑150 (MQ=255)

                                                                                                 aGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGcc                                                  <  2:965184/150‑1 (MQ=255)

                                                                                                  gatgatCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAGGGCCGCCATGACGGCGAATCGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGACGGCCGTGGCCt                                                 >  2:82421/1‑150 (MQ=255)

                                                                                                  gatgatCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCt                                                 <  2:164108/150‑1 (MQ=255)

                                                                                                  gatgatCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCt                                                 <  2:289192/150‑1 (MQ=255)

                                                                                                      ctCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTg                                             <  2:1046756/149‑1 (MQ=255)

                                                                                                       tCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGc                                                   <  1:1396101/143‑1 (MQ=255)

                                                                                                       tCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGc                                                   <  1:1399081/143‑1 (MQ=255)

                                                                                                       tCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGc                                                   >  2:1396101/1‑143 (MQ=255)

                                                                                                       tCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGc                                                   >  2:1399081/1‑143 (MQ=255)

                                                                                                       tCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGc                                            <  2:1919705/150‑1 (MQ=255)

                                                                                                       tCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGAAGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGc                                            <  1:2366562/150‑1 (MQ=255)

                                                                                                         cTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTc                                          <  1:2228892/150‑1 (MQ=255)

                                                                                                         cTGATCATGCCCGACAACTCCACCGCCGAGAGCAAAGCCGCCATNACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGcag                                          >  1:1046756/1‑148 (MQ=255)

                                                                                                          tGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCg                                         <  1:925120/150‑1 (MQ=255)

                                                                                                          tGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCg                                         >  2:2324050/1‑150 (MQ=255)

                                                                                                           gATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGa                                        <  1:2227996/150‑1 (MQ=255)

                                                                                                              cATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCa                                     <  2:2418823/150‑1 (MQ=255)

                                                                                                                tGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGt                                   <  1:1948998/150‑1 (MQ=255)

                                                                                                                  cccGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTc                                 <  2:695092/150‑1 (MQ=255)

                                                                                                                   ccGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAaggaggagg                                                                                                        >  1:2452696/1‑78 (MQ=255)

                                                                                                                   ccGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCg                                >  2:922623/1‑150 (MQ=255)

                                                                                                                    cGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGc                               <  1:2135749/150‑1 (MQ=255)

                                                                                                                    cGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGc                               <  2:2431954/150‑1 (MQ=255)

                                                                                                                     gACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGcc                              >  1:713912/1‑150 (MQ=255)

                                                                                                                           tCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACg                          <  1:560864/148‑1 (MQ=255)

                                                                                                                           tCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACg                          >  2:560864/1‑148 (MQ=255)

                                                                                                                           tCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGc                        >  1:565863/1‑150 (MQ=255)

                                                                                                                           cccACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGc                        <  1:1277934/149‑1 (MQ=255)

                                                                                                                           cccACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGc                        <  1:1934208/149‑1 (MQ=255)

                                                                                                                             cACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATTGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGa                      <  2:791024/150‑1 (MQ=255)

                                                                                                                                 gccgAGCGCAAGGCCGAGATGAGGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAAc                  >  2:2041102/1‑150 (MQ=255)

                                                                                                                                  ccgAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAAcc                 <  2:1377474/150‑1 (MQ=255)

                                                                                                                                      gccataGGCCGTCATGGCGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTg                                                     >  2:773340/6‑110 (MQ=255)

                                                                                                                                       cgcAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTg                                                     <  1:1726898/109‑1 (MQ=255)

                                                                                                                                       cgcAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAACCCGATc            <  1:922623/150‑1 (MQ=255)

                                                                                                                                       ccataGGCCGTCATGGCGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGGGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGAc                                       <  2:2414537/119‑1 (MQ=255)

                                                                                                                                       ccataGGCCGTCATGGCGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGGGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGAc                                       >  1:2414537/5‑123 (MQ=255)

                                                                                                                                           aGGCCGTCATGGCGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGa                                                                                     >  2:2274932/1‑73 (MQ=39)

                                                                                                                                             gccgccATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAg                                                                                               >  1:1456679/1‑61 (MQ=255)

                                                                                                                                             gccgccATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAg                                                                                               <  2:1456679/61‑1 (MQ=255)

                                                                                                                                             gccgccATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAACCCGATCGCCCAc      <  1:2324050/150‑1 (MQ=255)

                                                                                                                                             gccgccATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAACCCGATCGCCCAc      >  1:1607477/1‑150 (MQ=255)

                                                                                                                                               cgccATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAg                                                           <  1:1267882/95‑1 (MQ=255)

                                                                                                                                                 ccATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCTTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAACCCGATCGCCCACTaca  >  1:1806452/1‑150 (MQ=255)

                                                                                                                                                 ccATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGAGCTGGCTGAGAAGCTCATTGAAGAATTCCGTTGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAACCCGATCGCCCACTaca  >  1:168212/1‑150 (MQ=255)

TTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAACCCGATCGCCCACTACA  >  NC_002947/1848773‑1849067

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