BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,995,716	A→G	5.1%	H45R (CAT→CGT)	yigA →	DUF484 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,995,716	0	A	G	5.1%	372.6 / 7.0	118	H45R (CAT→CGT)	yigA	DUF484 family protein
Reads supporting (aligned to +/- strand): ref base A (67/45); new base G (3/3); total (70/48)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-01

Only 100 of 119 total aligned reads displayed.

ACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCG  >  NC_000913/3995567‑3995863

acgGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCa                                                                                                                                                     >  2:547566/1‑150 (MQ=255)

  ggATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATc                                                                                                                                                   >  2:1089438/1‑150 (MQ=255)

   gATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATcc                                                                                                                                                  >  1:776283/1‑150 (MQ=255)

     tttATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGCCCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCgt                                                                                                                                                >  1:331585/1‑149 (MQ=255)

     tttATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCgg                                                                                                                                                >  2:1230807/1‑150 (MQ=255)

       tATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTa                                                                                                                                              >  2:76736/1‑150 (MQ=255)

            aTCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCgg                                                                                                                                         <  2:612360/150‑1 (MQ=255)

                 tGAAGCACCCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCGTACCGGATGCGGCAcac                                                                                                                                    <  1:998615/150‑3 (MQ=255)

                 tGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCg                                                                                                                                    >  1:321134/1‑150 (MQ=255)

                 tGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCg                                                                                                                                    >  2:594148/1‑150 (MQ=255)

                 tGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGGCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCg                                                                                                                                    <  1:1642098/150‑1 (MQ=255)

                 tGAAGCAACCAGCGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCACTAGAAGCGAAACGTCTGCCGCATCCGGTACGCGCCAccc                                                                                                                                    >  1:1327905/1‑149 (MQ=255)

                  gAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGt                                                                                                                                   >  1:593428/1‑150 (MQ=255)

                   aaGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCGTACCGGATGCGGCAcacgt                                                                                                                                  <  1:1709631/150‑5 (MQ=255)

                    aGCAACCAGGGGAGGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGttt                                                                                                                                 >  2:570352/1‑150 (MQ=255)

                    aGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGttt                                                                                                                                 >  2:1321291/1‑150 (MQ=255)

                     gCAACCCGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGCCCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTc                                                                                                                                <  1:1713529/150‑1 (MQ=255)

                          cAGGGGAAGAACTGCAGGCAACACTCAAGGAGCTTGATGACCGGGCGGTTGTCGATTATCTCATTAAAAATCCTGAGTTTTTTAACCTTAATGCGCGCGCAGTAGAAGCGATACGTGAGACGCATCCGGTACGCGGCACCGTTTCGTTgg                                                                                                                           >  1:1546535/1‑150 (MQ=255)

                          cAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTgg                                                                                                                           >  2:1697734/1‑150 (MQ=255)

                          cAGGGGAAGAACTGCAGGAAACACTCACGGAGCGTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTACCGCATCCGGTACGCGGCACCGTTTCGTTgg                                                                                                                           <  2:524593/150‑1 (MQ=255)

                                aagacCTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGGGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGt                                                                                                                     <  1:628853/150‑1 (MQ=255)

                                aactaCTGCAGGAAACACTCACGGAGGTTGAAGACCGGGCGGTTGTCGATTATCTGATTAAAACTCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGt                                                                                                                     <  1:1633502/146‑1 (MQ=255)

                                 agaaCTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTg                                                                                                                    >  2:1734696/1‑150 (MQ=255)

                                  gaaCTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTgg                                                                                                                   <  1:409896/150‑1 (MQ=255)

                                    aCTGCAGGAACCACTCACGGCGCTTGATGACCGGGCGGTTGTCGAGTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGca                                                                                                                 <  2:1006665/150‑1 (MQ=255)

                                      tGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGcaca                                                                                                               >  2:45398/1‑150 (MQ=255)

                                       gCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACAt                                                                                                              >  2:1894560/1‑150 (MQ=255)

                                       gCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACAt                                                                                                              >  1:230913/1‑150 (MQ=255)

                                       gCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACAt                                                                                                              >  2:783847/1‑150 (MQ=255)

                                        cAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATg                                                                                                             <  1:594148/150‑1 (MQ=255)

                                          ggAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGc                                                                                                           >  2:1749240/1‑150 (MQ=255)

                                               cacTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCcgcg                                                                                                      >  1:1496114/1‑150 (MQ=255)

                                                ccTCACGGAGCGTGATGACCGGGCGGGTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCgcgc                                                                                                     <  1:1653208/149‑1 (MQ=255)

                                                 cTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCTTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGTTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCa                                                                                                    >  1:169467/1‑150 (MQ=255)

                                                 cTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCa                                                                                                    >  1:837631/1‑150 (MQ=255)

                                                    aCGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGt                                                                                                 >  1:861909/1‑150 (MQ=255)

                                                      ggAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGAACATGGCCCGCTCACGTaa                                                                                               >  1:439158/1‑150 (MQ=255)

                                                       gAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAAt                                                                                              <  1:744332/150‑1 (MQ=255)

                                                       gAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAAt                                                                                              <  1:1614579/150‑1 (MQ=255)

                                                         gCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCa                                                                                            >  1:1008055/1‑150 (MQ=255)

                                                           ttGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCata                                                                                          >  1:22246/1‑150 (MQ=255)

                                                            tgatgaCCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCatat                                                                                         >  2:572682/1‑150 (MQ=255)

                                                             gatgaCCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATAtt                                                                                        <  2:1012521/150‑1 (MQ=255)

                                                              atgaCCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTc                                                                                       <  1:534220/150‑1 (MQ=255)

                                                                gaCCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCAt                                                                                     <  2:1531012/150‑1 (MQ=255)

                                                                      gCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTg                                                                               <  1:182598/150‑1 (MQ=255)

                                                                      gCGGTTGTCGAGTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTg                                                                               <  1:156541/150‑1 (MQ=255)

                                                                         gTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGaa                                                                            >  1:48649/1‑150 (MQ=255)

                                                                          ttGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATTATATTCATGTTCTGGAag                                                                           <  1:76736/150‑1 (MQ=255)

                                                                          ttGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAag                                                                           <  1:1734882/150‑1 (MQ=255)

                                                                              cGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGaa                                                                       >  1:1202423/1‑150 (MQ=255)

                                                                              cGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGaa                                                                       <  1:1900317/150‑1 (MQ=255)

                                                                                 ttATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACAt                                                                    >  1:645918/1‑150 (MQ=255)

                                                                                     cTGATTAAAAATCCTGAGTATTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACAGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAACAGAACATGgcg                                                                >  2:290042/1‑150 (MQ=255)

                                                                                      tGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGgcgc                                                               >  1:1720201/1‑150 (MQ=255)

                                                                                      tGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGgcgc                                                               >  1:1343882/1‑150 (MQ=255)

                                                                                       gATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGCACATGGCGCt                                                              >  2:144171/1‑150 (MQ=255)

                                                                                       gATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCt                                                              >  2:1058988/1‑150 (MQ=255)

                                                                                       gATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCt                                                              >  1:683407/1‑150 (MQ=255)

                                                                                            aaaaTCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGa                                                         >  2:26435/1‑150 (MQ=255)

                                                                                            aaaaTCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGa                                                         <  2:1564059/150‑1 (MQ=255)

                                                                                                tCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGa                                                     <  1:386864/150‑1 (MQ=255)

                                                                                                  cTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGATGAGAACATGGCGCTGTTGATGGAAc                                                   >  1:1463905/1‑150 (MQ=255)

                                                                                                  cTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAAc                                                   <  2:780914/150‑1 (MQ=255)

                                                                                                   tGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACa                                                  <  2:562072/150‑1 (MQ=255)

                                                                                                    gAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAg                                                 <  2:17918/150‑1 (MQ=255)

                                                                                                      gTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGc                                               <  2:594648/150‑1 (MQ=255)

                                                                                                      gTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGGGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGc                                               <  1:370511/150‑1 (MQ=255)

                                                                                                         ttttATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTAt                                            >  1:429267/1‑150 (MQ=255)

                                                                                                           ttATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCg                                          <  2:1494027/150‑1 (MQ=255)

                                                                                                                cgtcaggcgcgcgcAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAAc                                     <  1:1383278/144‑1 (MQ=255)

                                                                                                                cgaaaTGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAAc                                     <  1:1447907/147‑1 (MQ=255)

                                                                                                                cGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAAc                                     <  2:189366/150‑1 (MQ=255)

                                                                                                                       cgcgcgcAGTAGAAGCGATCCGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGcc                              <  1:1833089/150‑1 (MQ=255)

                                                                                                                       cgcgcgcAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCGTATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGcc                              >  1:1297658/1‑150 (MQ=255)

                                                                                                                       cgcgcgcAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGcc                              <  1:866839/150‑1 (MQ=255)

                                                                                                                        gcgcgcAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGATCAGGCTATCGCCAACGAAGGCCt                             >  2:143784/1‑150 (MQ=255)

                                                                                                                         cgcgcAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTg                            >  2:1325285/1‑150 (MQ=255)

                                                                                                                           cgcAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGtt                          >  2:1514557/1‑150 (MQ=255)

                                                                                                                            gcAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGGTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGtgt                         >  1:1085956/1‑148 (MQ=255)

                                                                                                                              aGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCTAGTGTCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTa                       <  2:774677/150‑1 (MQ=255)

                                                                                                                                  gAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGc                   >  1:1231192/1‑150 (MQ=255)

                                                                                                                                  gAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGc                   >  1:806266/1‑150 (MQ=255)

                                                                                                                                   aaGCGATACGTGTTCCGCATCCGGTACGCGTCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTTATCATATTCATGTTCTGGAAGAGAAAATGGCGTTGTTGTTGTAACAGGCTTTCGCCAACGAAGGCCTGTTTTATCGcc                  <  2:1794889/150‑1 (MQ=255)

                                                                                                                                   aaGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGcc                  >  1:925180/1‑150 (MQ=255)

                                                                                                                                   aaGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGcc                  <  1:1342132/150‑1 (MQ=255)

                                                                                                                                   aaGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGcc                  <  1:1448250/150‑1 (MQ=255)

                                                                                                                                    aGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCt                 <  1:570352/150‑1 (MQ=255)

                                                                                                                                      cGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTAc               >  2:859670/1‑150 (MQ=255)

                                                                                                                                       gATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTAct              >  1:1254853/1‑150 (MQ=255)

                                                                                                                                            cggtgCCGCGTACCGGATGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTAcc         >  1:1385799/3‑150 (MQ=255)

                                                                                                                                             tgtgCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGCACATGGCGCGGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCt        >  2:312387/1‑150 (MQ=255)

                                                                                                                                             tgtgCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCt        >  1:707070/1‑150 (MQ=255)

                                                                                                                                             ggtgCCGCGTACCGGATGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCt        >  1:146567/2‑150 (MQ=255)

                                                                                                                                               tgCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGc      <  1:519063/150‑1 (MQ=255)

                                                                                                                                                gCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCa     >  2:650019/1‑150 (MQ=255)

                                                                                                                                                gCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCAGCGCACGTAATCATATTCATGTTCTGGAAGAGAACTTGGCGGTGTTCATGGAACAGGCTATCGCCAACGCAGGGCTGTTTTATCGCATACTCTACCTGcc     >  2:1862371/1‑149 (MQ=255)

                                                                                                                                                gCCGCATCCGGTACGCGGCACCGTTTCGATGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCa     >  2:667897/1‑150 (MQ=255)

                                                                                                                                                   gCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTATGGAAGAGAACATTGCGCTGGTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATAGCCTACTCTACCTGCAgcg  >  2:1455030/1‑150 (MQ=255)

                                                                                                                                                   gCATCCGGTACGCGGCACCGGTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGATGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCCACGAAGGCCTGTGTTTTCGCCTACTCTACCTgcggcg  >  1:302472/1‑150 (MQ=255)

ACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCG  >  NC_000913/3995567‑3995863

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