BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	733,650	A→T	9.6%	intergenic (+21/‑26)	PP_t09 → / → PP_t10	tRNA‑Lys/tRNA‑Pro

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	733,650	0	A	T	9.6%	326.0 / 20.3	122	intergenic (+21/‑26)	PP_t09/PP_t10	tRNA‑Lys/tRNA‑Pro
Reads supporting (aligned to +/- strand): ref base A (55/53); new base T (7/5); total (62/60)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

Only 100 of 129 total aligned reads displayed.

GCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACC  >  NC_002947/733501‑733797

gCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCaaa                                                                                                                                                      <  1:3907900/150‑1 (MQ=255)

gCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCaaa                                                                                                                                                      <  1:3138947/150‑1 (MQ=255)

 cTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACACACCATGTAGAAAAGGGCAAATCAAag                                                                                                                                                     >  1:4916371/1‑150 (MQ=255)

  tttGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAAGGCAATCTCAAtga                                                                                                                                                    >  2:4460106/1‑150 (MQ=11)

    tGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAgaga                                                                                                                                                  <  1:2287727/150‑1 (MQ=255)

         gCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGccc                                                                                                                                             <  2:4344699/150‑1 (MQ=255)

         gCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGccc                                                                                                                                             >  2:1236349/1‑150 (MQ=255)

             gATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑tttt                                                                                                                                        <  2:2555077/150‑1 (MQ=255)

                gTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑aattttt                                                                                                                                     <  2:1663808/150‑1 (MQ=255)

                 tCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑tttttttt                                                                                                                                    >  1:3610009/1‑150 (MQ=255)

                    ggAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCCAAGAGATGCCC‑TTTTTTTTGaa                                                                                                                                 <  1:4109149/150‑1 (MQ=255)

                     gAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGaaa                                                                                                                                <  2:3299488/150‑1 (MQ=255)

                          cGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTc                                                                                                                           >  2:4627673/1‑150 (MQ=255)

                           gAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGTCC‑TTTTTTTTGAAACACTCg                                                                                                                          <  2:1982619/150‑1 (MQ=255)

                           gAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCg                                                                                                                          <  1:1909394/150‑1 (MQ=255)

                            caTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGACGTAGGGTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCgg                                                                                                                         <  2:663877/149‑1 (MQ=255)

                             aTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCCTGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCggg                                                                                                                        >  1:51849/1‑150 (MQ=255)

                               ccccTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCCAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGgt                                                                                                                      <  1:4120377/150‑1 (MQ=255)

                                cccTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGgtg                                                                                                                     <  1:1088137/150‑1 (MQ=255)

                                    cACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAgc                                                                                                                 <  2:2237066/150‑1 (MQ=255)

                                    cACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAgc                                                                                                                 <  2:4232964/150‑1 (MQ=255)

                                     cccccGACCTTTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCCAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAgcg                                                                                                                <  1:5033380/149‑1 (MQ=255)

                                       cccGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCa                                                                                                              >  2:3594651/1‑150 (MQ=255)

                                          gACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCCAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTc                                                                                                           <  1:347051/150‑1 (MQ=255)

                                          gACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTAAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTc                                                                                                           >  1:5268755/1‑150 (MQ=255)

                                             cATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCTTTTGGTCGTAGGTTCGAATCCTACACGACCCACCCTGTAGAAAAGGGCATCTCCAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCgg                                                                                                        <  1:70591/150‑1 (MQ=255)

                                             cATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTGTTTTGTAACGCTCGGGGTGTAGCGCAGTCCgg                                                                                                        >  2:2943090/1‑150 (MQ=255)

                                              aTTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGt                                                                                                       >  1:1771573/1‑150 (MQ=255)

                                              aTTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGt                                                                                                       >  1:1385086/1‑150 (MQ=255)

                                                   ccGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAgcgc                                                                                                  <  2:4980459/150‑1 (MQ=255)

                                                           ttAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCg                                                                                                         >  1:3866122/1‑135 (MQ=255)

                                                           ttAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCg                                                                                                         <  2:3866122/135‑1 (MQ=255)

                                                            tAGCTCAGTGGGTAGAGGAGGTGGCTTTGAACCAATTGGTCGTAGGTTCGACTCCTCCACGACCCACCCTGTAGAAAAGGGCAGCTCCAAGCGATGCCG‑TTCTTTTTGACACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCttt                                                                                         <  1:5160033/150‑1 (MQ=25)

                                                            tAGCTCAGCGGGCAGAGCAGTTGGCTTTTAACAAATTGGTCGTAGGTTCGAATCCTACACCACCCACCATGTATAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCttt                                                                                         <  2:5345218/150‑1 (MQ=25)

                                                              gCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                       >  1:995135/1‑150 (MQ=255)

                                                              gCTCAGGTGGTAGAGCAGTTGGCTTTTCACCACGTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGACAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTTAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                       <  2:3846703/150‑1 (MQ=25)

                                                                 cAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGAAg                                                                                    >  1:4458500/1‑150 (MQ=255)

                                                                  aGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAATGGCCAGAACACAGAGATGCCC‑TTTTTATTGACACACTCGCGGAGACGCCCAGACCGGTCGCGCGCCAGCTTGGGGAGc                                                                                   >  2:3007971/1‑150 (MQ=17)

                                                                  aGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGc                                                                                   <  1:4627673/150‑1 (MQ=255)

                                                                   gTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCttt                                                                                         <  2:4146245/143‑1 (MQ=255)

                                                                   gTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCCCAGTCCGGTAGCGCGCCTGCTTTGGGAGCa                                                                                  <  2:939710/150‑1 (MQ=255)

                                                                    ttGGGAGCGCAGTTGGCTTTCAACCAATTGGTCGTAGGTGCGAATCCTACACGACCCACCCTGTAGAAAAGGGCATCTCAACGAGCTTCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAg                                                                                 <  1:2777417/150‑1 (MQ=38)

                                                                      ggTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGCCAACTCAAAGAGAAGACC‑TTTTTATAGAAAAACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGa                                                                               >  2:2875860/1‑150 (MQ=25)

                                                                            gCAGTTTGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCCAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAg                                                                         <  1:2347099/150‑1 (MQ=255)

                                                                            gCAGTTGGCTATTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCCTGTAGAAAAGGGCATCTCAAAGATATGCAC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAg                                                                         <  2:3411629/150‑1 (MQ=255)

                                                                                  ggCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTc                                                                   <  2:5337408/150‑1 (MQ=255)

                                                                                     ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGCCC‑t                                                                                                                                           >  1:4183527/1‑75 (MQ=17)

                                                                                     ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGCCC‑t                                                                                                                                           <  2:4183527/75‑1 (MQ=17)

                                                                                     ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTTAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGaa                                                                >  2:1154751/1‑150 (MQ=255)

                                                                                     ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGaa                                                                >  2:3793453/1‑150 (MQ=255)

                                                                                     ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGaa                                                                >  1:4686271/1‑150 (MQ=255)

                                                                                     ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGaa                                                                >  1:4534871/1‑150 (MQ=255)

                                                                                      tttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAAt                                                               <  1:712759/150‑1 (MQ=255)

                                                                                        tacccAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGTGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATcc                                                             <  1:1176480/147‑1 (MQ=255)

                                                                                           ccAATTGGTCGTAGGTTCGAATCCTACACGACCCAGCATGTAGCAAAGGGCAGCTCCAAGCGTTGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCt                                                          <  2:564423/150‑1 (MQ=25)

                                                                                               ttGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCAcc                                                      <  2:3610009/150‑1 (MQ=255)

                                                                                                tggtggTAGGTTCGTATGCTACACGTCCCTCCTTTTATAATAGGGCATCTGAAAGTGTTGCAC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCAccc                                                     <  2:2770837/150‑1 (MQ=18)

                                                                                                   tCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGa                                                  >  1:1284091/1‑150 (MQ=255)

                                                                                                    cGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGCCAACACAAAGAGATGCCC‑TTTTTTTTTAAACACTCTGTGTGTAGCCCAGTCCGGTAGCGCGCCTGCTTGTGGAGCAGGATTTCATGAGGGAGAATCCCCTCACCCCGAc                                                 >  1:1752730/1‑150 (MQ=25)

                                                                                                     gTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGAc                                                 >  1:903873/1‑150 (MQ=255)

                                                                                                     gTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGAcc                                                >  1:691070/1‑150 (MQ=255)

                                                                                                     gTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGAcc                                                >  2:463369/1‑150 (MQ=255)

                                                                                                       aGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCAt                                              >  1:4067318/1‑150 (MQ=255)

                                                                                                       aGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCAt                                              >  2:637216/1‑150 (MQ=255)

                                                                                                         gTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGgt                                                                                                                      >  1:2912470/1‑76 (MQ=255)

                                                                                                          ttCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCAtttt                                           >  1:1916486/1‑150 (MQ=255)

                                                                                                             gAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTgg                                        >  1:1771850/1‑150 (MQ=255)

                                                                                                              aaTCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTggg                                       >  1:4039064/1‑150 (MQ=255)

                                                                                                                tCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTTAAACACTCGGGGTGAAGCGCAGTCCGGTAGCGCGCCTGCGTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATGTTTGGGTc                                     >  2:2674416/1‑150 (MQ=255)

                                                                                                                tCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTc                                     >  2:3718807/1‑150 (MQ=255)

                                                                                                                   tACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCa                                                                                  <  2:3069745/102‑1 (MQ=255)

                                                                                                                    acacGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTa                                 <  1:463369/150‑1 (MQ=255)

                                                                                                                        gACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTc                             >  1:5222505/1‑150 (MQ=255)

                                                                                                                         aCCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCa                            <  2:1771850/150‑1 (MQ=255)

                                                                                                                           ccaccaTGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGt                          >  1:1023322/1‑150 (MQ=255)

                                                                                                                           ccaccaTGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGt                          >  1:1263632/1‑150 (MQ=255)

                                                                                                                             cccaTGTAGAAAAGGGCTTCTCACTGAGATGCCC‑TTATTTTTGCAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTg                        <  1:4893144/149‑1 (MQ=18)

                                                                                                                              ccaTGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTgg                       >  1:3365758/1‑150 (MQ=255)

                                                                                                                               caTGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGt                      >  2:3601820/1‑150 (MQ=255)

                                                                                                                               caTGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGt                      <  2:1385086/150‑1 (MQ=255)

                                                                                                                                aTGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTa                     <  1:3594651/150‑1 (MQ=255)

                                                                                                                                  gTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTaga                   <  2:4686271/150‑1 (MQ=255)

                                                                                                                                  gTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTaga                   <  2:4534871/150‑1 (MQ=255)

                                                                                                                                  aaaaaaaaGGGCATCTCTTTGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGGTGGTaga                   >  2:2751736/5‑150 (MQ=18)

                                                                                                                                   tAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTagag                  <  2:1023322/150‑1 (MQ=255)

                                                                                                                                    aaaaaaGGGCATCTCTTTGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGAGCAGTCCGGTAGCGCGCATGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGc                 >  2:1308251/3‑150 (MQ=18)

                                                                                                                                      aaaaGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTCGAGCAg               >  2:4370740/1‑150 (MQ=255)

                                                                                                                                      aaaaGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAg               >  2:1411834/1‑150 (MQ=255)

                                                                                                                                        aaGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCAttttt                                          >  2:3414344/1‑121 (MQ=255)

                                                                                                                                          gggCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGGTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTgg           >  1:3219313/1‑150 (MQ=255)

                                                                                                                                            gCATCTCTTTGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCt         >  2:54954/1‑150 (MQ=255)

                                                                                                                                            gCATCTCTTTGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGGTCGAATCCCCTCACCCCGACCATGTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCt         >  2:54980/1‑150 (MQ=17)

                                                                                                                                            gCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCt         >  1:1018391/1‑150 (MQ=255)

                                                                                                                                             cATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCtt        <  1:3573579/150‑1 (MQ=255)

                                                                                                                                             cATCACACAGAGTTTTCC‑TATTTTTTGAAACTCTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCtt        <  2:4990153/150‑1 (MQ=14)

                                                                                                                                              aTCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACTATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCttt       >  1:3854103/1‑150 (MQ=255)

                                                                                                                                              aTCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCttt       <  2:3365758/150‑1 (MQ=255)

                                                                                                                                                 tcAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTaa    >  2:1813874/1‑150 (MQ=255)

                                                                                                                                                   aaaGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGGTGGTAGAGCAGTTGGCTTTTAAcc  >  1:1864015/1‑150 (MQ=255)

                                                                                                                                                     agagATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCtt        >  1:1509890/1‑142 (MQ=255)

GCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACC  >  NC_002947/733501‑733797

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