BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,523,472	A→C	5.1%	F208C (TTC→TGC)	yjhG ←	putative dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,523,472	0	A	C	5.1%	445.5 / 9.4	149	F208C (TTC→TGC)	yjhG	putative dehydratase
Reads supporting (aligned to +/- strand): ref base A (72/69); new base C (2/6); total (74/75)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 151 total aligned reads displayed.

GGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCGCAACTCGCCACA  >  NC_000913/4523323‑4523621

gggCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGa                                                                                                                                                       <  2:1425294/150‑1 (MQ=255)

     gAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGt                                                                                                                                                  >  1:194840/1‑150 (MQ=255)

     gAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGt                                                                                                                                                  <  2:1302866/150‑1 (MQ=255)

     gAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGt                                                                                                                                                  >  1:330549/1‑150 (MQ=255)

       aTGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTgc                                                                                                                                                <  2:1070318/150‑1 (MQ=255)

          tgGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGcc                                                                                                                                             <  2:1377731/150‑1 (MQ=255)

          tgGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGcc                                                                                                                                             <  1:1680755/150‑1 (MQ=255)

          ggggATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGcc                                                                                                                                             <  1:1674339/149‑1 (MQ=255)

            ggATTGCCATTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCTCCGGAAGATGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCaa                                                                                                                                           >  2:943263/1‑150 (MQ=255)

            ggATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACCGCGTGCGTCATGTAGAGAAAATTCTCCATTGGCGAAGCGAGCGCCaa                                                                                                                                           >  2:813301/1‑150 (MQ=255)

               ttGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGCACGCTTCGCCAATgg                                                                                                                                        <  1:460500/150‑1 (MQ=255)

                tcccAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGt                                                                                                                                       <  1:11648/148‑1 (MQ=255)

                 gCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGCACGCTTCGCCAATggcg                                                                                                                                      <  1:1003970/150‑3 (MQ=255)

                  ccAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGttt                                                                                                                                     <  1:5290/150‑1 (MQ=255)

                   cAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTg                                                                                                                                    >  1:1657517/1‑150 (MQ=255)

                    aGTCCCAATCCTTCGGCCACCACCTGAGATGGCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGc                                                                                                                                   <  2:806560/150‑1 (MQ=255)

                     gTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCTCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCa                                                                                                                                  <  2:364583/150‑1 (MQ=255)

                          aaTCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTg                                                                                                                             >  1:1513633/1‑150 (MQ=255)

                            tCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGcc                                                                                                                           >  1:450168/1‑150 (MQ=255)

                              cTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGt                                                                                                                         >  2:778954/1‑150 (MQ=255)

                                tCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTg                                                                                                                       <  2:1535419/150‑1 (MQ=255)

                                  ggCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGtc                                                                                                                     >  1:1418977/1‑150 (MQ=255)

                                   gCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGtct                                                                                                                    <  1:928741/150‑1 (MQ=255)

                                   gCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGtct                                                                                                                    <  2:563084/150‑1 (MQ=255)

                                    ccaccaCCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGtctt                                                                                                                   >  1:923111/1‑150 (MQ=255)

                                    ccaccaCCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGtctt                                                                                                                   >  1:502752/1‑150 (MQ=255)

                                     caccacCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGt                                                                                                                         >  2:904120/1‑143 (MQ=255)

                                     caccacCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGt                                                                                                                         <  1:904120/143‑1 (MQ=255)

                                     caccacCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGtcttc                                                                                                                  >  2:475376/1‑150 (MQ=255)

                                       ccaccTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTtctc                                                                                                                <  2:1259259/150‑1 (MQ=255)

                                         accTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCa                                                                                                              >  1:78181/1‑150 (MQ=255)

                                          ccTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCAt                                                                                                             >  1:1326280/1‑150 (MQ=255)

                                                   tCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCgg                                                                                                    <  2:653298/150‑1 (MQ=255)

                                                    cccGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCTGTTGTCTTCTCCATCCTTTGCggg                                                                                                   >  2:1097483/1‑150 (MQ=255)

                                                         cAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCg                                                                                              <  2:1195439/150‑1 (MQ=255)

                                                           gTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTc                                                                                            <  1:1433204/150‑1 (MQ=255)

                                                            tGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTcg                                                                                           <  1:1705410/150‑1 (MQ=255)

                                                             gCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTcgc                                                                                          >  1:1120154/1‑150 (MQ=255)

                                                              cccAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTcgcg                                                                                         >  1:358638/1‑150 (MQ=255)

                                                              cccAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTcgcg                                                                                         >  1:513368/1‑150 (MQ=255)

                                                              cccAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTcgcg                                                                                         <  2:210261/150‑1 (MQ=255)

                                                              cccAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTcgcg                                                                                         >  2:112479/1‑150 (MQ=255)

                                                                cAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGTCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGcgcc                                                                                       <  1:73215/150‑1 (MQ=255)

                                                                cAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGcgcc                                                                                       >  1:298841/1‑150 (MQ=255)

                                                                  aaaaTTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGcgccgc                                                                                     <  2:1659210/150‑1 (MQ=255)

                                                                       tGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCgggg                                                                                <  2:194840/150‑1 (MQ=255)

                                                                       tGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCgggg                                                                                <  2:330549/150‑1 (MQ=255)

                                                                          cAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGAcc                                                                             >  2:523509/1‑150 (MQ=255)

                                                                          cAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGAcc                                                                             >  1:1664754/1‑150 (MQ=255)

                                                                            gccgccgccGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAg                                                                           >  2:1596603/1‑150 (MQ=255)

                                                                             ccgccgccgGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGTGTCGCGCCGCCGGGGACCAGc                                                                          >  2:1161596/1‑150 (MQ=255)

                                                                             ccgccgccgGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGc                                                                          <  1:112479/150‑1 (MQ=255)

                                                                              cgccgccgGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTAGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGcc                                                                         >  2:1141353/1‑149 (MQ=255)

                                                                                   ccgGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGtt                                                                    >  1:897587/1‑150 (MQ=255)

                                                                                   ccgGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGtt                                                                    >  2:1568747/1‑150 (MQ=255)

                                                                                     ggAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGc                                                                  >  1:878356/1‑150 (MQ=255)

                                                                                     ggAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGc                                                                  >  1:162409/1‑150 (MQ=255)

                                                                                      gAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCg                                                                 <  1:661831/150‑1 (MQ=255)

                                                                                       aaGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGa                                                                >  1:706971/1‑150 (MQ=255)

                                                                                         gagGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGAtg                                                              >  2:1186265/1‑150 (MQ=255)

                                                                                           ggCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGAtgtt                                                            <  1:1302087/150‑1 (MQ=255)

                                                                                            gCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCTCGCCGCCGGGGACCAGCACGGTTGCGAtgttg                                                           <  2:923111/150‑1 (MQ=255)

                                                                                             cacaGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGGCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGAtgttgt                                                          <  1:1237084/150‑1 (MQ=255)

                                                                                               caGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGtgc                                                        <  2:544663/150‑1 (MQ=255)

                                                                                                   cTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTgcg                                                    <  2:792544/150‑1 (MQ=255)

                                                                                                      tACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCgcgg                                                 <  1:1258102/150‑1 (MQ=255)

                                                                                                       aCAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCgcggc                                                >  1:1204551/1‑150 (MQ=255)

                                                                                                        caggacGCACGCCGTGCGGGCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCgcggcg                                               >  1:87459/6‑150 (MQ=255)

                                                                                                        caggacGCACGCCGTGCGGGCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCgcggcg                                               >  1:696500/6‑150 (MQ=255)

                                                                                                        caggacGCACGCCGTGCGGGCTGTAAAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCgcggcg                                               >  1:1475472/6‑150 (MQ=255)

                                                                                                          gCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAg                                             >  1:26865/1‑150 (MQ=255)

                                                                                                           cccGCACGGCGTGCGTCCTGTAGAGATAATTCGCCCTGGGCGAAGCGTGCGCCAATGGTGTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGGTGCGATGTTGTGCTGCGCGGCGAGt                                            <  1:102353/150‑1 (MQ=255)

                                                                                                                 cGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCTCGGCGAGTGCcatc                                      >  2:1255643/1‑150 (MQ=255)

                                                                                                                   gcgtgcgtCCTGTAGAGATCATACGCCATTGGCGAAGCGCGCGCCACTGGATTGCAACTTGCCGTTGGCTGCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGAGATGTTGTGCTGCGCGGCGAGTGCcttcat                                    <  1:1074251/150‑1 (MQ=255)

                                                                                                                     gtgcgtCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCAtgg                                  <  2:1322075/150‑1 (MQ=255)

                                                                                                                      tgcgtCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCAtggt                                 <  2:1204551/150‑1 (MQ=255)

                                                                                                                       gcgtCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCAtggtg                                <  1:219415/150‑1 (MQ=255)

                                                                                                                         gtCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGc                              >  2:681717/1‑150 (MQ=255)

                                                                                                                          tCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCTGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGcc                             <  1:1097483/150‑1 (MQ=255)

                                                                                                                           ccTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGCTGTTGTGCTGCGCGGCGAGTGCGATCATGGTGGCCg                            >  1:1566110/1‑150 (MQ=255)

                                                                                                                           ccTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCg                            <  2:1625891/150‑1 (MQ=255)

                                                                                                                               tAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAg                        <  2:1657517/150‑1 (MQ=255)

                                                                                                                                agagATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGc                       <  2:513368/150‑1 (MQ=255)

                                                                                                                                     tAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCtt                  <  2:162409/150‑1 (MQ=255)

                                                                                                                                      aaTTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTa                 <  1:1488315/150‑1 (MQ=255)

                                                                                                                                         tCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCg              >  2:402508/1‑150 (MQ=255)

                                                                                                                                         tCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCg              <  2:1504249/150‑1 (MQ=255)

                                                                                                                                           gCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCGCa            <  2:567411/150‑1 (MQ=255)

                                                                                                                                            ccATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCGCaa           >  2:798035/1‑150 (MQ=255)

                                                                                                                                             cATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCGCAAc          <  2:1353968/150‑1 (MQ=255)

                                                                                                                                             cATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCGCAAc          >  1:386405/1‑150 (MQ=255)

                                                                                                                                              aTTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGACGCCCCTTATCGCAACt         >  2:640077/1‑150 (MQ=255)

                                                                                                                                              aTTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCGCAACt         <  1:523509/150‑1 (MQ=255)

                                                                                                                                              aTTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCGCAACt         >  2:713293/1‑150 (MQ=255)

                                                                                                                                              aTTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCGCAACt         <  2:298841/150‑1 (MQ=255)

                                                                                                                                                 ggCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCGCAACTCGc      >  1:1625777/1‑150 (MQ=255)

                                                                                                                                                  gCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCGCAACTCGcc     >  1:974538/1‑150 (MQ=255)

                                                                                                                                                  gCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCAGCATGGTGGCCGGAAGCCCCTTATCGCAACTCGcc     >  1:777554/1‑150 (MQ=255)

                                                                                                                                                   ggAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCGCAACTCGCca    <  2:1418977/149‑1 (MQ=255)

                                                                                                                                                     aaGCGTGCGCCAATGGGTTGCACCTTGCCGTTGGCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCGCAACTCGCcaca  <  2:757951/150‑1 (MQ=255)

GGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCGCAACTCGCCACA  >  NC_000913/4523323‑4523621

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