BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,750,775	C→G	6.9%	A31G (GCC→GGC)	PP_1551 →	replication protein O

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,750,775	0	C	G	6.9%	578.0 / 24.5	175	A31G (GCC→GGC)	PP_1551	replication protein O
Reads supporting (aligned to +/- strand): ref base C (97/66); new base G (8/4); total (105/70)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 197 total aligned reads displayed.

ACGAAAAAAAACCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCG  >  NC_002947/1750626‑1750922

aCGAAAAAAAACCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTcccggc                                                                                                                                                     >  2:2172598/1‑150 (MQ=255)

   aaaaaaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTcccggcccg                                                                                                                                                  <  2:265633/150‑1 (MQ=255)

   aaaaaaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTcccggcccg                                                                                                                                                  <  1:4273485/150‑1 (MQ=255)

    aaaaaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTcccggcccgg                                                                                                                                                 >  2:2309457/1‑150 (MQ=255)

    aaaaaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTcccggcccgg                                                                                                                                                 <  1:1395424/150‑1 (MQ=255)

    aaaaaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTcccggcccgg                                                                                                                                                 <  2:5375524/150‑1 (MQ=255)

    aaaaaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTcccggcccgg                                                                                                                                                 >  1:4694161/1‑150 (MQ=255)

    aaaaaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTcccggcccgg                                                                                                                                                 <  1:2369475/150‑1 (MQ=255)

       aaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTcccggccc                                                                                                                                                   <  1:5225849/145‑1 (MQ=255)

       aaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGgag                                                                                                                                              >  1:3092874/1‑150 (MQ=255)

          aCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGGCCGGGAGAtc                                                                                                                                           <  1:4460374/150‑2 (MQ=255)

           ccGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGa                                                                                                                                          >  2:2828251/1‑150 (MQ=255)

           ccGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGa                                                                                                                                          >  1:4852462/1‑150 (MQ=255)

            cGCCTGGCAGGGCGGCTCTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTGGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGaa                                                                                                                                         <  2:2366286/150‑1 (MQ=255)

                tGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGAtt                                                                                                                                     <  2:131205/150‑1 (MQ=255)

                tGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGAtt                                                                                                                                     <  1:3103232/150‑1 (MQ=255)

                tGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGAtt                                                                                                                                     >  2:4872506/1‑150 (MQ=255)

                       gcggcTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGt                                                                                                                              >  1:4046403/1‑150 (MQ=255)

                         ggcTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTAc                                                                                                                            <  1:2491226/150‑1 (MQ=255)

                            tttCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTg                                                                                                                         >  2:3194891/1‑150 (MQ=255)

                                  tACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAg                                                                                                                   >  2:4602135/1‑150 (MQ=255)

                                  tACAGCATTAGAACGGGGTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAg                                                                                                                   >  2:1218141/1‑150 (MQ=255)

                                   aCAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAgg                                                                                                                  >  1:5177012/1‑150 (MQ=255)

                                              aCGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCTTGGccacgtacagc                                                                                                       <  2:3515276/150‑8 (MQ=255)

                                                 ggTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTc                                                                                                    >  2:381373/1‑150 (MQ=255)

                                                     aaaGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTg                                                                                                >  1:5024917/1‑150 (MQ=255)

                                                       aGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCg                                                                                              >  2:140815/1‑150 (MQ=255)

                                                             aCAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGGGGCCAAGGCCACCATCAACTTCGGTGCGGgcgcc                                                                                        <  1:285519/150‑2 (MQ=255)

                                                               aaaCATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGga                                                                                                                                               >  2:5131036/1‑93 (MQ=255)

                                                               aaaCATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGt                                                                                                                          <  2:340817/114‑1 (MQ=255)

                                                               aaaCATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGt                                                                                                                          >  1:340817/1‑114 (MQ=255)

                                                                 aCATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAAc                                                                                    <  2:3262650/150‑1 (MQ=255)

                                                                 aCATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAAc                                                                                    <  2:4422098/150‑1 (MQ=255)

                                                                 aCATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAAc                                                                                    >  1:119743/1‑150 (MQ=255)

                                                                 aCATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGAGCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAAc                                                                                    >  2:3316227/1‑150 (MQ=255)

                                                                       tCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAAAGCATcc                                                                              >  2:31341/1‑150 (MQ=255)

                                                                        cccACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATccc                                                                             >  2:261999/1‑150 (MQ=255)

                                                                          cACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCgg                                                                           >  1:2721922/1‑150 (MQ=255)

                                                                           aCTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGt                                                                                                                              <  1:4218092/98‑1 (MQ=255)

                                                                           aCTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGt                                                                                                                              >  2:4218092/1‑98 (MQ=255)

                                                                              tGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTAc                                                                       <  1:4167843/150‑1 (MQ=255)

                                                                               gACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTGTCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCCAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTAcc                                                                      >  1:2199737/1‑150 (MQ=255)

                                                                                  aaGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGAc                                                                   >  2:4284776/1‑150 (MQ=255)

                                                                                    gTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACAt                                                                 >  2:1609184/1‑150 (MQ=255)

                                                                                    gTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACAt                                                                 >  2:1911686/1‑150 (MQ=255)

                                                                                       cAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATcgc                                                              >  1:4830030/1‑150 (MQ=255)

                                                                                           gggTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGTTGTGGCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGaaa                                                          <  2:2699814/150‑1 (MQ=255)

                                                                                           gggTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGaaa                                                          <  2:4694161/150‑1 (MQ=255)

                                                                                           gggTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGaaa                                                          <  1:4354499/150‑1 (MQ=255)

                                                                                              ttCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGcc                                                       >  1:4568207/1‑150 (MQ=255)

                                                                                                 aCCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATc                                                    >  1:765481/1‑150 (MQ=255)

                                                                                                 aCCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATc                                                    <  2:2256368/150‑1 (MQ=255)

                                                                                                  cccGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATTAAGATTTTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATcc                                                   >  2:566244/1‑150 (MQ=255)

                                                                                                  cccGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATcc                                                   >  1:1377260/1‑150 (MQ=255)

                                                                                                     ggATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCACGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACg                                                >  1:4258735/1‑150 (MQ=255)

                                                                                                      gATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGc                                               >  1:1902252/1‑150 (MQ=255)

                                                                                                      gATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGc                                               <  1:2172598/150‑1 (MQ=255)

                                                                                                           aCAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCAcc                                          <  2:1434507/150‑1 (MQ=255)

                                                                                                           aCAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCAcc                                          >  2:1729677/1‑150 (MQ=255)

                                                                                                           aCAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCAcc                                          <  1:2999398/150‑1 (MQ=255)

                                                                                                           aCAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCAcc                                          <  2:52804/150‑1 (MQ=255)

                                                                                                           aCAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCAcc                                          <  1:1151602/150‑1 (MQ=255)

                                                                                                             aaCCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCt                                        >  2:4224740/1‑150 (MQ=255)

                                                                                                               ccAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGa                                      >  1:976257/1‑150 (MQ=255)

                                                                                                               ccAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGa                                      <  2:1874114/150‑1 (MQ=255)

                                                                                                               ccAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGa                                      >  1:1086959/1‑150 (MQ=255)

                                                                                                                  gCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTgtctt                                   >  1:1771913/1‑146 (MQ=255)

                                                                                                                  gCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGacac                                   >  1:5446291/1‑150 (MQ=255)

                                                                                                                      atggatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCgtg                               >  1:1437645/1‑150 (MQ=255)

                                                                                                                      atggatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCgtg                               >  1:3722475/1‑150 (MQ=255)

                                                                                                                      atggatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCgtg                               >  1:747798/1‑150 (MQ=255)

                                                                                                                      atggatggCCTGCTGGCTATCGAACTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACAACGCGAAAGCCATCCACGCTCACCCTGACACCgtg                               >  2:2676900/1‑150 (MQ=255)

                                                                                                                       tggatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCgtgt                              >  1:1872806/1‑150 (MQ=255)

                                                                                                                        ggatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTc                             <  1:4296988/150‑1 (MQ=255)

                                                                                                                         gatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTcc                            >  1:4875788/1‑150 (MQ=255)

                                                                                                                          atggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCgg                                                                           <  1:3230736/102‑1 (MQ=255)

                                                                                                                           tggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCAGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCaa                          <  1:2698811/150‑1 (MQ=255)

                                                                                                                            ggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAg                         >  1:4362612/1‑150 (MQ=255)

                                                                                                                             gCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAgg                        >  2:777987/1‑150 (MQ=255)

                                                                                                                              ccTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGc                       <  2:4468210/150‑1 (MQ=255)

                                                                                                                              ccTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGc                       >  2:1055106/1‑150 (MQ=255)

                                                                                                                               ctgctgGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCg                      >  2:2789494/1‑150 (MQ=255)

                                                                                                                                     gCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTcc                >  1:1548895/1‑150 (MQ=255)

                                                                                                                                      cTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCa               <  2:2930733/150‑1 (MQ=255)

                                                                                                                                      cTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCa               <  1:1729677/150‑1 (MQ=255)

                                                                                                                                        aTCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGc             <  2:4041184/150‑1 (MQ=255)

                                                                                                                                         ttcaTCTCCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGcc            >  2:4486948/4‑150 (MQ=255)

                                                                                                                                         ttcaTCTCCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGc             >  2:4784899/4‑149 (MQ=255)

                                                                                                                                         ttcaTCTCCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGc             <  1:4784843/146‑1 (MQ=255)

                                                                                                                                          tcaTCTCCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCct           >  2:2583684/3‑150 (MQ=255)

                                                                                                                                           caTCTCCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCctg          >  2:952196/2‑150 (MQ=255)

                                                                                                                                            aTCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGc             >  1:226971/1‑146 (MQ=255)

                                                                                                                                            aTCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCctgc         >  2:3755358/1‑150 (MQ=255)

                                                                                                                                             tctcCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCctgct        >  1:3226230/1‑150 (MQ=255)

                                                                                                                                             tctcCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCctgct        <  2:3537670/150‑1 (MQ=255)

                                                                                                                                              ctcCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCctgctg       >  2:5004763/1‑150 (MQ=255)

                                                                                                                                                cccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCg     >  1:3763343/1‑150 (MQ=255)

                                                                                                                                                 ccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGt    >  1:4145715/1‑150 (MQ=255)

                                                                                                                                                 ccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGt    >  2:3427583/1‑150 (MQ=255)

                                                                                                                                                 ccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGt    >  1:5411371/1‑150 (MQ=255)

ACGAAAAAAAACCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCG  >  NC_002947/1750626‑1750922

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