BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,848,923	A→G	8.4%	T108A (ACG→GCG)	cysM →	cysteine synthase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,848,923	0	A	G	8.4%	332.5 / 14.0	119	T108A (ACG→GCG)	cysM	cysteine synthase B
Reads supporting (aligned to +/- strand): ref base A (48/60); new base G (4/6); total (52/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01

Only 100 of 120 total aligned reads displayed.

GGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAACCCGATCGCCCACTACAACAGC  >  NC_002947/1848778‑1849072

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

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

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

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

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

        tCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGTCATGGCGGCCtt                                                                                                                                                <  1:2423408/145‑2 (MQ=255)

        tCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGTCATGGCGGCCtt                                                                                                                                                >  2:2423408/1‑144 (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)

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

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

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

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

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

                                   ccACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCCCCGCCGAGCGCACTGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAagacgg                                                                                                                >  1:523957/1‑146 (MQ=255)

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

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

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

                                       cTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAggg                                                                                                            >  1:965184/1‑150 (MQ=255)

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

                                         ccGGCAACGCGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGa                                                                                                          <  1:60376/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:1480647/1‑150 (MQ=255)

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

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

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

                                                                      gcggcggcgATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAg                                                                             <  1:1719055/150‑1 (MQ=255)

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

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

                                                                         gcggcgATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGc                                                                                                                                     <  1:1386006/91‑1 (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)

                                                                            gcgATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAg                                                                       <  1:1203354/150‑1 (MQ=255)

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

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

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

                                                                                 cAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGa                                                                  >  1:1919705/1‑150 (MQ=255)

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

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

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

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

                                                                                          cAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTgg                                                         >  1:1558157/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)

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

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

                                                                                                                gACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGcc                                   >  1:713912/1‑150 (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                                                          <  1:773340/105‑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)

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

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

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

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

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

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

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

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

                                                                                                                                                gACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAACCCGATCGCCCACTACAACAg   >  1:1423123/1‑150 (MQ=255)

                                                                                                                                                gACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAACCCGATCGCCCACTACAACAg   >  1:230115/1‑150 (MQ=255)

                                                                                                                                                 aCGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGACGAAGGCCGTGGCCTAGTGCTCGACAAGTTCGCCAACAGCAACAACCCGATCGCCCACTACAACAGc  <  2:1843885/150‑1 (MQ=255)

GGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGATGGAAGGCGCCCGTGACCTGGCCGAGAAGCTGCAGGCCGAAGGCCGTGGCCTGGTGCTCGACCAGTTCGCCAACGGCGACAACCCGATCGCCCACTACAACAGC  >  NC_002947/1848778‑1849072

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