BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,962,561	G→A	7.6%	N68N (AAC→AAT) ‡	ppdC ←	putative YgdA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,962,561	0	G	A	7.6%	326.8 / 10.9	105	N68N (AAC→AAT) ‡	ppdC	putative YgdA
Reads supporting (aligned to +/- strand): ref base G (45/52); new base A (2/6); total (47/58)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 106 total aligned reads displayed.

TGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCG  >  NC_000913/2962412‑2962709

tGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAg                                                                                                                                                      >  1:1614241/1‑150 (MQ=255)

  gTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGtt                                                                                                                                                    >  1:606961/1‑150 (MQ=255)

  gTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGtt                                                                                                                                                    >  1:17366/1‑150 (MQ=255)

     gACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGc                                                                                                                                                 <  2:1082566/150‑1 (MQ=255)

     gACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGc                                                                                                                                                 <  1:743440/150‑1 (MQ=255)

     gACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGc                                                                                                                                                 <  2:665171/150‑1 (MQ=255)

     gACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGc                                                                                                                                                 <  1:698789/150‑1 (MQ=255)

       cccTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAg                                                                                                                                               >  1:152498/1‑150 (MQ=255)

       cccTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAg                                                                                                                                               >  1:182790/1‑150 (MQ=255)

          ttAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGtc                                                                                                                                            <  2:1128459/150‑2 (MQ=255)

           tAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTgg                                                                                                                                           <  2:1626408/150‑1 (MQ=255)

             aCATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCg                                                                                                                                         >  2:44192/1‑150 (MQ=255)

             aCATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGtcaac                                                                                                                                         <  1:1157574/150‑5 (MQ=255)

             aCATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGtcaac                                                                                                                                         <  2:661569/150‑5 (MQ=255)

              cATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGa                                                                                                                                        <  1:1326530/150‑1 (MQ=255)

              cATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGtcaacc                                                                                                                                        <  1:1032085/150‑6 (MQ=255)

              cATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGtcaacc                                                                                                                                        <  2:375626/150‑6 (MQ=255)

               aTAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGaa                                                                                                                                       <  1:1120539/150‑1 (MQ=255)

                tAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGaaa                                                                                                                                      >  1:1126331/1‑150 (MQ=255)

                 aGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAAt                                                                                                                                     <  2:1311397/150‑1 (MQ=255)

                  gCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATc                                                                                                                                    >  2:1004350/1‑150 (MQ=255)

                  gCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATc                                                                                                                                    >  2:1010277/1‑150 (MQ=255)

                    ggCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATcgc                                                                                                                                  <  1:730161/150‑1 (MQ=255)

                     gCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATcgcg                                                                                                                                 <  1:1234700/150‑1 (MQ=255)

                     gCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATcgcg                                                                                                                                 <  1:1228900/150‑1 (MQ=255)

                        ccTGACTACTGACGATTCGGGCAATGCAGGCGGGGCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCgcgcgc                                                                                                                              <  1:171676/150‑1 (MQ=255)

                          cgctagtttcacccgggggcagagaaggcgagatgacccgcctgcattgcccCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAg                                                                                                                            >  1:1117995‑M2/53‑150 (MQ=255)

                           gACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGt                                                                                                                           >  1:1256194/1‑150 (MQ=255)

                           gACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCGTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGt                                                                                                                           >  2:649782/1‑150 (MQ=255)

                               actGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCg                                                                                                                       <  2:921654/150‑1 (MQ=255)

                                    acccgggggcagagaaggcgagatgacccgcctgcattgcccCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTtgc                                                                                                                  >  2:331959‑M2/43‑150 (MQ=255)

                                       cgggggcagagaaggcgagatgacccgcctgcattgcccCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCTTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTtgctgc                                                                                                               >  2:808272‑M2/40‑150 (MQ=255)

                                         cGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgcca                                                                                                             >  2:497566/1‑150 (MQ=255)

                                          GGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTTCATTCTGTTGACCTGCCAGTTGGCATGTGGCTAAATCGCGCGTAGTTTCGTTTTCTgacag                                                                                                            >  2:1480822‑M1/1‑150 (MQ=255)

                                          GGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgccag                                                                                                            <  2:1435414‑M1/150‑1 (MQ=255)

                                               gagaaggcgagatgacccgccttctctgcccCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATg                                                                                                       >  1:523637‑M2/32‑150 (MQ=255)

                                                 gCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATgcc                                                                                                     >  2:1513167/1‑150 (MQ=255)

                                                     gCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCa                                                                                                 <  2:1127965/150‑1 (MQ=255)

                                                      cGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCag                                                                                                <  1:1214026/150‑1 (MQ=255)

                                                         ggcagctcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCGGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGGTGACCTGCCAGTGGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGc                                                                                             <  2:59315/145‑1 (MQ=255)

                                                         gTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGc                                                                                             <  1:1010277/150‑1 (MQ=255)

                                                         gTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGc                                                                                             <  1:1004350/150‑1 (MQ=255)

                                                             tctcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGtt                                                                                         >  2:611566/1‑150 (MQ=255)

                                                              ctcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTg                                                                                        >  2:221228/1‑150 (MQ=255)

                                                              ctcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTg                                                                                        >  2:173652/1‑150 (MQ=255)

                                                              ctcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTg                                                                                        <  1:243317/150‑1 (MQ=255)

                                                               tcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTgg                                                                                       >  1:352429/1‑150 (MQ=255)

                                                                cGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGt                                                                                      >  2:851209/1‑150 (MQ=255)

                                                                cGCCTTCTCTGCCCCCGGGGGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGt                                                                                      <  2:1249104/150‑1 (MQ=255)

                                                                  ccTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTAc                                                                                    >  2:1041246/1‑150 (MQ=255)

                                                                     tctcTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTgg                                                                                 <  1:946111/150‑1 (MQ=255)

                                                                     tctcTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTgg                                                                                 <  2:668599/150‑1 (MQ=255)

                                                                      ctctGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGt                                                                                <  2:1126331/150‑1 (MQ=255)

                                                                       tctGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGtt                                                                               <  2:292234/150‑1 (MQ=255)

                                                                        ctGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGttt                                                                              >  2:762000/1‑150 (MQ=255)

                                                                         tgcccCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGGAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTc                                                                             <  1:865464‑M2/145‑1 (MQ=255)

                                                                               cGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCTGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTgg                                                                       <  1:1533597/150‑1 (MQ=255)

                                                                               cGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTgg                                                                       >  2:1226078/1‑150 (MQ=255)

                                                                               cGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTgg                                                                       >  2:292823/1‑150 (MQ=255)

                                                                               cGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTgg                                                                       <  2:603280/150‑1 (MQ=255)

                                                                                gggTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGc                                                                      <  2:1483056/150‑1 (MQ=255)

                                                                                  gTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCaa                                                                    <  1:1513167/150‑1 (MQ=255)

                                                                                   tGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAg                                                                   <  1:851209/150‑1 (MQ=255)

                                                                                       aCTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTg                                                               <  1:1083588/150‑1 (MQ=255)

                                                                                          aGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTc                                                            <  2:473514/150‑1 (MQ=255)

                                                                                           gCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCa                                                           >  2:1121230/1‑150 (MQ=255)

                                                                                             gTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCAtt                                                         <  1:278826/150‑1 (MQ=255)

                                                                                             gTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCAtt                                                         >  1:706987/1‑150 (MQ=255)

                                                                                               aaCGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGTCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTaa                                                       <  1:633880/150‑1 (MQ=255)

                                                                                                    tGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTc                                                  <  1:1215221/150‑1 (MQ=255)

                                                                                                      aTGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGc                                                <  2:875350/150‑1 (MQ=255)

                                                                                                         cTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTgg                                             <  2:182790/150‑1 (MQ=255)

                                                                                                         cTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTgg                                             <  1:331959/150‑1 (MQ=255)

                                                                                                             cacaTCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAAc                                         <  2:17366/150‑1 (MQ=255)

                                                                                                              acaTCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAAcc                                        <  1:1112984/150‑1 (MQ=255)

                                                                                                                aTCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCg                                      <  2:606961/150‑1 (MQ=255)

                                                                                                                  cccGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGAt                                    >  1:283837/1‑150 (MQ=255)

                                                                                                                  cccGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGAt                                    <  1:1041246/150‑1 (MQ=255)

                                                                                                                       cTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGc                               >  2:1579259/1‑150 (MQ=255)

                                                                                                                        tGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCa                              >  1:1042450/1‑150 (MQ=255)

                                                                                                                         gCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAg                             <  1:1257951/150‑1 (MQ=255)

                                                                                                                           gATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTg                           <  1:997078/150‑1 (MQ=255)

                                                                                                                           gATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTg                           <  2:632762/150‑1 (MQ=255)

                                                                                                                             tGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGAc                         <  2:1117995/150‑1 (MQ=255)

                                                                                                                             tGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGAc                         <  1:497566/150‑1 (MQ=255)

                                                                                                                                cTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAat                      >  2:803991/1‑150 (MQ=255)

                                                                                                                                cTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAat                      >  2:1456637/1‑150 (MQ=255)

                                                                                                                                 tGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAatc                     <  1:1435414/150‑1 (MQ=255)

                                                                                                                                  gCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAatca                    >  2:541079/1‑150 (MQ=255)

                                                                                                                                    aTTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAatcatc                  <  2:1614241/150‑1 (MQ=255)

                                                                                                                                        ggTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGTCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGCCTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCa              >  1:882604/1‑150 (MQ=255)

                                                                                                                                         cgccaCCTGCCAACTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCAt             >  1:855058/5‑150 (MQ=255)

                                                                                                                                           tGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCa           >  1:1503547/1‑150 (MQ=255)

                                                                                                                                           tGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCa           >  2:436764/1‑150 (MQ=255)

                                                                                                                                            gACCTGCCAGTTGGAAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCaa          >  2:621939/1‑150 (MQ=255)

                                                                                                                                             aCCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGATACCCGATAACGCAGTGACAATCATCACCATCAac         >  2:146737/1‑150 (MQ=255)

                                                                                                                                             aCCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAac         >  2:1248833/1‑150 (MQ=255)

                                                                                                                                              ccTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAaca        <  1:762000/150‑1 (MQ=255)

                                                                                                                                                 gCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCa     >  2:1069979/1‑150 (MQ=255)

                                                                                                                                                    aGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCg  >  2:1394185/1‑150 (MQ=255)

TGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCG  >  NC_000913/2962412‑2962709

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