BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,413,440	T→A	7.3%	N155Y (AAC→TAC) ‡	PP_1235 ←	AhpC/TSA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,413,440	0	T	A	7.3%	349.3 / 14.8	109	N155Y (AAC→TAC) ‡	PP_1235	AhpC/TSA family protein
Reads supporting (aligned to +/- strand): ref base T (47/54); new base A (5/3); total (52/57)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

Only 100 of 109 total aligned reads displayed.

GGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCAGC  >  NC_002947/1413292‑1413585

gggTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTc                                                                                                                                                  >  2:415038/1‑150 (MQ=255)

  tttgtgtgggATATTCTACGCCATCCCGCTTGAGACGCTGATCAAGGACGTGCTGGATGCGTGGTCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAAGGTAATTGGGGCTGCAAGGCTGCCCCATTCAATTTCAAGCCTTGTTCaa                                                                                                                                                <  2:777209/143‑1 (MQ=255)

   ttttGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAg                                                                                                                                               <  2:907237/150‑1 (MQ=255)

    tttGGGGGATATTCTTCGCCATCCCGTTGGCGGTTCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAgg                                                                                                                                              <  2:193927/150‑1 (MQ=255)

       gggggATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGAACAAGGCtt                                                                                                                                           <  1:2375693/150‑1 (MQ=255)

        ggggATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTg                                                                                                                                          >  1:727817/1‑150 (MQ=255)

         gggATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTgg                                                                                                                                         <  1:1512599/150‑1 (MQ=255)

          ggATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTggg                                                                                                                                        >  2:1448109/1‑150 (MQ=255)

              atTCTTCGCCATCCCGCTGGCGACGCTGATCAAAGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCAACGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGggcggc                                                                                                                                    >  1:1770987/1‑150 (MQ=255)

                tcttcGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGgccg                                                                                                                                  >  2:632074/1‑150 (MQ=255)

                  ttcGCCATCCCGCTGGCGACGCTGATCNAGGCCGTGCTGGNTGCGTGNCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGgccgcc                                                                                                                                >  1:1421069/1‑150 (MQ=255)

                  ttcGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGgccgcc                                                                                                                                >  2:1831314/1‑150 (MQ=255)

                      ccATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATa                                                                                                                            <  1:2408790/150‑1 (MQ=255)

                      ccATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATa                                                                                                                            <  1:632074/150‑1 (MQ=255)

                            cGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCaa                                                                                                                                                <  2:1071876/124‑1 (MQ=255)

                            cGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCaa                                                                                                                                                >  1:1071876/1‑124 (MQ=255)

                             gCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATc                                                                                                                     <  2:1884811/150‑1 (MQ=255)

                                ggCGACGCTGATCAAGGCCGTGCTGGATGCGCGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCcac                                                                                                                  <  1:995040/150‑1 (MQ=255)

                                  cGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACAt                                                                                                                <  1:1170458/150‑1 (MQ=255)

                                  cGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACAt                                                                                                                <  1:1831314/150‑1 (MQ=255)

                                      gCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGccc                                                                                                            <  2:1740577/150‑1 (MQ=255)

                                        tGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCgg                                                                                                          <  1:109389/150‑1 (MQ=255)

                                          aTCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGc                                                                                                        >  2:553395/1‑150 (MQ=255)

                                           tCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCa                                                                                                       <  1:1266507/150‑1 (MQ=255)

                                              aGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGTCCGCCAATACGGCATCCACATGCCCGGGCACtt                                                                                                    >  1:1413949/1‑150 (MQ=255)

                                                 ccGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTca                                                                                                 <  1:415038/150‑1 (MQ=255)

                                                  agagatGGATGGGGGGCCTAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCGAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTcac                                                                                                <  2:1341474/145‑1 (MQ=255)

                                                     gCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTcacgcc                                                                                             >  2:2400602/1‑150 (MQ=255)

                                                      cTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTcacgcca                                                                                            <  1:1335872/150‑1 (MQ=255)

                                                              gTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTc                                                                                    >  2:722162/1‑150 (MQ=255)

                                                                   ccAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTGAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGgc                                                                               >  1:1958152/1‑150 (MQ=255)

                                                                   ccAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGgc                                                                               >  1:1519202/1‑150 (MQ=255)

                                                                   ccAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGgc                                                                               <  1:2172835/150‑1 (MQ=255)

                                                                   ccAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCAAATACGGCATCCACAGGCCCGGGCACTTTCACGCCACGACCTTCATGgc                                                                               >  2:2201379/1‑150 (MQ=255)

                                                                   ccAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGATGAAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCATGGGCGAGCGCGAATACGGCATCCACATGCCCGATCACTTTCCCGACATGCAAtaacaaaa                                                                               >  1:593766/1‑143 (MQ=255)

                                                                     aaGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACCTGCCCGGGCACTTTCACGCCACGCCATTCCTGgcgc                                                                             <  1:1448109/150‑1 (MQ=255)

                                                                        cAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCgcagc                                                                          >  2:427233/1‑150 (MQ=255)

                                                                          ggAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCAc                                                                        >  1:997805/1‑150 (MQ=255)

                                                                           gAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTgg                                                                                                                                         <  2:1833926/84‑1 (MQ=255)

                                                                           gAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTgg                                                                                                                                         >  1:1833926/1‑84 (MQ=255)

                                                                           gAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCt                                                                                  <  1:189196/139‑1 (MQ=255)

                                                                           gAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCt                                                                                  >  2:189196/1‑139 (MQ=255)

                                                                              ccTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATcc                                                                                                                    >  1:1435347/1‑102 (MQ=255)

                                                                              ccTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATcc                                                                                                                    <  2:1435347/102‑1 (MQ=255)

                                                                                tAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTc                                                                  >  2:509035/1‑150 (MQ=255)

                                                                                 aGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTcc                                                                 <  1:894539/150‑1 (MQ=255)

                                                                                       tCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCtt                                                                                                                                           <  2:1822146/70‑1 (MQ=255)

                                                                                       tCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCtt                                                                                                                                           >  1:1822146/1‑70 (MQ=255)

                                                                                        cACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTcc                                                                                   <  1:435519/125‑1 (MQ=255)

                                                                                           cGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCa                                                       <  2:588948/150‑1 (MQ=255)

                                                                                           cGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCa                                                       >  1:588948/1‑150 (MQ=255)

                                                                                           cGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCa                                                       <  1:277659/150‑1 (MQ=255)

                                                                                           cGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCa                                                       <  1:330602/150‑1 (MQ=255)

                                                                                                 tttGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACg                                                                                                                          <  2:821944/77‑1 (MQ=255)

                                                                                                 tttGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACg                                                                                                                          >  1:821944/1‑77 (MQ=255)

                                                                                                 tttGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGt                                                             <  2:210144/138‑1 (MQ=255)

                                                                                                 tttGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGt                                                             >  1:210144/1‑138 (MQ=255)

                                                                                                 tttGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAgg                                                 >  1:270156/1‑150 (MQ=255)

                                                                                                 tttGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAgg                                                 >  1:1808126/1‑150 (MQ=255)

                                                                                                     aaTGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGtgct                                             >  1:36405/1‑150 (MQ=255)

                                                                                                     aaTGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGtgct                                             >  2:412176/1‑150 (MQ=255)

                                                                                                      aTGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGtgctg                                            <  1:722162/150‑1 (MQ=255)

                                                                                                       tgtgATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGtgctgc                                           <  1:310574/150‑1 (MQ=255)

                                                                                                       tGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTc                                                                  >  1:36500/1‑127 (MQ=255)

                                                                                                       tGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTc                                                                  <  2:36500/127‑1 (MQ=255)

                                                                                                                  cTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGccc                                <  2:1808126/150‑1 (MQ=255)

                                                                                                                   tGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCa                               >  2:659365/1‑150 (MQ=255)

                                                                                                                    gCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGTCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCAt                              <  2:1413949/150‑1 (MQ=255)

                                                                                                                        ggCAGCACGCTTCCACTTCTAGCTTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTAc                          <  2:1268058/150‑1 (MQ=255)

                                                                                                                          cAGACCCATACAATGTTAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCCCGCCACGCCATTCCTAGCGCCGCTCCCCGTCCTTGTCGATCAGGAACGTACTGCGTTCAACGCCCCTGTCCTc                        <  2:945515/150‑1 (MQ=255)

                                                                                                                             ccccATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCtt                     >  1:1668997/1‑150 (MQ=255)

                                                                                                                             ccccATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCtt                     >  1:16858/1‑150 (MQ=255)

                                                                                                                               ccATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGc                   <  2:1828531/150‑1 (MQ=255)

                                                                                                                                 aTTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCa                 <  2:1668997/150‑1 (MQ=255)

                                                                                                                                 aTTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCa                 <  1:1842395/150‑1 (MQ=255)

                                                                                                                                 aTTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCa                 <  2:16858/150‑1 (MQ=255)

                                                                                                                                  ttCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGt                                                >  1:1856405/1‑118 (MQ=255)

                                                                                                                                  ttCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGt                                                <  2:1856405/118‑1 (MQ=255)

                                                                                                                                  ttCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCAt                >  2:842361/1‑150 (MQ=255)

                                                                                                                                  ttCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGAAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCAt                >  1:1306425/1‑150 (MQ=255)

                                                                                                                                     aaTTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCaa                                     <  2:1213888/126‑1 (MQ=255)

                                                                                                                                     aaTTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCaa                                     >  1:1213888/1‑126 (MQ=255)

                                                                                                                                     aaTTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACa             >  1:1334543/1‑150 (MQ=255)

                                                                                                                                      aTTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCCTCGACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAg            >  1:791988/1‑150 (MQ=255)

                                                                                                                                      aTTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCATGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAg            <  2:1421069/150‑1 (MQ=255)

                                                                                                                                       tttCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCtt                                                               <  1:1924967/98‑1 (MQ=255)

                                                                                                                                       tttCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGc           <  1:659365/150‑1 (MQ=255)

                                                                                                                                       tttCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGc           <  2:1334543/150‑1 (MQ=255)

                                                                                                                                        cccAAGCCTTGAACAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGtgctgc                                           <  1:248200/115‑1 (MQ=255)

                                                                                                                                        cccAAGCCTTGAACAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGtgctgc                                           >  2:248200/3‑117 (MQ=255)

                                                                                                                                         ccAAGCCTTGAACAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGctt         >  1:1047217/2‑150 (MQ=255)

                                                                                                                                         ccAAGCCTTGAACAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGctt         >  2:1445758/2‑150 (MQ=255)

                                                                                                                                          cAAGCCTTGAACAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGcttc        <  1:1445758/150‑1 (MQ=255)

                                                                                                                                          cAAGCCTTGAACAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGcttc        >  2:434961/1‑150 (MQ=255)

                                                                                                                                          cAAGCCTTGAACAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGcttc        >  1:1200644/1‑150 (MQ=255)

                                                                                                                                              ccTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTg                                                                                 >  1:1068170/1‑73 (MQ=255)

                                                                                                                                               cTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGct          >  1:1219833/1‑143 (MQ=255)

                                                                                                                                               cTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGct          <  2:1219833/143‑1 (MQ=255)

                                                                                                                                                ttGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCAGc  >  1:1799432/1‑150 (MQ=255)

                                                                                                                                                ttGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCAGc  >  2:2238708/1‑150 (MQ=255)

GGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCAGC  >  NC_002947/1413292‑1413585

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