BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,008,433	T→G	6.3%	F1017V (TTC→GTC)	PP_2627 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,008,433	0	T	G	6.3%	518.4 / 18.8	158	F1017V (TTC→GTC)	PP_2627	membrane protein
Reads supporting (aligned to +/- strand): ref base T (60/88); new base G (8/2); total (68/90)
Fisher's exact test for biased strand distribution p-value = 2.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

Only 100 of 159 total aligned reads displayed.

GCTGACGTGTACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCA  >  NC_002947/3008286‑3008580

  tGACGTGTACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCt                                                                                                                                                 <  2:1092526/150‑1 (MQ=255)

    aCGTGTACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGa                                                                                                                                               <  2:2576379/150‑1 (MQ=255)

    aCGTGTACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGa                                                                                                                                               <  1:134771/150‑1 (MQ=255)

     cGTGTACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGaa                                                                                                                                              <  2:3758999/150‑1 (MQ=255)

      gtgtACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAg                                                                                                                                             >  1:3642721/1‑150 (MQ=255)

         tACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGAc                                                                                                                                          >  2:3862852/1‑150 (MQ=255)

          aCAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACa                                                                                                                                         >  2:3246115/1‑150 (MQ=255)

             gTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAAcc                                                                                                                                      <  2:1357031/150‑1 (MQ=255)

              ttttttCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCt                                                                                                                                     >  2:1468569/1‑150 (MQ=255)

               tttttCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCtt                                                                                                                                    <  2:2304853/150‑1 (MQ=255)

                 tttCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGa                                                                                                                                  <  2:2167385/150‑1 (MQ=255)

                 tttCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGa                                                                                                                                  <  1:938937/150‑1 (MQ=255)

                 tttCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGa                                                                                                                                  <  1:3681561/150‑1 (MQ=255)

                 tttCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGa                                                                                                                                  <  1:2967997/150‑1 (MQ=255)

                  ttCCAGCCACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGCTGCTTCTTTGGAGGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAg                                                                                                                                 <  1:5120502/150‑1 (MQ=255)

                   tCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAgg                                                                                                                                <  2:4783634/150‑1 (MQ=255)

                    ccAGCAACGCCTGGCAGGCCGCTACACCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGc                                                                                                                               >  1:3375769/1‑150 (MQ=255)

                      aGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCt                                                                                                                             >  2:4152687/1‑150 (MQ=255)

                      aGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCt                                                                                                                             >  1:4184657/1‑150 (MQ=255)

                      aGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCt                                                                                                                             <  2:3256944/150‑1 (MQ=255)

                      aGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCt                                                                                                                             >  2:4994288/1‑150 (MQ=255)

                       gCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCtt                                                                                                                            <  1:876300/150‑1 (MQ=255)

                        cAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTg                                                                                                                           <  2:4731491/150‑1 (MQ=255)

                        cAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTg                                                                                                                           <  1:4897543/150‑1 (MQ=255)

                        cAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTg                                                                                                                           <  1:4599476/150‑1 (MQ=255)

                          aCGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCa                                                                                                                         <  1:151760/150‑1 (MQ=255)

                                    ggCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAg                                                                                                               >  2:2588062/1‑150 (MQ=255)

                                       cGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCAc                                                                                                            >  1:3465949/1‑150 (MQ=255)

                                          tACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACggg                                                                                                         <  2:2665943/150‑1 (MQ=255)

                                            ccccTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGTACGGGCa                                                                                                       >  1:470129/1‑150 (MQ=255)

                                            ccccTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGTACGGGCa                                                                                                       >  1:470110/1‑150 (MQ=255)

                                            ccccTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCa                                                                                                       >  1:1150704/1‑150 (MQ=255)

                                             cccTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCaa                                                                                                      >  1:1422248/1‑150 (MQ=255)

                                               cTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATc                                                                                                    >  2:5089761/1‑150 (MQ=255)

                                                   gtggtgAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTg                                                                                                <  2:316838/150‑1 (MQ=255)

                                                    tggtgAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGa                                                                                               <  2:4075393/150‑1 (MQ=255)

                                                      gtgAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGAtt                                                                                             >  1:1421706/1‑150 (MQ=255)

                                                      gtgAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGAtt                                                                                             >  1:1724087/1‑150 (MQ=255)

                                                      gtgAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGAtt                                                                                             >  1:289139/1‑150 (MQ=255)

                                                      gtgAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTGTCCTTCAGGAAAACCTTGAg                                                                                                                                 <  1:4818988/114‑1 (MQ=39)

                                                      gtgAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTGTCCTTCAGGAAAACCTTGAg                                                                                                                                 >  2:4818988/1‑114 (MQ=39)

                                                       tgAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCTGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTc                                                                                            >  1:4265813/1‑150 (MQ=255)

                                                       tgAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCTGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTc                                                                                            >  1:1518819/1‑150 (MQ=255)

                                                              aCCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGa                                                                                     <  2:1229066/150‑1 (MQ=255)

                                                              aCCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGa                                                                                     <  2:1229071/150‑1 (MQ=255)

                                                              aCCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGa                                                                                     >  2:5284366/1‑150 (MQ=255)

                                                                  gATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTg                                                                                 >  1:1470876/1‑150 (MQ=255)

                                                                     gCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATc                                                                              <  1:3167192/150‑1 (MQ=255)

                                                                     gCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATc                                                                              <  1:3411715/150‑1 (MQ=255)

                                                                       ttcttTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGa                                                                            <  1:626386/150‑1 (MQ=255)

                                                                         cttTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGagc                                                                          >  1:4067351/1‑150 (MQ=255)

                                                                            tGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGagcagc                                                                       <  2:650182/150‑1 (MQ=255)

                                                                              gATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTc                                                                     >  1:1953574/1‑150 (MQ=255)

                                                                                   cTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCt                                                                <  2:3737033/150‑1 (MQ=255)

                                                                                    ttCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATc                                                                              >  2:383885/1‑135 (MQ=255)

                                                                                    ttCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATc                                                                              <  1:383885/135‑1 (MQ=255)

                                                                                    ttCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTg                                                               >  1:4258517/1‑150 (MQ=255)

                                                                                    ttCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTg                                                               >  2:1615911/1‑150 (MQ=255)

                                                                                         ggCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGa                                                          <  2:2289429/150‑1 (MQ=255)

                                                                                         ggCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGa                                                          <  2:1601127/150‑1 (MQ=255)

                                                                                         ggCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGa                                                          <  2:1563381/150‑1 (MQ=255)

                                                                                           cGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAgc                                                        <  2:814333/150‑1 (MQ=255)

                                                                                                ttcGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATc                                                   >  2:3187495/1‑150 (MQ=255)

                                                                                                ttcGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATc                                                   >  1:5275858/1‑150 (MQ=255)

                                                                                                   gggCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGaa                                                                                                                                              <  2:1985737/56‑1 (MQ=255)

                                                                                                   gggCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGaa                                                                                                                                              >  1:1985737/1‑56 (MQ=255)

                                                                                                   gggCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCcgc                                                <  2:2859452/150‑1 (MQ=255)

                                                                                                   gggCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCcgc                                                >  2:3074029/1‑150 (MQ=255)

                                                                                                     gCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGa                                              <  2:433676/150‑1 (MQ=255)

                                                                                                      ccGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGaa                                             >  2:4597316/1‑150 (MQ=255)

                                                                                                      ccGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGaa                                             <  1:3110482/150‑1 (MQ=255)

                                                                                                      ccGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGaa                                             <  1:5358351/150‑1 (MQ=255)

                                                                                                         aaaGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGt                                                                                                                       >  2:5256374/1‑73 (MQ=255)

                                                                                                         aaaGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGt                                                                                                                       <  1:5256374/73‑1 (MQ=255)

                                                                                                               cgcCTGCAGCAGTTTAATGATCGGGATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCttcttc                                    <  2:5117065/150‑1 (MQ=255)

                                                                                                                         agTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTg                          <  1:1478271/150‑1 (MQ=255)

                                                                                                                          gTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTgg                         >  1:1232420/1‑150 (MQ=255)

                                                                                                                             tAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCaa                      <  1:1338606/150‑1 (MQ=255)

                                                                                                                              aaTGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAg                     >  2:4984230/1‑150 (MQ=255)

                                                                                                                                  aTCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTaa                 <  1:146627/150‑1 (MQ=255)

                                                                                                                                  aTCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTaa                 <  2:4829444/150‑1 (MQ=255)

                                                                                                                                    cGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAgt               <  2:648005/150‑1 (MQ=255)

                                                                                                                                      gTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAgtgt             >  1:2645206/1‑150 (MQ=255)

                                                                                                                                       tATCTCCAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGtgtg            <  1:4301915/150‑1 (MQ=255)

                                                                                                                                       tATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGtgtg            <  2:2402634/150‑1 (MQ=255)

                                                                                                                                       tATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGtgtg            <  2:3467984/150‑1 (MQ=255)

                                                                                                                                       tATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGtgtg            <  2:1724087/150‑1 (MQ=255)

                                                                                                                                         tctcAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCa          <  2:496695/150‑1 (MQ=255)

                                                                                                                                          ctcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAg         >  2:4062105/1‑150 (MQ=255)

                                                                                                                                           tcAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGt        >  2:5096562/1‑150 (MQ=255)

                                                                                                                                           tcAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGt        >  2:2439511/1‑150 (MQ=255)

                                                                                                                                           tcAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGCGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGt        >  2:4355541/1‑150 (MQ=255)

                                                                                                                                           tcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGt        >  1:2331835/1‑150 (MQ=255)

                                                                                                                                            cAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGtt       >  1:5376760/1‑150 (MQ=255)

                                                                                                                                            cAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGtt       <  2:2897082/150‑1 (MQ=255)

                                                                                                                                             aaGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACgg                                                                                                          <  1:645168/50‑1 (MQ=255)

                                                                                                                                             aaGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACgg                                                                                                          >  2:645168/1‑50 (MQ=255)

                                                                                                                                              aGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCa     <  1:2976846/150‑1 (MQ=255)

                                                                                                                                              aGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCa     <  1:3862852/150‑1 (MQ=255)

                                                                                                                                               ggTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATccc                                                                                                  >  2:542375/1‑56 (MQ=255)

GCTGACGTGTACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCA  >  NC_002947/3008286‑3008580

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