BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,802,077	T→C	5.1%	I110V (ATC→GTC)	PP_4227 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,802,077	0	T	C	5.1%	634.1 / 13.7	195	I110V (ATC→GTC)	PP_4227	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (101/84); new base C (7/3); total (108/87)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 203 total aligned reads displayed.

TTAATACTACTTGCACAGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATACGCGCCGACGCCGCCAGGC  >  NC_002947/4801928‑4802226

ttAATACTACTTGCACAGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGAt                                                                                                                                                       >  1:5351898/1‑150 (MQ=255)

ttAATACTACTTGCACAGATCTCCACCCTCGAACCCTAGAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGAt                                                                                                                                                       <  1:513981/150‑1 (MQ=255)

 tAATACTACTTGCACAGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCCTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGAtc                                                                                                                                                      <  1:5184976/150‑1 (MQ=255)

   aTACTACTTGCACAGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGAtctc                                                                                                                                                    <  2:120385/150‑1 (MQ=255)

     actactTGCACAGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGc                                                                                                                                                  >  1:4283232/1‑150 (MQ=255)

      ctactTGCACAGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGcc                                                                                                                                                 <  2:3851351/150‑1 (MQ=255)

       tactTGCACAGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGATCCTCGCCa                                                                                                                                                <  2:4492630/150‑1 (MQ=255)

           tGCACAGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCa                                                                                                                                            <  2:1911145/150‑1 (MQ=255)

                aGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGAGGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAg                                                                                                                                       >  1:132679/1‑150 (MQ=255)

                   tctcCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGctg                                                                                                                                    >  2:5265336/1‑150 (MQ=255)

                     tcCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGctgac                                                                                                                                  >  2:3326969/1‑150 (MQ=255)

                       cACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGctgacct                                                                                                                                >  1:4257914/1‑150 (MQ=255)

                        aCCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGctgacctg                                                                                                                               >  1:227140/1‑150 (MQ=255)

                        aCCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCATGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGctgacctg                                                                                                                               >  1:1174909/1‑150 (MQ=255)

                           cTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACa                                                                                                                            >  2:1985157/1‑150 (MQ=255)

                           cTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACa                                                                                                                            <  2:2209581/150‑1 (MQ=255)

                            tCGACCCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCCTGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGTCAg                                                                                                                           <  1:5227079/150‑1 (MQ=255)

                            tCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAg                                                                                                                           >  1:2541845/1‑150 (MQ=255)

                            tCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAg                                                                                                                           >  1:570742/1‑150 (MQ=255)

                            tCGAACCCTATAGCTGTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAg                                                                                                                           <  2:3495655/150‑1 (MQ=255)

                             cGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGGGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAgg                                                                                                                          <  1:929429/150‑1 (MQ=255)

                                  ccTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGgccatg                                                                                                                     <  1:5179250/150‑1 (MQ=255)

                                  ccTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGgccatg                                                                                                                     <  1:5088735/150‑1 (MQ=255)

                                    tataGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGgccatgcc                                                                                                                   >  2:3452208/1‑150 (MQ=255)

                                     ataGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGgccatgcca                                                                                                                  >  2:5381088/1‑150 (MQ=255)

                                       aGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATc                                                                                                                <  1:3326969/150‑1 (MQ=255)

                                       aGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATc                                                                                                                <  1:83000/150‑1 (MQ=255)

                                            ttCTGTCTTCTAACCACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTa                                                                                                           <  2:4573886/150‑1 (MQ=255)

                                               tcacatatAAACACAGTATTGGCCGGGAAAAAGCGGCTGAAGAGATGCACGGCATAGTCCCGGCGTCCATGACGCCAACTTGCTTTTGGTTGCTCGCGATGATCTCGCCACGCAGCGAGCTGACCAGACAGGCCATGCCATCCGGAACGc                                                                                                        >  2:1258650/8‑150 (MQ=21)

                                               tGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGc                                                                                                        >  2:3842991/1‑150 (MQ=255)

                                                   ttCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGct                                                                                                                                     >  2:2484676/1‑117 (MQ=255)

                                                   ttCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGaa                                                                                                    >  2:143572/1‑150 (MQ=255)

                                                   ttCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGaa                                                                                                    >  2:2895729/1‑150 (MQ=255)

                                                   ttCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGaa                                                                                                    >  2:4164675/1‑150 (MQ=255)

                                                   ttCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGaa                                                                                                    >  2:4330353/1‑150 (MQ=255)

                                                   ttCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGaa                                                                                                    >  1:109065/1‑150 (MQ=255)

                                                   ttCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGAGGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGaa                                                                                                    >  1:761837/1‑150 (MQ=255)

                                                    tCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAAc                                                                                                   >  2:118752/1‑150 (MQ=255)

                                                        aaCACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGt                                                                                               >  1:208713/1‑150 (MQ=255)

                                                        aaCACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGt                                                                                               >  2:16429/1‑150 (MQ=255)

                                                        aaCACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGt                                                                                               >  2:16125/1‑150 (MQ=255)

                                                           acaGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAAc                                                                                                   <  2:3964859/143‑1 (MQ=255)

                                                           acaGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAAc                                                                                                   >  1:3964859/1‑143 (MQ=255)

                                                             aGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTcggcg                                                                                          <  2:1459083/150‑1 (MQ=255)

                                                                aTTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGaa                                                                                       <  1:1860752/150‑1 (MQ=255)

                                                                    gCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAAt                                                                                   <  2:2115422/150‑1 (MQ=255)

                                                                      cGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATgc                                                                                 >  1:4210391/1‑150 (MQ=255)

                                                                       gggAAAAAGCGGCGGATGAGATGTACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATgcg                                                                                >  2:3399016/1‑150 (MQ=255)

                                                                       gggAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATgcg                                                                                <  1:2620800/150‑1 (MQ=255)

                                                                       gggAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATgcg                                                                                <  2:2764187/150‑1 (MQ=255)

                                                                        ggAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATgcgc                                                                               >  1:2454989/1‑150 (MQ=255)

                                                                                cggcggATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGc                                                                       >  1:455026/1‑150 (MQ=255)

                                                                                 ggcggATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGcc                                                                      >  1:1032788/1‑150 (MQ=255)

                                                                                     gATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTg                                                                  >  1:4077774/1‑150 (MQ=255)

                                                                                       tGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGc                                                                >  2:4667430/1‑150 (MQ=255)

                                                                                        gagaTGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCg                                                               <  1:152746/150‑1 (MQ=255)

                                                                                           aTGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTc                                                            >  1:3874011/1‑150 (MQ=255)

                                                                                              cACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGc                                                         >  1:5005508/1‑150 (MQ=255)

                                                                                               aCGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCg                                                        >  1:3234763/1‑150 (MQ=255)

                                                                                                cGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCgg                                                       >  2:1033653/1‑150 (MQ=255)

                                                                                                  gCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCt                                                     >  2:3215074/1‑150 (MQ=255)

                                                                                                    ttCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTcggcg                                                                                          >  2:2847697/1‑111 (MQ=255)

                                                                                                    ttCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTcggcg                                                                                          <  1:2847697/111‑1 (MQ=255)

                                                                                                    ttCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGc                                                   >  1:1841958/1‑150 (MQ=255)

                                                                                                        tCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCa                                               >  1:2050624/1‑150 (MQ=255)

                                                                                                          ccGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGc                                                         >  1:4568128/1‑138 (MQ=40)

                                                                                                          ccGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAAc                                             <  1:2906867/150‑1 (MQ=255)

                                                                                                          ccGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAAc                                             >  2:4296583/1‑150 (MQ=255)

                                                                                                            ggCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTg                                           <  2:4890035/150‑1 (MQ=255)

                                                                                                              cGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGt                                         <  2:4543152/150‑1 (MQ=255)

                                                                                                              cGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGt                                         <  1:2895729/150‑1 (MQ=255)

                                                                                                              cGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGt                                         <  2:3011202/150‑1 (MQ=255)

                                                                                                               gTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTc                                        <  2:4220916/150‑1 (MQ=255)

                                                                                                                  cATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGAt                                     >  1:4014333/1‑150 (MQ=255)

                                                                                                                  atggACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGAt                                     >  1:3586471/4‑150 (MQ=255)

                                                                                                                   aTGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTc                                                            >  1:5217063/1‑126 (MQ=38)

                                                                                                                   aTGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTc                                                            <  2:5217063/126‑1 (MQ=38)

                                                                                                                   aTGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAAc                                             >  1:1627589/1‑141 (MQ=38)

                                                                                                                   aTGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAAc                                             <  2:1627589/141‑1 (MQ=38)

                                                                                                                   aTGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATc                                    >  1:3172647/1‑150 (MQ=255)

                                                                                                                              gTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGAc                         >  1:3403223/1‑150 (MQ=39)

                                                                                                                               tGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGAcc                        <  1:3259626/150‑1 (MQ=39)

                                                                                                                               tGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGAcc                        >  2:1941908/1‑150 (MQ=39)

                                                                                                                               tGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGACAACTGGTCGATCAGCCCGGCGAcc                        <  2:570742/150‑1 (MQ=255)

                                                                                                                                gCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCa                       >  2:1368636/1‑150 (MQ=39)

                                                                                                                                 cTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCAt                      >  1:4728955/1‑150 (MQ=39)

                                                                                                                                  tctttgatGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATa                     >  2:173829/8‑150 (MQ=255)

                                                                                                                                   tgtgGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATAc                    <  2:4588004/150‑1 (MQ=39)

                                                                                                                                     tgGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATAcgc                  <  2:761867/150‑1 (MQ=39)

                                                                                                                                     tgGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATAcgc                  <  1:1918403/150‑1 (MQ=39)

                                                                                                                                       gTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATACgcgc                >  2:2392554/1‑150 (MQ=39)

                                                                                                                                       ctgcgtGGCGAGATCCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATACgcgc                >  1:5229943/6‑150 (MQ=255)

                                                                                                                                        ttGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATACGCGcc               <  1:274739/150‑1 (MQ=39)

                                                                                                                                        tgcgtGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATgc                                                                                 <  2:3789142/80‑1 (MQ=255)

                                                                                                                                        tgcgtGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATgc                                                                                 >  1:3789142/5‑84 (MQ=255)

                                                                                                                                        tgcgtGGCGAGATCCTCGCCA‑GCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCACCTGGTCGATCAGCCCGGCGACCATACGCGCCg              >  2:5147562/5‑150 (MQ=11)

                                                                                                                                           gtGGCGAGATCCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATACGCGCCGAc            >  1:4774805/2‑150 (MQ=255)

                                                                                                                                            tGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCg                           >  2:382250/1‑134 (MQ=40)

                                                                                                                                              gCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATACGCGCCGAcgcc         <  1:3215074/150‑1 (MQ=39)

                                                                                                                                                  ggATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATACGCGCCGACGCCGCCa     >  1:2151635/1‑150 (MQ=255)

TTAATACTACTTGCACAGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATACGCGCCGACGCCGCCAGGC  >  NC_002947/4801928‑4802226

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