BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,893,446	T→A	37.9%	K63* (AAA→TAA)	ygcW ←	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,893,446	0	T	A	37.9%	327.3 / inf	327	K63* (AAA→TAA)	ygcW	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (114/89); new base A (60/64); total (174/153)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 337 total aligned reads displayed.

CCAGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAA  >  CP009273/2893299‑2893595

 cAGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATAACCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTc                                                                                                                                                    <  2:305868/150‑1 (MQ=255)

  aGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCa                                                                                                                                                   >  2:4018486/1‑150 (MQ=255)

    cACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAAt                                                                                                                                                 <  2:103959/150‑1 (MQ=255)

     aCCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATc                                                                                                                                                >  1:2303505/1‑150 (MQ=255)

      ccTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCa                                                                                                                                               <  2:126552/150‑1 (MQ=255)

         tATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCTATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCAttt                                                                                                                                            <  2:4345765/150‑1 (MQ=255)

         tATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCAttt                                                                                                                                            <  1:4516399/150‑1 (MQ=255)

         tATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTAACGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCAttt                                                                                                                                            >  1:3371725/1‑150 (MQ=255)

            tCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCt                                                                                                                                         >  1:3997962/1‑150 (MQ=255)

            tCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCt                                                                                                                                         <  1:1405163/150‑1 (MQ=255)

            tCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCt                                                                                                                                         >  1:3434755/1‑150 (MQ=255)

             cAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATAACCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCtt                                                                                                                                        >  2:1599690/1‑150 (MQ=255)

              aGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCttt                                                                                                                                       >  2:2195236/1‑150 (MQ=255)

               gCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTg                                                                                                                                      >  1:1000397/1‑150 (MQ=255)

               gCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTg                                                                                                                                      <  2:1496180/150‑1 (MQ=255)

                cTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGt                                                                                                                                     <  2:1625370/150‑1 (MQ=255)

                 ttACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGtt                                                                                                                                    >  2:4139041/1‑150 (MQ=255)

                  tACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGttt                                                                                                                                   <  1:1972003/150‑1 (MQ=255)

                   aCAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTc                                                                                                                                  >  1:5077013/1‑150 (MQ=255)

                   aCAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTc                                                                                                                                  <  1:4639029/150‑1 (MQ=255)

                   aCAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTc                                                                                                                                  <  2:3219031/150‑1 (MQ=255)

                   aCAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTc                                                                                                                                  <  1:1190812/150‑1 (MQ=255)

                    cAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCg                                                                                                                                 <  1:491617/150‑1 (MQ=255)

                       aTACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCg                                                                                                                              >  1:1762073/1‑150 (MQ=255)

                        tACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGAGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGt                                                                                                                             >  2:4881300/1‑150 (MQ=255)

                         aCCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGtt                                                                                                                            <  1:3681813/150‑1 (MQ=255)

                           cGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTAt                                                                                                                          >  2:4241905/1‑150 (MQ=255)

                           cGGCATTGTTAACCAGAAGATCAACTGGACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTAt                                                                                                                          <  1:28052/150‑1 (MQ=255)

                            ggCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATc                                                                                                                         >  2:5140152/1‑150 (MQ=255)

                             gCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCt                                                                                                                        >  1:1702107/1‑150 (MQ=255)

                               aTTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCttt                                                                                                                      >  1:2172014/1‑150 (MQ=255)

                                  gTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGAc                                                                                                                   <  1:4438373/150‑1 (MQ=255)

                                   ttAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACg                                                                                                                  <  1:1662274/150‑1 (MQ=255)

                                     aaCCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGaa                                                                                                                >  1:3069797/1‑150 (MQ=255)

                                        cAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACt                                                                                                             >  2:1182969/1‑150 (MQ=255)

                                           aaTATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAgg                                                                                                          <  2:3758678/150‑1 (MQ=255)

                                            atatCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGa                                                                                                         >  2:955402/1‑150 (MQ=255)

                                              atCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAAt                                                                                                       >  1:1606335/1‑150 (MQ=255)

                                                cAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATaa                                                                                                     >  1:643145/1‑150 (MQ=255)

                                                 aaCTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATaaa                                                                                                    <  1:4034792/150‑1 (MQ=255)

                                                 aaCTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATaaa                                                                                                    >  1:4287667/1‑150 (MQ=255)

                                                  aCTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAg                                                                                                   <  1:3064363/150‑1 (MQ=255)

                                                    tGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGACACTCGGACTAAAGat                                                                                                 >  2:2830350/1‑150 (MQ=255)

                                                      tACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGatat                                                                                               <  2:3968355/150‑1 (MQ=255)

                                                          gAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTgc                                                                                           >  2:907422/1‑150 (MQ=255)

                                                          gAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTgc                                                                                           >  1:2260000/1‑150 (MQ=255)

                                                              cGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCa                                                                                       <  2:4575840/150‑1 (MQ=255)

                                                               gCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAg                                                                                      >  1:1331066/1‑150 (MQ=255)

                                                                cTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGc                                                                                     >  1:2524435/1‑150 (MQ=255)

                                                                cTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGc                                                                                     >  2:5032216/1‑150 (MQ=255)

                                                                 tCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCt                                                                                    >  1:919084/1‑150 (MQ=255)

                                                                  cacaGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCtt                                                                                   >  1:4778885/1‑150 (MQ=255)

                                                                  cacaGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCtt                                                                                   >  1:1089002/1‑150 (MQ=255)

                                                                  cacaGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCtt                                                                                   >  1:29766/1‑150 (MQ=255)

                                                                   acaGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCttt                                                                                  >  2:1875111/1‑150 (MQ=255)

                                                                   acaGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTGTGACGAAACTAGGAATAAAGATATTTGCGCCAGCtgt                                                                                  >  2:3177597/1‑148 (MQ=255)

                                                                    caGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTg                                                                                 >  1:3144580/1‑150 (MQ=255)

                                                                     agcaagcaGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTgg                                                                                >  1:695387/1‑150 (MQ=255)

                                                                      gcaagcaGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGc                                                                               <  1:3927012/150‑1 (MQ=255)

                                                                       caagcaGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGcc                                                                              <  2:5111407/150‑1 (MQ=255)

                                                                          gcagcGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAAc                                                                           <  1:4987565/150‑1 (MQ=255)

                                                                          gcagcGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAAc                                                                           <  1:2046102/150‑1 (MQ=255)

                                                                             gcGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGcc                                                                        >  1:4158535/1‑150 (MQ=255)

                                                                             gcGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGcc                                                                        <  2:2050091/150‑1 (MQ=255)

                                                                                aTAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATg                                                                     <  2:4835254/150‑1 (MQ=255)

                                                                                     cttctGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCaaa                                                                >  2:3628355/1‑150 (MQ=255)

                                                                                     cttctGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCaaa                                                                <  2:795265/150‑1 (MQ=255)

                                                                                       tctGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATg                                                              <  2:4689010/150‑1 (MQ=255)

                                                                                               gcgcCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCt                                                      >  1:2741/1‑150 (MQ=255)

                                                                                               gcgcCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTATGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCt                                                      >  1:2606/1‑150 (MQ=255)

                                                                                                  ccTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTaaa                                                   >  2:736694/1‑150 (MQ=255)

                                                                                                    ttCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAAACATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAAcc                                                 <  2:1257610/150‑1 (MQ=255)

                                                                                                    ttCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAAcc                                                 >  2:4991587/1‑150 (MQ=255)

                                                                                                     tCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCg                                                <  1:4393960/150‑1 (MQ=255)

                                                                                                      cTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGc                                               <  1:2899901/150‑1 (MQ=255)

                                                                                                         cGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTAt                                            >  2:2471144/1‑150 (MQ=255)

                                                                                                          ggTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCCATCCTTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTAtt                                           <  2:3465936/150‑1 (MQ=255)

                                                                                                             gATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTccc                                        >  1:3943024/1‑150 (MQ=255)

                                                                                                              aTACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCa                                       >  1:1519634/1‑150 (MQ=255)

                                                                                                              aTACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCa                                       <  1:562905/150‑1 (MQ=255)

                                                                                                                aCCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCAcc                                     <  2:3355455/150‑1 (MQ=255)

                                                                                                                aCCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCAcc                                     >  2:2286685/1‑150 (MQ=255)

                                                                                                                   cACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGt                                  >  2:2316783/1‑150 (MQ=255)

                                                                                                                     ccTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGTAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTaa                                <  1:1871307/150‑1 (MQ=255)

                                                                                                                     ccTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTaa                                <  2:1341887/150‑1 (MQ=255)

                                                                                                                           tGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTg                          >  1:3109930/1‑150 (MQ=255)

                                                                                                                              aGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCgg                       >  2:4804723/1‑150 (MQ=255)

                                                                                                                                 ccACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGttt                    <  1:3559301/150‑1 (MQ=255)

                                                                                                                                  cACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGtttt                   >  1:3032916/1‑150 (MQ=255)

                                                                                                                                     cTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAcc                >  2:3716604/1‑150 (MQ=255)

                                                                                                                                       cAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCtt              <  2:1787464/150‑1 (MQ=255)

                                                                                                                                         acacCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTc            >  1:1594530/1‑150 (MQ=255)

                                                                                                                                         acacCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTc            >  2:2217836/1‑150 (MQ=255)

                                                                                                                                         acacCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTc            <  1:5064418/150‑1 (MQ=255)

                                                                                                                                         acacCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTc            <  1:2521926/150‑1 (MQ=255)

                                                                                                                                           acCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAg          >  1:1825029/1‑150 (MQ=255)

                                                                                                                                             ccTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAggg        <  2:2646868/150‑1 (MQ=255)

                                                                                                                                             ccTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAggg        >  2:4914648/1‑150 (MQ=255)

                                                                                                                                              cTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGga       >  1:1431386/1‑150 (MQ=255)

                                                                                                                                                  ttttCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGaaa   <  1:3581816/150‑1 (MQ=255)

CCAGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAA  >  CP009273/2893299‑2893595

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